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A method has been developed for studying the feasibility of 
reducing annoyance due to aircraft noise by modifying the flight 
trajectories over a community. Numerical optimization is used to com- 
pute the optimum flight paths, based upon a parametric form that 
implicitly includes some of the problem restrictions. The other con- 
straints are formulated as penalties in the cost function. Various 
aircraft on multiple trajectores (landing and takeoff) can be considered. 
The modular design employed allows for the substitution of alternate 
models of the population distribution, aircraft noise, flight paths, and 
annoyance, or for the addition of other features (e.g. , fuel consump- 
tion) in the cost function. A reduction in the reguired amount of 
searching over local minima has been achieved through use of the 
presence of statistical lateral dispersion in the flight paths. 

Three annoyance models have been studied: Noise Impact Index 

(Nil), Level-Weighted Population (LWP), and Highly Annoyed Population 
Number (HAPN). It is shown that Nil is not suitable for comparing 
different annoyance reduction strategies in any one community . Use of 
either LWP or HAPN as the criterion indicates that significant reductions 
in community annoyance are possible through this approach. 

The equivalent aircraft noise concept has been developed, by 
which the time needed for computation of the noise levels is reduced 
significantly when a large variety of aircraft is under consideration. 
This concept has also been shown to have applications in the problem of 
aircraft- trajectory assignments as it relates to noise reduction. 


It is found that, generally, the population distribution does not 
have a dominant influence in determining Uie optimum set of flight paths 
when the mixture of aircraft on those paths changes,* the optimum set 
may have to be recomputed for a change in the flight mix. 
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CHAPTER I 
INTRODUCTION 

The question to be considered is this: Is it possible to achieve a 

substantial reduction in community annoyance due to aircraft noise 
through the use of. trajectory modifications? While not a complete 
remedy for the problem of aircraft noise, such modifications will be . 
shown to aid significantly in its solution. What is particularly appealing 
about this technique of noise control is that its cost is relatively low 
compared with other methods. This study has produced a system that 
may be used to analyze the noise effect of aircraft operations on a 
given population and to minimize that effect (the community annoyance) 
to the extent possible via modifications in the flight paths. 

Historical Background 

The problem of severe aircraft noise began in the 1950's with the 
advent of the commercial jet airliners. Complaints and litigation pro- 
ceedings followed soon thereafter, with the result that numerous studies 
were begun. As the use of jet transportation became more extensive, 
so did the awareness of the effects of its noise upon the environment 
and people living near airport terminals. Today the aircraft noise issue 
is perceived by the airiine industry ab one of its most important prob- 
lems. 

Analyzing the noise effects of aircraft operations is not a wholly 
objective matter; the physical measurements of sound levels are only 
part of the analysis. Numerous physical measures of sound noise levels 
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have been defined, e.g. L^, Some of these attempt to include 

the frequency response of the human ear, while others try to indicate 
the degree of perceived annoyance. Ref. [1) is a useful compendium of 
these measures and rating schemes. 

Less easy to quantify, and open to many interpretations, are tlie 
subjective psychological and. sociological effects of the noise, collectively 
referred to as ”noise annoyance." Several attempts have been made to 
relate the physical measurement to the noise annoyance. Alexander [21 
examined some of the relationship models used in studies of the aircraft 
noise problem, Tarnopolsky [3] has studied the mental health effects 
from extended exposure to aircraft noise. While Borsky [4] showed 
that a linear relationship exists between tlie percentage of people 
"highly annoyed" and the noise exposure level in the range 55 to 

80 dB, Kryter [5, p.444] indicates tliat the relationship is nonlinear 
(power law) in terms of the maximum EPNL.^ More recent studies by 
von Gierke [6] and Goldstein [7] have used a single quantity to repre- 
sent the .average degree of noise annoyance to a community. This 
measure, the Noise Impact Index (Nil) has the properties of intensivity 
(includes the intensity of the annoyance) and extensivity (includes the 


^ Both of these are logarithmic measures of sound power intensity. L 
is the A -weighted sound level, with the weighting based upon thG 
frequency response of the human ear. is the average day-night 
sound level, based upon the average power density (A-weighted) dur- 
ing a 24-hour period. Nighttime sound levels are weighted more 
heavily in the average. Chapter III elaborates on these measures. 

2 The Effective Perceived Noise Level (EPNL) is a measure of the 
sound power density averaged over a short period of time (8 seconds), 
with spectral weighting, See Appendix A. 
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extent of the annoyance, i.e. the numbers of people affected). The 
Noise Impact Index is defined as the weighted sum of all the people 
affected, divided by the total affected population. Schultz [81 has 
shown that the proper weighting to be used for the Nil is a nonlinear 
function of the noise exposure level (in 

Methods for reducing * the annoyance due to aircraft operations 
generally lie in one of three categories: 1) noise source modifications, 

2) land use control, and 3) aircraft operational control. Noise source 
modification is more readily used in the design stage of future aircraft 
than it is in existing ones, where "retrofitting” the engine nacelles with 
sound absorption material is quite costly. For example, tlie estimated 
cost of retrofitting a Boeing 737 with sound absorption material is 
$202,000, while the cost of replacing the engines' front fans witlt 
quieter ones is $1,000,000 [9J. Land use control involves areas exper- 
iencing the most noise annoyance near an airport. Such land areas may 
be purchased by the airport, or perhaps zoned as industrial areas, 
thus reducing the numbers of people affected. This remedy is not 
suitable for metropolitan areas where land for housing is in short sup- 
ply. 

The third metliod for reducing noise annoyance, airGraft opera- 
tional control, requires a method of ttredicting noise levels that will be 
experienced as a result of whatever operational changes are made. 
Considerable effort has been made in this regard. Padula and Liu [10] 
examined a mathematical model which treated an aircraft as a moving 
prolate spheroid that scattered acoustic waves from a trailing point 
source. Barkhana and Cook [11] presented a model for determining the 
region on the ground that would experience noise above 70 EPNdB. 
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Dunn and Peart [12] studied five basic aircraft types; turbojet, turbo- 
fan, turboprop, V/STOL, and the helicopter, developing models of both 
the noise sources and the noise levels contours. These have been 
incorporated into compuiv/' programs for noise level estimation by 
Crowley et al. [13]. The statistical lateral dispersion of flight paths 
about the nominally specified trajectories has been analyzed by Filotas 
[14] , who finds that this has a significant effect on noise levels on the 
ground fairly near to the mean flight paths. The most sophisticated 
model devised for predicting sound levels is the FAA Integrated Noise 
Model. A variety of noise ratings may be estimated based upon actual 
atmospheric conditions and the mixture of aircraft landing and taking- 
off at a particular site. 

The more common methods used to reduce aircraft noise effects 
through operational controls are the procedures of two-segment 
approach and thrust cut-back on takeoff. Zalovcik [15] empirically 
determined that a combination of the two — a 6 degree slide slope 
(angle of ..the flight path with respect to the ground) followed by the 
standard 3 degree slope, and an associated thrust reduction -- could 
yield a decrease in sound pressure levels (SPL) of up to 14 dB. Click 
et al. [16], simulating aircraft operations near the San Jose Municipal 
Airport, studied similar operational modifications, including alternate 
routes for approaches and landings; however, no systematic method for 
varying the flight paths was employed. 

Jacob [17] presented an analytic solution to the problem of mini- 
mizing the average effective perceived noise level for landing aircraft 
confined to a vertical plane. Barkhana et al. [18] showed the 


feasibility of computing trajectories in three dimensions, with noise 

impact as part of the performance measure being minimized. Their 

♦ 

method, a combination of dynamic simulation and the use of steepest 
descents to determine the aircraft control inputs, required large 
amounts of computer time. A somewhat different metliod, parameter 
optimization, hai^ been shown by Rader and Hull (19] to be relevent to 
the computation of optimum aircraft trajectories. They employed a 
gradient method to find the angle-of-attack history for supersonic 
aircraft that results in the minimum time-to-climb trajectory at full 
power. A fifth-order Tschebycheff polynomial series was used to 
approximate the angle-of-attack history, with the polynomial coefficients 
as the parameters to be determined by the optimization scheme. 

Jacobson et al, [20] developed a system for determining optimum 
aircraft landing trajectories (via parametric optimization) with respect to 
the Nil, for Boeing 707 and 727 aircraft. In the work reported here, 
the system has been expanded to include both takeoff and landing 
paths, and a variety of commercial aircraft. 

Using this system, the current study has been directed to certain 
essential problems of noise impact reduction through the use of opera- 
tional controls. First is the matter of which has the dominant influence 
in determining the noise impact on a community : the particular assign- 

ment of aircraft to the various trajectories, or the population distribu- 
tion. This has an important bearing on the work of the air traffic 
controllors, who must make the flight path assignments; recomputing 
the optimum trajectories for a change in the aircraft mixture is, at 
present, a lengthy process. 



Also examined has been the problem of aircraft takeoff. Because 

they are less severely constrained, trajectories used in taking off are 
« 

more easily modified to reduce the noise impact than are landing trajec- 
tories. The endpoint of a takeoff path is not restricted to be a specific 
point in space, and the rate of climb is allowed to vary more than the 
rate of descent ih the landing case. 

Another problem inherent to trajectory modications is that of 
restricting the fuel consumed and/or the time of flight. Although 
reduced noise effects are desired, it would not be acceptable to achieve 
them at great expense in terms of fuel or time. A study has been 
conducted on the effects of including the fuel and time of flight in the 
performance measure. 

Previous work [21] relied on Nil as the measure of noise annoy- 
ance to a community . While Nil acts as a good indication of annoyance 
per person in comparing communities of different sizes, it does not 
reliably measure the change in annoyance for modified trajectories in 
any one community. This has prompted the examination of alternate 
measures and their effects on the optimum trajectories thus computed. 

The parameter optimization approach works well for airports of 
modest size, where the number of trajectories and hence, the number of 
parameters needed to specify them, is small. For larger airports, the 
number of parameters required can result in memory requirements that 
exceed the capabilities of the computer being used. An alternative 

scheme has been devised for determining the solution in this situation. 
In addition, the concept of "equivalent aircraft type" has been devel- 
oped, in which the noise characteristics of all the aircraft on a given 
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trajectory arc lumped together into a single "equivalent aircraft." This 
transformation has yielded a considerable reduction in the time required 
to calculate the noise levels experienced by the populace, and as a 
byproduct, some insight into the sensitivity of these levels to changes 
In the numbers of aircraft on the flight paths. 

Lastly, the* penalty function method of including the various con- 
straints in the problem formulation has been examined to determine the 
correct weighting for the penalties. 

Chapter II is an outline of the problem solution method. Following 
It are detailed descriptions in Chapters III to VII of the models devel- 
oped and employed in the optimization software. Chapter VIII contains 
the results, and Chapter IX, the conclusions and recommendations. 
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CHAPTER II 

THE APPROACH TO THE SOLUTION 


To examine the feasibility of reducing noise impact on a community 
using operational changes in aircraft flight paths, several major compo- 
nents of the system. must be considered. These include: specifying 

aircraft trajectories, incorporating passenger, pilot, aircraft, and regu- . 
latory restrictions/, and searching for "best" paths that reduce the 
noise impact. In this chapter, each of these areas is addressed in 
order to show the rationale behind the approach developed for solving 
the problem. 

Although specific mathematical models and functions have been 
employed in the development, it is not claimed that they are the sole 
choices. The modularity of the approach allows for perhaps more 
accurate replacements to be made in any of the areas subsequently 
discussed. 

Any attempt to reduce noise in a community requires some quanti- 

»>* 

tative evaluation (i.e. model) of the effects of the noise. From the 
model may be extracted a measure, a numerical function, that quantifies 
the noise effect.^ Many such models (and measures) have been estab- 
lished, emphasizing different noise effects. For example, the Speech 
Interference Level (SIL) attempts to quantify the interference of noise 
with speech communication , while the Percent Hearing Loss (PHL) 


^ It is preferred that the measure be a scalar quantity, a single value, 
for purposes of ranking different solutions to a particular noise 
problem. A unique method of ranking different collections of values 
does not exist, 


9 


measure is intended to predict that physiological effect for a given 
range of average noise levels. For the model of overall annoyance of a 
community because of noisCy there are two major considerations; the 
Intensity of the sound levels present and the exti-nt of those levels in 
terms of the numbers of people exposed to them. ’Among the measures 
that attempt to indicate both* the intensity and exter sity of the noise is 
the Noise Impact Index (Nil), which has been the principal i.iessure 
used in this study. 

The model of annoyance upon which Nil is based treats extensity 
and intensity in the following ways: For annoyance intensity, it is 

assumed the psychological response to noise depends upon the average 
noise levels experienced. Nighttime levels are weighted more heavily 
than are daytime levels. The intensity of annoyance is represented by 
a nonlinear weighting function of the average sound level. For aVmoy- 
ance extensity, the Nil model attempts to indicate the annoyance of the 
entire community exposed, by including the number of people exposed 
at different levels. Intensity and extensity are combined in the Nil by 
summing the number of people exposed to a given noise level, multiplied 
by the intensity weighting of that level. This sum is then normalized 
by the total population exposed to annoying levels, giving an indication 
of the average annoyance per person in* the community. 

Once the annoyance measure has been chosen, a practical way to 
evaluate it must be devised. In the use of Nil, this entails the use of 
two mathematical models: a noise level model and a population distribu- 

tion model. Under ideal conditions, it may be a'ssumed with reasonable 
accuracy that the noise levels from aircraft far away (several hundred 
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meters) depend only upon the distance from the aircraft. ^ Since the 
details of that dependence (and the flight paths) are known for all the 
aircraft operating in the area, it is possible to compute the average 
noise level at any point on the ground. The spatial distribution of 
people around an airport may be approximated by dividing the region 
into a number of smaller areas (cells). United States Census Data can 
then be used to obtain the number of people In each cell. For given 
sets of trajectories and the associated aircraft on them, the average 
sound level in each cell can be calculated. It is then a straightforward 
matter to calculate the Nil for those given conditions. 

With a practical implementation of the annoyance measure estab- 
lished, it is then feasible to begin the search for trajectories that 

result in the minimum possible value of the measure. Traditionally, 
aircraft trajectory modification has been performed using the calculus of 
variations or dynamic programming to solve the equations of motion of 
the aircraft, subject to certain constraint and minimize some criterion 
referred to as the cost function. Typically the cost function has a 
physical significance, such as the total fuel consumed by the aircraft 
flying these trajectories. While this method of dynamic simulation works 
well for aircraft confined to a vertical plane (see for example Erz- 
berger [22]), the Complexity of the equations in three dimensions 

results in a large increase in the amount of time needed for computa- 

tion . 


2 This also assumes given power and control settings on the aircraft. 




An alternative technique has been applied here. Each trajectory Is 
modelled as a specific functional in parametric form. The optimization 
process then consists of searching for the set of parameters that speci- 
fies the optimum set of trajectories, i.e., the set that results in the 
minimum value of the cost function. In this transformed version of the 
problem, the number of variables to be determined, and hence, the 
computation time, is smaller than in the original form, but the presence 
of various restrictions on the trajectories now creates some difficulty. 
Whereas these restrictions may be readily handled in the simulation 
technique and in fact are often automatically built into the problem 
formulation, they must be explicitly included when employing the 
parameter optimization scheme. 

Two classes of optimization methods exist for solving this type of 
problem. The first, referred to as primal methods, searches for the 
set of parameters that will result in the optimum solution, but the 
search is conducted over only those sets of parameters that specify 
acceptable*- solutions, i.e., no restriction is violated. The second class, 
approxiniating methods, searches over any sets of parameters without 
restriction. Inclusion of the problem constraints is accomplished by 
adding terms (penalty or barrier functions) to the cost function. The 
penalty functions are such that the cPst function increases when con- 
straints are violated; if no violation occurs, then no penalties are 
added. The unrestricted optimization will then, ideally, avoid 
restricted sets of parameters in its search for the lowest value of the 
cost function. For reasons given in Chapter V, an approximating 
method was chosen for this work. The unconstrained optimization 
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method used was the one of Davldon, Fletcher, and Powell, modified by 
Luenberger [233 to Improve its performance In a variety of problems. 

As for the formulation of the cost function, clearly the minimization 
of annoyance Is the true objective, given that the presence of the 
penalties is a necessity; however, other concerns may be included, 
weighted in prdportion to *thelr importance. For example, this work 
contains a study of the effects upon the solution if the total fuel con- 
sumed and time of flight are included. 

A suitable parametric form for mathematically describing the tra- 
jectories is one that requires a small number of parameters, but that 
still has sufficient flexibility to satisfy all of the constraints. A good 
choice for such a trajectory model is a set of linear segments joined by 
smooth arcs. The types of flight paths that result from this model are 
similar to those currently flown by commercial airliners. This is 
desirable from the perspective of the pilots, who must not be burdened 
with overly complicated flight paths. A trajectory is completely speci- 
fied once -the coordinate of the intersections of Its linear segments are 
determined, This is the result of the optimization process. 

Included in the constraints are; the dynamic capabilities of the 
aircraft, considerations of passenger comfort, maintenance of safe 
separation distances between two trajectories, avoidance of frequent 
high noise level flights over any one area, and a limit to the number of 
maneuvers required of the pilots. If some constraints place different 
requirements upon some property of the trajectory set, then the most 
restrictive of the constraints prevails. 


In the following chapters, the models and measures used in the 
solution method are treated in detail. Chapter III develops the annoy- 
ance model and considers alternative to Nil. Chapter IV derives the 
trajectory model. The optimization method Is described in Chapter V. 
Finally in this development, the cost function, including the formulation 
of the constraint* penalties , is discussed in Chapter VI. 


CHAPTER III 

THE ANNOYANCE MODELS 


As discussed in the previous chapter, a desirable community noise 
annoyance model includes both the intensity and extensity of the annoys 
ance and has associated with it a scalar measure that quantifies the 

m 

annoyance. Three measures are examined in the following sections; . 
Noise Impact Index (Nil), Level- Weighted Population (LWP), and Highly 
Annoyed Population Number (HAPN). First, however, models of air- 
craft noise level and population distribution, which arc required by the 
annoyance models, are developed, 


Aircraft Noise Model 

For most types of aircraft, the sound level or intensity (power per 
unit area) in dB units may be expressed as a function of the distance r 
from the plane; 

L Cj - c^logjQt (3-1) 


where c^ and Cg arr constants for a given type of aircraft operating 
with given power and control settings. 

The human ear converts sound energy into nerve impulses with 
some filtering of certain frequency ranges, In order to account for 
this, a measure of sound level that includes an approximation of the 
ear's frequency response has been devised (see Ref. {11). Known as 
the A-weighted level, L^ indicates the "Instantaneously" perceived ^ 


^ Actually, the time required to perceive a change in sound level is 
about 0.5 sec. See Kryter [5, p.3). 
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sound level. There are, however, other measures that incorporate the 
duration of the noise or the presence of high amplitude pure tones in 
predicting the judged sound level. Of particular interest here is a 
scale known as Effective Perceived Noise Level (EPNL), whose units are 
EPNdB. In developing the noise model for this work, data from the 
Federal Aviation* Administration's Integrated Noise Model (see Ref. [24]) 
has been employed. These relate the noise level in EPNdB to the 
distance from a particular type of aircraft. Figure 3.1 shows an 
example of the data for a DC-8-30 jet approaching for landing on a 3® 
glide slope. The solid line is the least squares fit of the data. It has 
the relation 

EPNdBjjj,.g_3Q = 155.74 - 28.32‘logjQj:. (3-2) 

Similar fits were obtained for other types of aircraft to be studied. 
The results are shown in Table 3.1. 

All three of the annoyance indicators to be examined use for 
calculating the annoyance intensity. The approximate conversion is 

: EPNL - 13dB. (3-3) 

This relation is derived in Appendix A, along with the procedure for 
calculation noise levels when multiple sources are present. 

The Population Model 

In order to assess the extensity of noise annoyance, a model of 
the population distribution is required . Some annoyance measures are 
cast in terms of population distribution as a function of noise level; 
however, for computational purposes, it is better to start with the 
spatial distribution of people. This is compatible with the noise model 


Bata from INM Bata Base [24 
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Figure 3.1 EPNL vs. Slant Range for DC-8-30 ,3° Approach 
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already developed (the two are appropriately combined later in this 
chapter) and has the advantage of being independently calculable (i.e. 
no noise calculations are required), 

As described in Chapter II, the distribution is approximaled by 
determining the number of people inside each cell of a grid that over- 
lays a map of the community, The grid should contain a sufficient 
number of cells to insure good resolution of the distribution, but not so 
many cells that large numbers of calculations occur in the annoyance 
measure evaluation. A scheme that satisfies both of these conditions is 
one in which the size of the cells grows with distance from the airport. 
The small cells near the airport provide better resolution in a region 
where the noise levels are changing more rapidly over the ground. 

As an example, Figure 3.2 shows the grid used for the Patrick 
Henry Airport at Hampton, Virginia. (Some of the grid lines in the 
center are not shown, so that the runways will be visible.) Here, the 
cells adjacent to the runways are 400 meters on a side. The length of 
a cell side grows both horizontally and vertically, according to the 
relation 

^k+l “ 

where Ij^, is the length of cell k in either direction. Covering a square 
region 32,000 meters on a side, the grid contains 576 cells, approxi- 
mately 11 times fewer than if all the cells were 400 meters on a side. 

Alternative grid schemes also may be used. Some examples are: 
concentric rings divided into sectors by radial lines (with the ring radii 
either constant, or growing with distance), and square cells of constant 
dimensions (if the number required is acceptable). 





riquro 3.2 F’opul.ition (»rid at Patrick Henry Airport U’artial) 
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Once the grid scheme is chosen, the number of people in each grid 
cell may be determined using United States Census Data in a computer 
program calied SITE II [25]. It requires as input a reference point in 
latitude and longitude, and coordinates with respect to the reference 
that specify the boundaries of each cell. The population in each cell is 
output, along with a variety* of demographic information. Although only 
the population data have been used in this study, the other data could 
be useful in an annoyance measure that includes social and economic 
information , 

In each grid call, it is presumed that the population is concen- 
trated at the geometric centroid of the cell. This is done to simplify 
the noise level calculations, i.e., all the people will receive the noise 
level present at the centroid. Chang [26, p.29] has shown that this is 
a good approximation in that the aircraft trajectories determined using 
this method are not significantly different from those calculated using 
more exact methods. 

-%» 

Level-Weighted Population (LWP) 

This indicator of community annoyance, devised by the National 
Academy of Sciences Committee on Hearing, Bioacoustics and Biomech- 
anics [27] assumes that the intensity of annoyance due to noise varies 
with the noise level. The relation between annoyance and noise level is 
characterized by the intensity weighting function W(L^j^) , where is 
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the average day-night noise level. is defined as^ 


^dn ” ^ 


H 

1 w^lO 
t=i 




10 


N 


^ j 1 for noise between 7 A.M. and 10 P.M. 
'^t 10' for noise between 10 P.M, and 7 A.M., 


( 3 - 5 ) 


N = Number of noise level samples taken in 24 hours, 
t-th sample of the noise level, A-weighted, 

A plot of is shown in Figure 3.3. The analytic expression 

for it is; 


. 0 . 103 * 1 .^ 

3.36 X 10“*’ -10 

" 0 . 2 . 10 “-«lcl„ + 1.43 X 


(3-6) 


This function is based upon a collection of social surveys of annoyance 
caused by various types of noise [8]. Although this synthesis respre- 
sents possibly the best available data for annoyance prediction, its 
matliematical model has the property of being unbounded. This is 
acceptable, though, if there is a bound to the largest value of 
likely to be encountered. 

Level-Weighted Population is defined by 


Lw = ; wa^„) d(i,,„) (3-7) 

where P(I^jjj^) is the number of people receiving noise between L^^^ and 
'^dn ^ *^^^dn^‘ The integration is performed over the range of 
considered to be annoying. By convention, the lower limit is 55 dB; 


2 relies upon the concept of energy addition, which is developed in 
Appendix A. 
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the upper limit is the highest level present. Evaluation of this integral 
is facilitated by approximating it as a descrete sum 

IWP (3-8) 

where i is the index of a population group, is the number of 

people in groupvi, is the noise level they receive and is 

the intensity weighting for that 


Noise Impact Index (Nil) 

This measure, also developed by the National Academy of Sciences, 
is defined as 


Nil = 


LWP 


(3-9) 


Where the symbols and limits of integration are the same as in Equation 
(3-7). Again, a discrete sum is usually used to approximate the integ- 
ral; 


= total population exposed to annoying noise levels. 
The purpose of the denominator is to normalize LWP , making Nil an 
indicator of the average annoyance per person in the community . Such 
an indicator is quite useful when comparing noise problems in communi- 
ties with widely different populations; (LWP is not a good choice for 
use in such comparisons because it tends to grow with population size.) 
There is, however, a serious problem with Nil regarding its use in 
comparing different strategies for noise reduction in any one community. 
As an example of the problem, consider a simple population distribution 
of just two clusters of people, as shown in Figure 3.4. The people 
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receive noise from aircraft flying on one of the tracks (ground projec- 
tlons) shown. Table 3.2 gives the results for each of the tracks. 

Tracks 2 , 3, and 4 all have lower bWP than Track 1, but the Nil 

values rank differently. Track 2 has the lowest Nil value and Track 
4, with the lowest LWP value, has the highest Nil. The problem occurs 
because both LW^* and the denominator in Equation (3-9) both change, 
but by different factors. Clearly, the use of Track 4, where 2,000 
people receive an of 67dB is superior (in terms of total community 

annoyance) to the use of Track 3, where 8,000 people receive an 

of 64dB. The LWP measure indicates this intuitive choice, while the Nil 
shows the opposite. In Chapter VIII, it will be demonstrated that a 
similar situation occurs for a real population distribution. 


Highly Annoyed Population Number (HAPN) 

Perhaps a more tangible indicator of community annoyance would be 
desirable, An example of such an indicator is HAPN, which attempts to 
measure the total number of people who are highly annoyed. 

The model of annoyance intensity used by HAP differs from that of 
LWP and Nil in two ways; the sound level scale used is maximum 
dB(A), and the intensity weighting function is:'"* 


„ ^ j 6.5885 . 10-'2 . (1^„)5-«7603 

6.5885 X i 75.3 dBCA) 


( 3 - 10 ) 


^ This weighting was obtained by a least square-error fit of data in 
Kryter [5, p.444]. 
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where 

f _ t maximum A-level received, 7 A.M. to 10 P.M* 

*^Amax “ ‘ maximum A-level received + lOdB, 10 P.M. to 7 A.M. 

%AP interpreted as the fraction of people receiving the noise who 

will be highly annoyed. 

The definition is 

I 

^'^HAP^^A,max,day^i * ^ '^HAP^^A,wax,night^i 

^Total (3-11) 

where the day-night notation is as indicated in Equation (3-5), Pj is 
the i-th population group, n is the number of groups, and P^otal 
the total population being considered. Note that daytime and nighttime 
levels are treated separately, so that it is possible for some population 
groups to be included twice. Dividing by the total population does not 
present the same problem here as does the denominator of Equation 
(3-5), since Pfj.Qj.ai is a constant for a given community. 






CHAPTER IV 

THE TRAJECTORY MODEL 

Having established in Chapter II that each trajectory should be 
determined by specifying the unknowns in a general parametric express 
Sion, the next matter is to choose that expression, using a form that is 
compatible with the restrictions Involved in the problem. In developing 
a suitable expression, or model, it is helpful first to examine the sorts 
of flight paths that are currently used in the region near an airport 
terminal (within 30,000 meters, where the noise levels are significant). 
Next, all of the restrictions that pertain directly to the trajectories 
must be inspected and cast in mathematical terms, Finally, a choice of 
the parametric form of the fUght paths may be made, and the restric-* 
tlons restated in terms specific to that form. 

Typical Flight Paths 
A. Landing 

'in general, an airport will have several runways and be 
receiving flights from a number of directions, For reasons of safety 
and organization (from the aspect of air traffic controllers), the number 
of patha available into the runways is limited py requiring that each 
aircraft first fly to a reporting point (also referred to as an arrival fix 
or an entry point), From there it proceeds to its assigned runway, 
with the requirement that as it passes over a point called the inner 
marker, its heading must be toward the runway and remain so until 
touchdown. The direction of all traffic, taking off and landing, on a 


given runway is fixed at any given time, determined chiefly by the 
wind velocity in the area. Figure 4.1 is a reproduction of an approach 
’’plate" used by pilots with instrument landing systems (ILS) in the 
aircraft. The projection of the trajectory onto the ground is referred 
to as the "ground track." This particuiar plate, for Patrick Henry 
International Airport at Hampton, Virginia, indicates two initial ap- 
proach fixes (lAF’s), Swing and Franklin, for use with the major run- 
way. The entry points and inner markers conveniently serve as end- 
points^ for the landing trajectories that are to be computed. 

B. Takeoff 

Takeoff paths are considerably less restricted than are land- 
ing paths. An aircraft taking off may execute a turn before crossing 
over the inner marker. It may then proceed along a path that is not 
required to pass through a reporting point (or "exit fix"). Figure 4.2 
shows a typical takeoff path, (ground track) for Lambert- St. Louis 
International Airport. 

..Except for their different endpoint conditions, landing and 
takeoff ground tracks are seen to be rather similar : each consists of a 

combination of straight segments and smooth turns. As for the vertical 
projection of each type of path (called the "profile"), landing aircraft 
usually do not descend with a glide slope of more than three degrees , 
while those taking off may climb at angles as great as fifteen degrees. 

i A fixed entry point is not required for aircraft with more sophisti- 
cated navigation equipment, such as VOR/DME (VHP Omni-Range/ 
Distance Measuring Equipment). The flight path model should there- 
fore allow for moveable entry points . 
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Restrictions 

The restrictions that pertain directly to the feasibility of the 
trajectories can be grouped into four categories: aircraft dynamics, 

passenger comfort, statutory regulations, and piloting considerations. 
A fifth constraint, concerning the maximum allowable number of high- 
level noise events for a given population grid cell, has been included in 
this work. 


A. Aircraft Dynamic Constraints 

The aircraft dynamic constraints limit a trajectory to one 
which given types of aircraft are capable of flying. The restrictions 
may be stated in terms of fixed properties of the aircraft and maximum 
allowable control surface deflections. In the event of several restric- 
tions on any one quantity, the most severe will prevail. The following 
expressions, developed from linear equations of aircraft motion [28], 
are divided into lateral (parallel to tlie locally flat ground) and longi- 
tudinal (in a vertical plane) constraints . 

4. Lateral Constraint 


where 
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The right-hand side of Eqn. (4-1) will be a constant for a given 
type of aircraft, and the left side may be written in terms of the func- 
tion that specifies the trajectory. In Figure 4.3a, the ground track is 
described in Cartesian coordinates as y = f(x). With respect to the 
arbitrary x-axis, the slope of the flight path is dy/dx, with corres- 
ponding angle * 

f = tan"^dy/dx (4-2) 

however the quantity of interest is the rate at which that angle is 
changing (also the yaw rate of the aircraft): 


V = gf tan’^ dy/dx = gfCtan'^ g)g (4-3) 

But dx/dt is the x-component of the aircraft's velocity (with respect to 
the ground-fixed coordinate system): 


V = ~ - V cox V = V( 

X dt ' 


[1 + (dy/dx) ^1"^ 


-) 


(4-4) 


where V is the average, total velocity, so that 


= = J- (4-5) 

[1 t (dy/dx) R(x) 

where R(x) is the radius of curvature. The constraint (4-1) may now 
be written 

d2y/dx2 (cj^ + C2C3)min(drj,6r2,6r2) 

2. Longitudinal Constraint 

Yd 1 Y 1 Y(, ' (4-7) 

max max 

where y is the angle of the flight path with respect to the horizontal, 

and Yj and Vp si's the maximum allowable angles of descent and 
^ax ^ max 
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climb, respectively. Referring to Figure (4-3b), it is seen that, in 
cartesian coordinates 

tan Y = ^ (4-8) 


so that Eqn, (4-7) may be rewritten 


dz 

Vd 5 di - Vc • 

« max max 


(4-9) 


B. Passenger Ride Quality Constraints 

The passenger comfort, or ride quality, constraints limit the 
operation of an aircraft to those maneuvers found statistically to be 
acceptable to the passengers. An expression of these limitations (see 
Ref. [301) for an aircraft executing a coordinated turn® with climb or 
descent is 


n - f 0.76 for 90% passenger satisfaction 
^ ‘ 1.76 for 80% passenger satisfaction, 


where “ maximum roll angle (degrees) 

‘Pmax “ maximum roll rate (degree/sec) 

’^max “ maximum pitch angle (degrees) 

^max “ maximum flight path angle (degrees). 

For a coordinated turn and considering only the roll angle 
part of Eqn. (4-10), the constraint may be written 


R 




(«-n) 


? A coordinated turn is one in which no sideslip occurs, l.e. the 
component of gravitational force along a wing provides the centripetal 
force necessary to maintain the particular turning radius being used. 
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or, using derivatives of the trajectory function: 

[1 t 


d^y/dx^ 




(A-12) 


where 


r0.19 for 90% passenger satisfaction 
‘0.50 for 80% passenger satisfaction 


g 

<!> 


5! average total velocity 

to 

= acceleration due to gravity 
= roll angle. 

C. Regulatory Constraints 

Regulations affecting allowable flight paths may pertain either 
to safety, governmental security, or noise annoyance. Safety restric- 
tions involve single trajectories (e.g. avoidance of low altitude obsta- 
cles) and multiple paths (e.g. maintenance of safe separating dis- 
tances); military installations often prohibit general and commercial 
aviation flights from passing too closely for security reasons; and there 
may be certain areas (e.g. schools and hospitals) which are particularly 
sensitive to noise from nearby aircraft. 

’Two safety restrictions that apply to every airport have been 
included in this study . They are: 

1. A safe distance should be maintained between any two 
trajectories as far as is possible. Because several 
trajectories may share a single runway, the region in 
which this restriction applies often excludes an area 
within several miles of the airport. The minimum allowed 
separating distance used in this work was 800 meters . 


2. No landing aircraft should be required to execute a turn 
of more than 45 degrees, as it crosses the inner marker, 
in order to come into alignment with the runway. This 
is referred to here as the final heading requirement. 

Restrictions that are related to special areas to be avoided, 
for security or Annoyance reasons, are particular to the given airport 
community being considered. Although not included in this study, 
these constraints could be formulated and contained in the problem as 
described in Chapter VI. 

D, Pilot Operating Constraints 

It is possible that the optimum trajectories computed would 
require more maneuvering of an aircraft than a pilot can or is willing to 
perform, especially in the near terminal area, where his attention is 
required on a number of matters. The constraints then, are to limit 
the number and complexity of the turning maneuvers on the flight 
paths . 

E. threshold Noise Constraint 

A difficulty with ail three of the annoyance measures dis- 
cussed in Chapter III (and other ratings) is that an "improvement,” as 
determined by these measures, can be made by exchanging a large 
number of people exposed to medium noise levels for a small group 
exposed to extremely high levels. To avoid this possibility, a con- 
straint is imposed that prohibits any area* ** (grid cell) from receiving 


* For some cells near an airport, the high levels may be unavoidable 
due to runway locations . These cells must be excluded from the 

restriction. 
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maximum noise levels (L^ during a flyover, at or above a '’thres- 
hold" level, more than times per day, The threshold level and 

value of chosen should be as small as possible, without overly re- 
stricting the possible trajectories. 

Flight Path Model 

m 

In selecting an appropriate parametric form for specifying a tra- 
jectory, it is necessary to weigh the benefits of a large number of 
parameters —more degrees of freedom and hence, more flexibility — 
against the attendant increase in computation needed to determine the 
values of those parameters. Chang [26], for example, reported that 
using only five terms in a Fourier series representation of aircraft 
landing trajectories (the unknown parameters were the Fourier coef- 
ficients) led to overly complicated flight paths, while Increasing the 
number of terms would have resulted in unacceptable computational 
requirements without necessarily simplifying the paths. 

The model chosen for a flight path in this work avoids that prob- 
lem by directly relating the maximum number of turning maneuvers to 
the number of parameters to be determined. Each trajectory is repre- 
sented by a chain of linear segments joined by helical arcs (circular in 
the ground plane, linear in the profile). Figures 4.4a and 4.4b show 
an example of a path with three segments. The two intersections of the 
linear segments in Figure 4,4a are called "corner points;" the cartesian 
coordinate of these points, along With those of the other trajectories, 
make up the set of parameters to be determined. In the ground plane, 
a circular arc with minimum radius (determined by the dynamic and 
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Figure 4.4a Ground Track, Linear Segment Representation 



Figure 4.4b Profile, Linear Segment Representation 
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passenger comfort constraints) is inscribed in the angle between any 
two linear segments. The arc, tangent to both segments, describes the 
ground track of the path flown during a turn at that corner. Between 
tangent points, as shown in Figure 4.4b, the profile of the path is 
approximated by a linear segment. The resulting turn, in three di* 
mansions, is a helical arc. • Since each corner point (requiring three 
coordinates) may have at most one turn, the complexity of the flight 
path is limited in a natural way by the allowed flexibility, i.e. the 
number of corner points used. 

Also treated in an implicit way are the endpoint restrictions. 
Since only the corner points' coordinates are variables, the endpoints of 
a trajectory may be placed in the appropriate (constant) positions 
without need of an explicit restriction to. hold them constant. Por 
example, in Figure 4.4a, a landing.^ path ijegins "at the entry point 
(arrival fix) and ends at the inner marker, 'after which the aircraft 
flies straight to the runway. FoCr.tskcof(8* (amd iandings .not required 
to report to the entry point), the endpoint* away from the runway may 

• * • • . e • 

• •• * V** • • •• 

be made variable as shown in Figure 4.5. Here; a pseudo-segment has 


been added so that the true endpoint becomes a variable corner point. 

The pseudo-segment must be long^ enough to* prcvjent a violation al the . 
dynamic and passenger comfort c ong^ aihf s Jj^subse^uentT y^'^ ^ ) , 

• • • • •• • # s • • # 

and no noi.se calculations should inVoWiP this ■^ftUousu'iifegmertf.t' . 

explicit statements of the remaning constraints i are ^ Wfuirisi lor 
this particular trajectory model. The lateral dynmttc' restriction, 'liqn. 

t «.v -SM V. ^ • - • : 

(4-6) and passenger comfort constt^t (4-^2)jnrerboth ‘itw^f^ 

.-F.X A ^ *7 tif -A- ^ , 

ties involving the radius of curvat^n^etof •i8e'groimd7^c^l^'*tJiiwTOT 

combined to yield 

“rr-'A 


1 


Corner Point #3 



Figure 4.5 Variable Endpoint 
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. r V2 , V } , (4-13) 

R _ max { (C J+C2C3 )mitv (6r j , ) 

where R is the radius of curvature of any point on the ground track. 

For the linear segments of a path, R -> »; for the turns at the corner 

point, R is just the radius of the circular arc used to describe a turn 

in the ground track. Generally, the passenger comfort term (first on 

the right in Eqn. (4-13)) is the larger, so that a constant value of R 

(Rmin) that satisfies Eqn. (4-13) may be employed for all turns once 

the maximum value of 'Q is known for a trajectory. 

Incorporation of the longitudinal dynamic constraint is equally 

straightforward. Eqn. (4-9) involves only the vertical slope (the slope 

of the profile with respect to the ground plane), a quantity that is 


readily calculated once the coordinates of the tangent points are known 
(see Figures 4.4a,b). Appendix C gives a derivation of the formulae 
used to obtain these coordinates. 

Because of the form of the flight path model used here, some 
additional '^restrictions on the feasibility of any trajectory must be in- 
cluded. As shown in Figure 4.6, it is possible to specify corner points 
that result in nonsensical trajectories for a given value of Rj^^jj^- 
Figure 4.7 indicates the pertinent quantities to use in the restrictions. 
The first restriction prevents a path’ from becoming disconnected be- 
tween any two corner points by requiring that 


R . 


R . 

min 

4. 

mxn 

tan^ 

r 

tan^ 


(4-14) 




OortUM’ Point //3 



Gonun’ Point GoviuMr Point Gornor Point 

n H 


Figure 4.6 Corner Points that Specify an Impossible Trajectory 



Figure 4.7 Rciations Needed to Specify Feasible Flight Paths 
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Second, for the outermost segment in the chain, an explicit requirement 
that the circular arc have a tangent point on the segment^ must be 
stated, i.e.. 



^ a - ‘r ' ' 

tanj 

k 

Maintaining a safe disfance between any two trajectories may be 
accomplished by sampling the distance between their ground tracks at a 
number of locations, as shown in Figure 4.8a. This test makes two 
assumptions j that the profiles of the two paths are .sufficiently similar 
to allow these sampled distances in ^ ground plane to be good approx- 
imations of the actual distances, and that enough samples can be taken 
to obtain an accurate representation of spacing between the two paths. 
The first assumption is not an object of concern because it leads, at 
most, to an underestimation of the spacing. As for the number of 
samples needed, it has been found in this work that as few as ten 
measurements will suffice; however, as shown in Figure 4.8b, the 

successive*' values y. must be tested not only for minimum magnitude, 

but also for a change in sign, indicating a crossing of the paths. The 
two constraints are then; 

n i '■ml,.- t = 1. s M-16) 

^ i = 1, s-1 (4-17) 

where s is the number of samples. As stated previously, when two or 

more flight paths share the same runway, or use runways whose ends 


®The assumption of tangency is necessary for determining the 
expression for the arc, as developed in Appendix C. On the other 
segments, Eqn. (4-14) or (4-18), (4-19) assure this condition. 


< > 


Figure 4.8a Sampling the Distances Between Two Paths 



Figure 4.8b Two Crossing Trajectories with y: > d . 

i ~ min 





are close together, the separation constraint, Eqn. (4-16), should not 
be applied near these runways where the close spacing between trajec- 
tories is unavoidable. Furtlier, at larger airports with sophisticated 
navigational and air traffic conti'ol systems, both the separation and 
crossover restrictions, Eqns. (4-16) and (4-17) may not be required. 

The final Beading constraint is modelled by requiring tliat the 
angle 6, and length 1^^, shown in Figure 4*9, satisfy the following: 

6 > 135®, (4-18) 

Ij < ~1. (4-19) 

Together, these guarantee that a landing aircraft will have a path 
aligned witli the runway before reaching it, and that any final turn will 
be no greater than 45®. 

Finally, the threshold noise constraint (which is independent of 
the tr'ajectory model) requires an analytic form* In each grid cell to 
which the constraint applies, tlie maximum A-level (%^nax^ exper- 
ienced during each aircraft flyover may not exceed a threshold level 

(L^ ThreslP ^max expressed as 


where 


i=N>V 
N -t s 


. - ^A, Thresh^ ® 

index of the first flyover after flights with L 
- Thresh 


(4-20) 
A , max* 


N| = number of flights during a 24 hour period, 

] = index for the grid cells to be included in this constraint, 
i.e., a set denoted J. 





CHAPTER V 

THE OPTIMIZATION SCHEME 

In Chapter II, the notion of representing each aircraft trajectory 
in parametric form was introduced, the idea being to reduce the amount 
of computation required for determining the optimum set of flight paths, 
This chapter develops the numerical algorithm that will be used to . 
determing the optimum parameter values. In fact, several algorithms 
are described in the development, giving a good representation of 
available methods that are compatible with the solution approach adopted 
in this work. 

All of the optimization methods to be discussed require a criterion, 
referred to here as the cost function, by which to Judge the improve- 
ment being made in a system. The cost function must depend, perhaps 
implicitly, only upon the parameters needed to describe the system (in 
this case, those that specify the trajectories). As formulated in Chap- 
ter VI, the cost will be a function (explicitly) of the noise annoyance, 
the solution restrictions (constraints), and if desired, other attributes 
of the system, All of these components, however, depend upon the 
trajectories, which in turn, are functions of the parameters. 

Most optimization algorithms are iterative descent methods. They 
approximate the solution to a problem usually in several steps, with 
each successive step making use of previous approximations. This 
search progresses in such a way that the value of the cost function at 
each iteration is lower the previous value; the search is then said to be 
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descending.^ All of the optimization methods suitable for this work 
involve derivatives of the cost function with respect to the unknown 
parameters. Because the cost function includes the annoyance model, 
which lacks a tractable analytic expression , derivatives of the cost are 
best determined numerically. This strongly influences the choice of 
which class of optimization methods to employ. 

Of the two suitable for use in this work, primal and approximating 
methods, the former relies more upon extremely accurate derivative 
calculations. These methods search over only those parameter combina- 
tions not prohibited by the restrictions. Altliough primal algorithms 
offer the advantage of "suboptimal” solutions at each iteration^ (reduced 
cost and all constraints satisfied), the necessity of large numbers of 
calculations to compute the derivatives accurately outbalances the at- 
tractive features. 

The approximating methods begin with a problem in the form 

minimize f(x) 

X6E" (5-1) 

subject toj h (x) < c 

where x is an n- vector whose components are the parameters to be 
determined f(x) is Ute cost function, h(x) is an m- vector- valued func- 
tion that includes all the constraints, c is a constant m-vector, and E*' 


^ It is always possible to cast an optimization problem in a form that 
seeks a minimum value of the cost. 

^ If for some reason the iterative process must stop before reaching the 
optimal solution , these sub-optimal solutions offer at least an improve- 
ment over the initial condition. 
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is Euclidean n-space. With the appropriate transformation, the problem 
becomes unconstrained : ^ 


minimise ?(x) = f(x) t w'*^$(x) (5-2) 

xeE*' 

where w is a m- vector of constants, called the weighting or controlling 


parameter, and $(x) is an m-vector-valued function that depends upon 
the problem constraints, 

Two forms are popular for use in constructing $(x): barriers and 

penalties. Barrier functions are not useful in the work for two rea- 
sons. First is their requirement that the initial point in the search be 
a feasible point (satisfy all the constraints). Determining a feasible 
starting point can be quite difficult for cases with many trajectories 
having several segments. More important though, is a fundamental 
defflciency of barrier methods. Theoretically, the term w ^(x) would 
create a 'barrier" in the n-dimensional surface defined by the modified 
cost function in Eqn. (5-2). Presumably, the solution point Xj^ 


(obtained in the k-th iteration) would not move near the barrier, in- 
curring a high cost, in the next iteration. This, however, is not the 
case if large changes in x are allowed; it is possible for the solution to 
"jump" the barrier, giving an that has a lower cost, but that is 
infeasible, Without explicitly checking the solution at each iteration, 
this condition would go "undetected" by the optimization algorithm. 


2 Formally, the problem is constrained by the requirement xeE*^, but 
computationally , the solution is unrestricted in real n-space.’” 
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Fortunately, the penalty methods are somewhat more robust in this 
respect. Although they also permit infeasible solutions to pass for 
feasible ones, the condition is detectable and automatically corrected by 
the standard procedure for implementing the penalty methods. 

A popular form for the penalty function is (see Ref, [23]) 

' Pj(x) = 

with the weighting parameters chosen as 

= [Mi» M2> •••» Mi > 0- (5"^) 

The problem is then expressed as (c.f. Eqn. (4-2)); 

m 

minimize £(x) + Z p.P.(x). (5-5) 

XCE" " 

To remove totally the approximation in the constraints, introduced 
by the penalties, the problem must be solved a number of times, with 
increasing values of the Pj. A sequence with |i| > |i., is used to 

determine a sequence of solutions {xh * While some of the x^ may 
violate constraints, it may be shown (see Luenberger (23, pp. 278- 
280]) tliat any limit point** of the sequence {x^} associated with the 
sequence is a solution to the original problem, Eqn. (5-1). 

As a practical matter, though, it is necessary only to use a value 

of that results (at worst) in acceptably small violations of the con- 

« 

straints. Some experimentation with the particular problem being solved 
may be required to determine the components of p that will give accept- 
able results. Too large values tend to slow the convergence to the 


If these exists a subset K of the positive integers such that the sub- 

k k 

sequence { x ]|^ ^ converges to x, then x is a limit point of {x } . 
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sofutlon, while small values may result in constraint violations that are 
excessive (by whatever approximation standards apply), 

With the contrained problem formulated as an unconstrained one, a 
suitable optimization method must be selected, The following sections 
develop several traditional algorithms, leading to the technique chosen 
for this study. * ’ 

As with most nonlinear problems, the task of finding the minimum 
of a nonlinear function f(x) lacks a general, closed-form solution. 
Iterative technique, however, can be used to make successive approxi- 
mations of the solution, until the calculated value is judged to be rea- 
sonably close to the actual one, 

Once a starting point is chosen, new points are generated by a 
recursion of the form 


Xk+l = 2ik "k 4 


(5-6) 


where dj^ is the "direction'* vector and oij^ is a scalar that determines 
how far to move along the direction dj^ in order to reach x^+i’ 
values of and d|^ vary in the methods subsequently discussed. 


Steepest Descent 

Perhaps the most "short-sighted" of its type, this method often 
gives satisfactory results, but more important, it serves as a foundation 
for the more sophisticated algorithms. 

Given a function f(x) and a point the task is to find a good 
direction (in n-space) in which to search for a new point such that 
f(Xj^) < f(5^). The first approach is conveniently provided by a prop- 
erty of the gradient operator 
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(a row vector), 


(5-7) 


namely that £(x) ® y^f(x) Is a vector that points in the direction in 
which f(x) is increasing most rapidly. Since the problem is to find the 
minimum value of f(x), the direction in which to search is along -g(x), 
CChis will result ’in the most rapid decrease of f(x) for that iteration . 
(Later methods sacrifice this immediate improvement for other advant- 
ages , ) 

As stated earlier, It is not practical to obtain an analytic expres- 
sion for g( x) in the problem under consideration. An approximate 
numerical form 




( 


f(Xj, , Xj^ ^ Ax, > > . , x^) “ f(?l) 


Ax 


(5-8) 


i == l,n 

has been employed. In the limit Ax 0, Eqn. (5-8) yields the exact 
gradient; small values of Ax should therefore be used to obtain accurate 
approximations. 

The method of steepest descent is stated as the following algo- 
rithm: 


Steepest Descent Algorithm 
Step 0: Choose an initial point set k=0, 

Step 1: Compute dj^ = -aj^_(x^) 

Step 2: Minimize (with respect to 
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Figure 5.1 Steepest Descent Example 
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Iterative values of up to k»4 are shown as circles, connected 
by line segments to indicate the progression. This example demon- 
strates a major problem of steepest descent, namely, its slow conver- 
gence properties. Its initial step to is directly "downhill,” causing 
a large decrease in the value of f(x); however, once in the shallow 
region, it oscillates ineffectually, making slow progress. The curvature 
along the valley prevents the search direction, from pointing 

directly to the solution. A more rigorous explanation of the conver- 
gence properties of this algorithm may be found in Luenberger [23, pp. 
148-X55J , 

Newton's Method 

A more sophisticated approach to the general problem of minimizing 
a function f is to approximate that function with a quadratic form in the 
region of the solution: 

fCx) z ^ '*■ *4 2,^92,* (5-10) 

This approximate function has a minimum at 

'* it _i 

2 “ 9 E (5-11) 

provided that g is a positive definite symmetric matrix,® Taylor's 
Theorem® may be used to expand f(x) about x, keeping only up to 
quadratic terms: 

f(x) 2 : ■** SfCXkXx-x^) + %(x-Xk)^ ' (5-12) 

where £(Xj^) the Hessian of f, is the rvxn matrix of second-partial 
derivatives of f(x), evaluated at x^, 

® In R*', a matrix g, is positive definite if x"^(S2 > ^ 

® See Pierre [321. 
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ra2f(x)' 

i^J ■ 

Comparing t^rms of (5-12) with (5-11) gives a minimum’ at 

(r\) - -lF(x^)l"^ — 
from which follows the recursion 

' 4+1 = 5k ■ 

The algorithm for Newton's method is then: 

Newton's Method 


.f(Xk)^ 


(5-13) 


(5-14) 


(5-15) 


Step 0: 

Choose an initial point Xq, set 

Step 1: 

Compute fi(Xj^), [F(Xj^)]“'. 

Step 2: 

= X|, - (F(x^)l'‘j(x^). 


Set k=k+l 


Go to Step 1. 


As before, this procedure continues until a stopping criterion is met. 

No line search is required, but computing, inverting, and storing 
the inverse Hessian requires more computing resource, than does the 
steepest descent algorithm. Newton's method does, however, converge 
rapidly when is near the solution x*. This depends somewhat upon 
how nearly quadratic f(x) is; a scaling factor and a line search 
may be included for cases where the quadratic approximation is not 
sufficient. The recursion then becomes 


-1 


£k+l = ik ■ “k<?<Sk)‘ 4<Sk>' 


(5-16) 


’ The requirement that Q (and F) be positive definite is based upon 
the necessary and sufficient conditions for x* to be a minimum point. 
(See Luenberger [23, pp. 110-114]). 
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and a line search to minimize f(Xj^ + «]<^[E(Xk)r^aCXj^)) must be added 
to Stop 2 ol; the algorithm, 

Newton's Method: Example 2 

Figure 5.2 shows the results of applying Newton's method to the 
same problem as* in Example*!. The pure form of the method, without 
the search parameter was used with excellent results. The seem- 
ingly odd behavior exhibited in iteration 1 demonstrates how this 
method, in sacrificing the initial rapid progress of steepest descent, 
converges swiftly on the solution after only five iterations. 

Conjugate Gradient Methods 

Of the descent techniques examined thus far, one converges quite 
slowly when near the solution, and the other, quite rapidly, but at the 
expense of large increases in computation, A good compromise is pro- 
vided by the conjugate gradient methods. These are developed quali- 
tatively here; theoretical treatments may be found in Ref. [23] and 
[32], 

Conjugate gradient techniques are a subset of the more general 
conjugate direction method, As the name implies, gradient information 
is used in computing the search directions; it is the special conjugacy 
of these directions that distinguishes the methods. The quadratic 
problem considered previously serves to illustrate the notions involved. 

A quadratic function 

(5-17) 
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rigure 5.2 Newton's Method Example 
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has a minimum at 


X* = 2 


(5-18) 


if Q is positive definite. This solution may be expanded, using any 


n-1 


linearly independent set of vectors as 


n-l 

X* s 2 a . d . . 
" i=0 


(5-19) 


If the dj are Q-orthogonal (or Q-conjugate) , i.e. djQd. = 0 for 
i 7 ^ j, then the following results®: 


^-l. 


n-l 


q^b = 1 ad. , 
i=0 ^ 

-1 

2Q = I a Qd 
i=0 


m ^ n — 1 n . 

dTOO" b = I a.d!qd. 

- i_o 


T T 

d!b = a.dtQd. , 
-J- J-J^J ’ 


(5-20) 

(5-21) 

(5-22) 

(5-23) 


a 


d.b 

^ dTqd . 


The solution, in terms of the conjugate directions, is then 


T 

n- 1 d ;b 

X* = 1 -7ir — d. 

.T« , -1 


i =0 dtqd. 

-I 


(5-24) 


(5-25) 


® Proof that the g-conjugate set {d.} are linearly independent: 

Assume that the d. are not linearly independent. Then there exist 
n-l 

a. ^ 0 such that 2 a.d. = 0. Multiplying by d.Q yields (for each i) 

1 1 1 - J 

a.d.Qd. = 0, Since Q is positive definite, d.Qd. > 0 and a- = 0, 
contradicting the assumption. 




61 


All that remains is to specify the dj in terms of Q and b, the known 
quantities in the problem. 

What the conjugate gradient algorithm does is to compute the d. 
iteratively, using linear combinations of the gradients in all previous 
iterations . 

k ’ Conjugate Gradient Algorithm 


Step 0: 

Set k = 0, Xj^ = Xq, Uq = -Sq = b - QXq. 


Step 1: 

T., 

"Sk 4 

Compute a. = ■ 

(5-26) 





^^k+i = 4 ^ 44 

(5-27) 

Step 2: 

fik+iQ^k 
Pk “ T 

(5-28) 


444 



^ktl ' "^ktl ''' ^k-k 

(5-29) 


(Sk = 2Sic - i) 

(5-30) 


Step 3: Set k = k+1 and go to Step 1 


These iterative calculations are exact for a purely quadratic problem 

and the solution, Eqn. (5-25) is obtained in n steps. 

In general, the function being minimized is not quadratic. Near a 

minimum point, the Taylor expansion is 

f(x) = f(x*) + ^ t (x-x*) + 0(lx-x‘>''*l®) . 

(5-31) 

A necessary condition for x* to be an unconstrained minimum point is 
that 2^f(x*) = 0. The fouth term becomes negligible as x x*, leaving 
only the quadratic term in x. This leads to the following extension of 
the conjugate gradient algorithm in which Q is replaced by F. 
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General Conjugate Gradient Algorithm 


Step 0: 

Set k = Oj ^ “ ’gQ == SZ. ^(^0^ 

(5-31) 

Step 1: 

T. 

■&k 4 

Compute a. - — 

(5-32) 


Sk+1 = iSk * “k4 

(5-33) 

s£ep 2: 

Compate , 2^l£(xk+l> 

(5-34) 

If k ran-1 (m an integec ) , 



. o _ ^k+l£^4^4 

compute Pt = — ST-— — 

ii£(iSk>4 

(5-35) 


s=t t 

(5-36) 


k - k+1 



Else , 



4+1 " "Sk+1 

(5-37) 


k = k+1 



Go to Step 1. 

In Step 2, the "else" branch is called a spacer (or restart) step. 
It' is inserted to aid in the convergence to tlie solution by reinitiating 
the conjugate direction calculations every n steps. 

A disadvantage of the above algorithm is its requirement of calcu- 
lating F(x). Even if storage of the nxn matric is not an obstacle, 
numerical accuracy (along with that of still a concern. Errors 

may grow during the iterative construction of the dj^ to the point that 
even the spacer step cannot "recover" the search and no progress is 
made. 


I 


6 


Figure 5.3 
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Step 2: 

Set = -S^g^. 

(5-40) 

Step 3: 

Minimize f(xj^ + ®*k~k^ with respect to a > 0 

(5-41) 

Set 

^k+l " Xk ^ “k^k 

(5-42) 


^k “ “k-k 

(5-43) 


\ ~ fik-tl " h 

(5-44) 

(Note: cij^ must be selected so that > 0) 


Step 4: 

1£ k ^ mn (m an integer), set 


n _ f n -k^k^k-k , i^k^k . ^k^k 

SiC+1 Silr T * ^ 

alh\ 

(5-45) 


Set k = k+1. 



Go to Step 2. 


Else, 

Go to Step 1. 

A, Requirement for step-wise descent 

The above algorithm, along with steepest descent and New- 
ton's method, belongs to a class of methods that use the recursion 

-k+1 “ -k " ^kMk (5-A6) 

to compute the minimum of f(x), where a > 0, and is an n x n 

matrix that modifies the direction -gj^. For steepest descent, V j[, 
for Newton's Method, Mj^ = IF(xj^)]’^, and for this method, Mj^ = 

point which remains to be proved). A desirable 
property of such iterative methods is that each step produces a lower 
value of f(x). The requirement on Mj^ that insures this is now derived, 


Using Taylor's theorem to expand f(x) 


about x^ gives: 

- - k+1 






(5-47) 
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into which Eqn. (5-46) may be inserted to yield: 

f(x^^j) = f(x^) - + 0(lal2). (5-48) 

For small a, the trailing term becomes negligible compared to the second 
term on the right. If is to be less than f(Xj^) for small Ou > 0, 

then 


* ’ * ^ (5-49) 

This can be insured for all f 0 by requiring the Mj^ be positive 

definite. (Clearly, for steepest descent this is the case, but for New- 

-1 

ton’s method, = [F(Xj^)l is not guaranteed to be positive definite 
except near a minimum. ) 

B. Positive Definiteness of Sj^ 

Given the algorithm for generating the Sj^, it will be shown 
that the positive definiteness of follows from that of Sj^. For x s 


_T„ .IV .. 

T 


S .§.|c+i- ” Sk- ■ 


k- -k- 


'k T„ 


where 


T 

V - ^k^k 
" Jo „ 
3k§kaic 

Substituting a = b = §'^5^ 9ives 


” “k+l- ~ ''^k ^ 




b^b 


1 + 


T 

^k^k 


(5-50) 


(5-51) 


By definition, so that 

T T T 

lik^k ~ ^kSk+1 ■ Hk^k • 


(5-52) 
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Since is the minimum of f(x) along = “ic^k' 


M(x + „d ) = VkiAc ^ 0 

iW sK 


Therefore Eqn. (5-52) yields 


T _ T 
~ ”^k^k 


By definition, £|^ ~ “k-k " "‘*k -k^k' ®° 

eK = = “ksfekSk' 

Eqn, (5-51) may now be written 


X S. . ,x = Yu 1 


(a'^a)Ck'^k) - 


(5-53) 


(5-54) 


(5-55) 


(5-56) 


By the Cauchy-Schwarz inequality,® the first term on the right of Eqn. 
(5-56) is positive semi-definite, as is the second term, by virtue of the 
definition of inner product. To demonstrate that both terms cannot 
vanish simultaneously, assume that a = |b, for which case the first 


term vanishes. Then x = 4gj^ aud 


EkS ' fek\ = 4“k2k§k% * " 


(5-57) 


sincoi Sj^ is positive definite. Therefore, x^Sj^^j^x > 0 for all x 0, 
Since the algorithm starts with Sq ~ I, all of the Sj^ will be positive 
definite. Q.E.D. 

An extremely useful piece of information in this proof regards 

T 

the determination of ot-^. The magnitude of in Eqn. (5-55) is not 

paramount; only its sign affects the definiteness of For this 


® Cauchy-Schwarz Inequality. 
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reason, the algorithm requires that the line search in Step 3 be accur- 
ate enough only to give > 0. (When applied to nonquadratic 

» 

problems, the requirement that minimizes f(x) along “jtdj^ 

becomes more stringent), 

.1 

C, Convergence of Sj^ to F 

For the modified ’ Da vidon-Pletcher-PoweU method, it may be 
shown that If f(x) is quadratic (F constant and positive definite), then 


E^jEj =0, 0 < i< j < k 

(5-58) 

§kn£Ei = Ei. 0 < 1 < E. 

(5-59) 


from which the convergence of the algorithm in a finite number of steps 

-1 

will be proved, along with the fact that = F . 

T 

A Taylor expansion of £(x) = _ f(x) gives 

fiCx) • (5-60) 

Using Eqn. (5-42) - (5-44) in the above yields 

Sk = Sk+i - Kk = lik+i - esk = £Ek' 

From (5-45) 

" «k+l^k " % (5-62) 

The proof of (5-58) and (5-59) is by induction: They are 

true for k = 0 by virtue of (5-61) and (5-62), Assuming that they are 
true for iteration k-1, it will be shown that they are true for iteration 
k. 

Using (5-61), may be constructed from previous itera- 
tions : 

Sk = * • • • - Ek-i) (5-63) 
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From (5-53) and (5-58) 

» 0 < i < k . (5-64) 

Transposing (5-62) nnd substituting into (5-64) yields 

ElE§kSk " ° ^ ^ 

From the definition of the above may be rewritten 

‘ ■ eIIEk = 0 I i (5-66) 

proving (5-58) for.k. 

Still assuming that (5-61) is true for k-1, 

§kEEi = Ei . 0 < i < k. (5-67) 

which combined with (5-62) gives 

• ' «k+l^i ~ % * 6 1^1 k* Q.E.D. 

Since the are F-orthogonal, they are linearly independent 
(see footnote 8)*, If the matrix P whose columns are the pj^ is con- 
structed 

' £ “ (EqI £ 2 } '*'1 Efi-i^ > (5-68) 

then (5-59) yields 


- F P - P . (5-69) 

The linear independence of the Pj^ insures that P exists, so that 

§n £ = PP'^ = I (5-70) 

and 

I = F"^ . (5-71) 

Finally , this method is seen to be a conjugate direction meth- 
od, by virtue of the F-conjugacy of the pj^, For the quadratic case 

used in the preceeding analyses, it will converge in n steps or less. 

-1 

The presence of the scaling factor alters the eigenvalues in Sj^ ~ F 
at each iteration such that they span a small range. This leads to a 




better rate of convergence than was generally obtained In the original 
method. For a development of this feature, see Ref. {23J. A pure 
steepest descent step every n steps insures that the conjugate direction 
process will be initiated in the vicinity of the minimum, where its con- 
vergence is fastest. 

P 

Modified D-F-P Method; Example 4 

Returning to the problem (5-9) used in the previous examples, the 
last algorithm is employed with excellent results, Figure 5.4 shows that 
after five iterations, the approximate solution is comparable to that of 
the conjugate gradient method in Example 3. 

The last algorithm was chosen for use in this study for 
several reasons. First, it has convergence properties superior to those 
of steepest descent. Second, it requires only first order derivative 
information, which may be generated practically. Finally, the search 
directions employed are guaranteed to decrease the cost function at 
each iteration, an assurance that the conjugate direction methods lack 
for nonquadratic problems. 

Until now, the problem of nonconvex (or multimodal) cost 
functions has not been considered . The preceeding search methods all 

f 

refer to finding *'a minimum'' of the cost, but in nonconvex cases, there 
is more than one minimum. The lowest valued of these so-called rela- 
tive or local minima is referred to as the absolute or global minimum. 
There is no way to guarantee that any of the standard iterative tech- 
niques will converge to the global minimum of a nonconvex function. 
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For simple functions, different starting points may be used to find all 
of the local minima and hence, the global one. In Chapter VII, it will 
be seen that the cost function employed in this work has far too many 
local minima for this technique to be practical. A scheme for examining 
ranges of the variables will be shown to be applicable. 


CHAPTER VI 
THE COST FUNCTION 


With analytic expressions for the noise annoyance (Chapter III) 
and the system constraints (Chapter IV), it is possible to formulate a 
mathematical criterion, the cost function, by which the amount of any 

It 

change may be judged. As explained in the preceeding chapter, in- . 
eluding the constraints as penalties in the cost function allows the use 
of standard, unconstrained optimization techniques, which are more 
practical for this work than are the so-called primal methods. While the 
annoyance is the primary quantity to be minimized (with the penalties 
as a necessary adjunct of the cost), it may be desired to include other 
properties of the system in the cost, such as the fuel and time of flight 
required for each trajectory. 

The Penalty Functions 

All of the constraints developed in Chapter IV, Eqns. (4-9) and 
(4-13) through (4-20) are in the form 

1 ( 6 - 1 ) 
as required by Eqn. (5-1), even though the unknown x does not ap- 
pear explicitly in any of them. Penalties in the form of Eqn. (5-3) may 
then be written for each constraint as ’follows : 

[Dynamic and Passenger Comfort Constraint, Eqns. (4-9), (4-13)] 


P, = Z Z (max(0,max{ ^ — > ■ j. v 
1 traj. corner ^2® 

points 


]-Rl) 


(6-2) 
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P, = X X («ax[0, tany. - - tany ])' 

^ traj. segments dl dX 

(Geometric Feasibility Constraints, Eqns. (4-14), (4-15)] 


^3 = ^ 

traj. corner 


Z (niax(0, 
points 


^min 


^min 

tan| 

T 

tan| 


- 12])^ 


= 1 (max(0, 


R 


min 


traj . tans 


a 


(6-3) 


(6-4) 

(6-5) 


[Separation and Crossover Constraints, Eqn. (4-16), (4-17)] 


P5 = 2; 1 (raax[0, d . - y. ])' 

^ all traj. i=l ^ 

pairs 

^ 2 (max[0, - 

® all traj. i=l ^ 

pairs 


( 6 - 6 ) 


(6-7) 


(Final Heading Constraint, Eqns. (4-18), (4-19)] 

= Z (max[0, 135° - 6]) 2 
traj . 

Pg =: Z (max( 0 , 1 - |l])^ 

traj . 

(Threshold Constraint, Eqn. (4-20)] 


P 


9 


Z (max( 0 , { Z max( 0 ,L. . 
jeJ i=N* A,ina>i,i 


^A,Thresh^^^^ 


(6-8) 

(6-9) 


( 6 - 10 ) 


The Total Cost 

Adjoining the penalties to the desired annoyance measures , say 
Nil, gives the total cost function: 
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^(S) - ( 6 - 11 ) 

where A(x) is the annoyance (Nil), is the annoyance weight, |j is 
the m- vector of penalty weights, and P(x) is the m- vector of penalties 
(there are m constraints). As before, the variable x is present only in 
a formal sense; it does not appear explicitly in these terms in this 
work. ‘ ’ 

They annoyance weight w^ governs the relative importance of the ' 
annoyance, keeping it from being dominated by the penalty terms. 
This is most important in cases where the optimal solution is on a 
constraint boundary, In such a region, |j P(x) will in general be 
nonzero, but as discussed in Chapter V, small values are acceptable. 
The purpose of w^ then is to insure that whatever reduction can be 
made in the total cost corresponds to a decrease in the annoyance, 
rather than a decrease in penalties that represent tolerable constraint 
violations. When the constraints are violated excessively, the attendant 
large penalties naturally should comprise the dominant portion of the 
cost. 

Experimentally, the following values of w^ and the have been 
found to give satisfactory results when Nil is used for the annoyance 


measure : 


= 104 . 106 

(la’teral dynamic, passenger comfort) 

M2 = 102 . 106 

(longitudinal dynamic) 

P3 = 104 . 106 

(geometric feasibility) 

= 1 Q 4 - 10® 

(geometric feasibility, first segment) 
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Mj = 10“ 

(separation) 

Me = 10“ 

(crossover) 

= lo’’' - 10“ 

(final heading) 

llg = lO'* - 10“ 

(final heading) 

Mq ~ 

(threshold noise). 


The rather wide range in (Jg'nierits some explanation. It has been often 
found ^ in this work that the optimum trajectories will have the maximum 
vertical slopes allowed by the longitudinal dynamic constraint, Eqn. 
(4-9) ♦ This is intuitively expected, since the maximum slopes give the 
highest average altitude, the largest average distances from all the 
populace, and hence, the lowest average noise levels. Because of this 
tendency in the solutions, the associated penalty (P 2 ) generally is non- 
zero. Starting with a small value of 1 J 2 allows the optimization to work 
on the ground tracks without hinderance from penalties on the profiles 
(vertical components). After several iterations, ^2 1^® increased if 

excessive longitudinal constraint violations are occurring. 

Modifying the Cost Function; Inclusion of Fuel and Time of Flight 

An example of altering the cost function to include other attributes 
is a study conducted to examine the relationship between the noise 


1 In cases where a trajectory segment passes over a region of low 
population density, increasing its slope does not significantly lower 
the noise levels and hence, the total community annoyance. The 
optimization stops when the change in cost between successive itera- 
tions is less than some criterion. Therefore such a segment may not 
be exactly optimum, but practically it is so. Fuel consumption rates 
that vary with climb and descent rates are not included in this analy- 
sis. 
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annoyance, and the fuel and time of flight needed. It is assumed that 
the latter two vary linearly with the flight path length only. That is, 
constant power settings and noise propogation characteristics are used. 
The resulting term to be added to the cost is 


T,T0F 


'^FJOF^T 


( 6 - 12 ) 


where Wp .pQp is the fuel and time of flight weight, and l.p is the 
trajectory length. Eqn. (6-11) becomes 


fj(x) = w^A(x) + ^f,TOF^-^ (6-13) 

Using f]^(x) as the cost function, a relation between the annoyance and 

Ik 

fp .pQp was obtained by varying Wp ^Qp. This investigation is re- 
ported in Chapter VIII . 


CHAPTER VII 

INTEGRATION OF THE MODELS 


The annoyance, flight path, and cost function models are combined 
with the chosen optimization method to create a system for analyzing the 
aircraft noise annoyance problem and the feasibility of remedying it 

i , 

through trajectory modifications. While the components of the system . 
are essentially modular (and hence, capable of being altered individually 
without replacing the entire system), some linkage among them is re- 
quired. Further, some special considerations must be given to larger 
airports and to some population distributions. 

The Flyover/Noise Simulation; Small Airports 

At any airport, there are, in general, a number of flight paths 
(both for takeoffs and landings) in use on a given day, with different 
types and numbers of aircraft using these trajectories, In this work, it 
is assumed that for the period of time under consideration, the type of 
each path., (takeoff or landing) and the runway to which it is assigned 
do not change. In addition, the noise level data used is that for 
landing aircraft with constant power and control settings.^ 

Calculation of the annoyance measure requires information about 
the distribution of the sound levels • (L^^^) on the ground; for the 
current population distribution model, only the sound levels at the cell 


^ These assumptions are not fundamental to the operation of the sys- 
tem. A more sophisticated noise model, for instance, that employs 
thrust dependent noise data could be used instead of the present 
one. 
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centroids are needed. First, the position of each type of aircraft, on 
each trajectory, is sampled at equal time intervals (10 sec.). Then 
Eqn. (3-1); along with data from Table 3.1, is used to compute the 
sound level (L^) at each centroid due to each type of aircraft , from 
every sample position and trajectory. This computational approach is 
practical so long as * the number of flight paths and aircraft types is 
small. Referring to Appendix A, these levels for all of the aircraft are 
added on an energy basis, ^ and the average day-night level at 

each centroid may then be calculated from Eqn. (3-5). 

Knowledge of the distribution, along with the population 
distribution, allows the calculation of either LWP (c.f. Eqn. (3-8)) or 
Nil (c.f. Eqn, (3-9)); however, evaluation of HAPN is not necessarily 
as straightforward. This measure requires that the sound level (L^) 
be calculated at each centroid, at each sampling instant, and that the 
maximum values (day and night) be retained. At smaller airports, 
where the flight operations are dispersed and non-overlapping in time, 
the values depend upon only one aircraft at any sampling instant. 
For a larger number of operations where the overlap (in time) is un- 
avoidable, the sequence of flight operations must be explicitly included 
in the calculations, since at every sampling instant the types of 


2 It is not necessary to compute (at each centroid) for every sam- 
pling instant, since these levels '%ill themselves be averaged (using 
energy addition) to get Only the energy term, Eqn, (A-4), 

must be computed and accumulated (correctly weighted, c.f. Eqn. 
(3-5)) into the total centroid energy for every sampling instant. If n 
aircraft of the same type fly on a give trajectory, then their total 
energy contribution at any centroid is n times the energy added by 
one aircraft. 




aircraft flying will determine the value of at any centroid, In this 
study, only a situation of the former type has been considered for the 
purpose of comparing HAPN with LWP and Nil. 

Large Airports and the Equivalent Aircraft Concept 

Summing (on an, energy basis) over all trajectories, sample points, 
and aircraft types may become impractical as the number of any of 
these grows, depending upon the computing resources available. At 
large airports, such as the Lambert-St. Louis facility in St. Louis, 
Missouri, the number of trajectories and aircraft (fifteen paths and five 
hundred operations per day) precludes the straightforward approach of 
the last section. In such a case, the summation over aircraft types 
may be removed by replacing the mixture of planes on each trajectory 
with a single equivalent aircraft , whose sound level distribution is the 
same as if all the aircraft assigned to that path were flying together on 
it. For this treatment to be valid, it is necessary that either Nil or 
LWP be used as the annoyance measure. The reason is that both of 
them are functions of the distribution, which depends only on the 
total sound energy distribution, integrated over time (one 24-hour 
period) and not on the time history of the energy contributions (except 
for the distinction between day and night occurrences). 

The exact expression for the total sound level (using any power 
intensity related scale) contributed from n sources is (see Appendix A) 

n L./IO 

I»Y0T - 10 • .5, 






(A-6) 


with, from Eqn. (3-1) 


‘■i = 'i,i ' 'a, I • 

where the subscript I is the aircraft index and n is the total number of 
aircraft. A simple, approximate expression for Eqn. (A-6) may be 
obtained in the following way. For all the aircraft on a given traject- 

f 

ory, the appropriate Lj's are added (on an energy basis, and using the 

Ldn convention which weights nighttime energy terms with a factor of 

♦ 

ten), i.e. an equivalent level L may be written 
* n L.(r)/10 

L^(r) = 10 • logjQ ^ (7-2) 


where the subscript k denotes that this is the L* contribution from the 
k-th trajectory, and w. is the same time-of-day factor as in Eqn. (3-5). 
If Lj^(r) vs. r is plotted for a variety of values of r, an excellent 
approximation of the resulting data is the expression 




‘^l,k,eq " 


logjQC*;) 


(7-3) 


where c^ and C 2 gq are the equivalent aircraft noise coefficients 
for the k-th trajectory. An example of this is shown in Figure 7,1. 
The largest value of r was chosen to be 60,000 meters, since that is 
the largest value of the range that would be encountered in the near 
2 one area (32,000 meter radius) of an airport. In this example, the 
estimated correlation coefficient squared in the least square-error fit 
(dotted line) was r^ = 0.9995. (Spacing the abscissa values such that 
*'j+l “ ®*'j' ^ ^ constant, results in their being equally spaced 

on a logarithmic scale, and gives a better curve fit.) 
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With the coefficients Cj^ and C2 for all the flight paths, 
the total L* at any point on the ground, at each sampling instant is 
Just the sum (on an energy basis): 

* N it Jio 

t * 10 • log., r 10 . (7-4) 

t 10 j 

where N,p is the number of trajectories, and t is the sampling index. 

Since the time-of-day weighting has already been included in the 

expression for becomes (c,f. En. (3-5)) 

N t*/io 
z 10 ^ 

'■dn = '» • 1“*10 

where N is the number of noise level samples. 

This method of using equivalent aircraft coefficients reduces the 
number of calculations needed to determine at each grid cell by a 
factor approximately equal to the number of aircraft types that are 
operating In an area. Since at many airports this number is frequently 
on the order of ten, a considerable savings in computation time may be 
achieved. (Minimal computation time is required to determine the equi- 
valent aircraft noise coefficients.) 

As it is the method relies upon the assumption that Cj^ . and C2 | 
are constant for each aircraft type. If varying thrust levels and 
control settings are to be considered, it might be possible to make 
these coefficients functions of time, and for a given airport, predict 
their values (and hence, those of the equivalent aircraft coefficients) 
by estimating the time histories of the thrust and controls. A second 
disadvantage of the method is that values are not computed at each 
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sampling instant. Thus, the threshold noise constraint either must be 
omitted or altered so that the limitation Is In terms of In each grid 
cell, rather' than the maximum A-level experienced, Finally, the errors 
Introduced In the values by the equivalent aircraft method are 
small (< 1 dB); however, the nonlinear intensity weighting function, 
W(L^j^), may amplify these, errors, giving Nil and LWP values that 
differ by several percent from those calculated exactly. The amount of 
error in Nil and LWP changes with the sets of trajectories, but the 
locations of the local minima do not differ significantly in comparisons of 
the exact and approximate calculations. ' 

Equivalent Aircraft Coefficients and the Assignment Problem 

An interesting feature of the equivalent aircraft method is that the 
noise level coefficient (k is the trajectory index) has a strong 

dependence only on the total number of aircraft on a particular flight 
path, as long as there is a uniform mixture of aircraft types present. 
A plot of k eq vs. the L^j^-equivalent number^ of aircraft per tra- 
jectory (for the Lambert/St. Louis International Airport) is shown In 
Figure 7.2. The dotted line is the least square-error fit and has the 
form 


^ For example, using the Phoenix Sky Harbor Airport case, the error 
in Nil for the nominal (currently specified) trajectories is 2.1% and 
for the optimum set, 0,2%. The corner pts, (equivalent aircraft 
method) differ from the optimum ones (exact calculations) by no more 
than 100 meters . 

2 The Lj convention of weighting a nighttime sound level with a factor 
ten is “equivalent to adding the sound levels , on an energy basis, of 
ten identical sound level occurences. Hence, one nighttime flight is 
equivalent (in L^ calculations) to ten daytime flights. 
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*^l,k,eq " « + b . logjo , 


(7-6) 


where N^^^'is the L^^^-equivalent number of aircraft per trajectory and 
in this case, a = 144.65 and b = 10.16. This fit has an estimated 
correlation coefficient squared, r^ = 0.993. As for the coefficient 
k eq' Lambert-St, Louis case) are not as 

strongly correlated with Figure 7.3 shows a log-log plot of 

^2 k eq ^a/c' square-error fit (dotted line) for 


> 250. The functional form is 


‘^2,k,eq " ^”a/c » 


(7-7) 


where a = 23.44 and b = 7.89 x lO"*^; the estimated correlation coeffi- 
cient squared is r^ = 0.894. Restricting to values greater than 
250 allows for a good fit of the data with a simple functional form. For 
Na/c < 250, Eqn. (7-7) will not accurately predict C2 g^, and an 
individual fit of the data , as in Figure 7.1, will be necessary . 

These relations between the coefficients and Ng/g# Eqns. (7-6) and 
(7-7), have a bearing on the problem of how to assign aircraft to the 
different trajectories, with the community annoyance from noise as the 
criterion . 

In this work, the assignment problem is not considered (the 
assignments are assumed to be known and remain constant for a given 
airport); however, it will be a necessary addition in subsequent refine- 
ments of this study. The aforementioned relations may be useful in 
formulating a performance criterion to be used in the initial assignment 
of aircraft, or in altering an existing plan. 


c - 
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The Integrated System 

A general summary of the systems operation is provided by the 

t 

flowchart in Figure 7.4. 

The input required is: 1) the population distribution data, 2) the 
number of each aircraft type and their time-of-day designations (day 
or night), 3) the noise level characteristics (either by individual air- 
craft type or equivalent aircraft noise coefficients), 4) the constraint 
data (minimum turning radius, threshold noise level, etc., and penalty 
weights), 5) an initial set of flight paths (specified by their end points 
and corner points), and 6) the stopping criteria (the maximum number 
of iterations, and the lower limit on relative changes in the cost func- 
tion). 

For the initial trajectories, the total cost must be evaluated, before 
starting the iterative search. (If the system is being used only for 
evaluation of the trajectories, then no search is performed, and the 
process stops after this step.) First, the annoyance measure (e.g., 
Nil, LWP, HAPN) is evaluated. Then the constraints are tested: a 

violated constraint results in a positive penalty; otherwise, the penalty 
is zero. Finally, the total cost is computed by adding the penalties to 
the annoyance. As shown in Chapter VI, other cost terms may be 
added to the total; the results of one Such modification appear in Chap- 
ter VIII. 

Now the optimization loop begins. After each search step, a test 
is made to determine if any of the stopping criterion is satisfied. If 
not, the iteration continues; if so, the process stops and the "minimum" 
annoyance trajectories are determined. In fact, the set determined may 
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Figure 7.4 System Flowchart 
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be only a locally minimum set (a problem discussed in the next section), 
or perhaps not even that (the maximum limit on the number iterations 
allowed may have been reached); however, because of the type of 
optimization method used (Chapter V), each successive set of paths will 
result in a lower total cost (and a lower annoyance, if the penalties are 
zero). * ’ 

Besides a limit on the number of iterations, the stopping criteria 
include the tests; 1) Is the gradient of the cost function zero? (i.e., 
the search has found a local minimum), and 2) Is the relative difference 
between successive values of the cost negligible? (i.e., the search has 
come close enough® to a local minimum to warrant no further search- 
ing). 


Local Minima and the Population Distribution 

A problem with the optimization, already expressed in Chapter V, 
is that of local minima in the cost function. As an example of the 
problem, consider a cost function of two variables, for which the con- 
tours are plotted in Figure 7.5. Points A, B, and C are the local 
minima; point B corresponds to the lowest value of the cost, and hence, 
is the global minimum. The difficulty associated with this situation is 
that of finding point B by use of any of the optimization metliods dis- 
cussed in Chatper V. Given an initial point from Which to search , each 


® The error in the solution, E - x - , where x^ is the k-th 

iterative solution and x* is the exact solution, may in fact be large. 
If is small, though, and the relative difference between suc- 

cessive costs is small, then further searching may not lead to a 
significant reduction in the cost. Lacking other information, it may 
be assumed that, practically, a local minimum has been located. 
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of those metliods will converge to one of the local minima, but will not 
automatically find the global minimum. Which local minimum is found 
depends upon the starting point in the search. The conventional 
technique used to find the global minimum is to perform the optimization 
a number of times, from different initial points. In cases where the 
cost function has only a few minima, this approach may indeed allow all 
of them to be found. If the total number of local minima is unknown, 
but judged to be small, it is often possible to estimate the probability 
of overlooking a number of minima. That probability may then be made 
acceptably small by repeating the optimization witl) a large enough 
number of starting points. Ref. [35] and [36] are good examples of 
this approach, 

The cost function used in this work has an extremely large number 
of local minima. This is demonstrated by considering the population 
distribution shown in Figure 7.6a along with the simple straight-path 
trajectory. This path is at a constant altitude (300 m.). There are 
1000 people at each dot on the map. In Figure 7.6b, the resulting Nil 
values versus 0 are plotted (a typical mixture of commercial aircraft is 
assumed for generating the noise levels). 

Even though Nil is computed from noise levels experienced by the 
entire community, the local minima exhibited in Figure 7,6b demonstrate 
the influence of high levels that occur locally within the community. 
For instance, at = 20“ , the trajectory passes close to one population 
center and almost directly over another; the resulting noise level at 
these points, and hence, their contributions to the Nil will be high: a 
peak occurs in Figure 7.6b at Oj^ = 20®, Conversely, at O2 “ 26®, the 






Figure 7.6a Uniform Discrete Population Distribution and 
Straight Flight Paths 





Figure 7.6b Nil vs. Trajectory Direction 



population points nearest to the trajectory are farther away from it than 
those in the first case. The noise levels they experience are lower, 
and a minimum in the Nil occurs at ©2 = 26®. (At 0 = 45®, the sharp 
peak is a local effect caused by the large number of population centers 
directly beneath the path.) 

This example demonstrates two important problems. First, there 
are many local minima in the cost function, caused mostly by the use of 
a particular representation of the population distribution.** In a realis- 
tic case, there are more population centers and trajectories, which lead 
to an even larger number of local minima. Second, the arbitrary lower 
limit on values to consider (55 dB) for use in the Nil has altered 
the results that would be expected intuitively. The envelope of Nil 
peaks in Figure 7.6b is seen to decrease as 6 increases, i.e., as more 
of the population receives higher noise intensities, the community annoy- 
ance decreases. This occurs because the denominator in Eqn. (3-9) 
increases, in this case, faster than LWP. Once again, the inapprop- 
riateness .of Nil for use in comparing the annoyance from different 
trajectories in a given community is in evidence. 

A dramatic change in the number of local minima in this example 
occurs if Nil is calculated without a lower limit on the interval 
(c.f. Eqn. (3-9)) but with the total population fixed. This modified 
version of Nil will be denoted Nil'. Figure 7.6c shows the results 


** In each population grid cell, the people are assumed to be concen- 
trated at the geometric centroid of the cell. The overall result is a 
discontinuous (discrete) density distribution. 



Figure 7.6c Nil' vs. Trajectory Direction 
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(solid line), The complete disappearance of all the local minima is 
fortuitous in this case (the influence of local effects still persists, 
though, as seen In the inset); it occurs as a consequence of the uni- 
form (discrete) distribution of people and also because of the particular 
degree of freedom given to trajectory variations here.® 

There remains the dilemma of local minima, caused by macroscopic 
variations in the population density. Again, referring to Figure 7,6a, 
assume that in Region I, there are 20,000 people at each population 
center; in Region II, 25,000 people at each center; in Region III, 
10,000 at each center; all other centers have 1,000 people. The result- 
ing values of Nil' vs. 0 appear in Figure 7.6c (dotted line), At 0 = 16® 
and 0-70®, there are local minima, caused by the lack of large num- 
bers of people close to these paths. Because the truly large clusters 
of people are in Regions I, II and III, it is reasonable to conclude that 
these minima are real effects of the population distribution, and not 
artifices of its approximate representation. Note that at 0 = 8°, 82°, 
42®, 48®, „ there are slight dips in Nil', caused by the flight path's 
avoidance of all the nearest population centers, 

The possibility of artificial minima remains, though. Note that at 
0 = 8®, 42®, 48® and 82®, there are inflections in Nil'. Although not 
minimum points, there correspond closely to minima in Figure 7.6b. 
Further, no such behavior is seen in Figure 7.6c (solid line) for the 


® Only 0 may vary , which limits the possible paths. In the next sec- 
tion, it will be seen that trajectories with several segments still have 
local minimum problems with this type of population distribution 
model. 
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uniform distribution. It may be suspected then that artificially pro- 
duced minima may occur even with the use of Nil', depending upon the 
actual distribution of people and the use of this population distribution 
model . 

Global Searching 

Regardless of the cause of the local minima, some procedure for 
determining the optimum solution must be developed. First, it must be 
recognized that finding the exact optimal set of flight paths is not as 
important as determining a set (if one exists) that produces a signifi- 
cant reduction in the annoyance measure, compared with the set pre- 
sently used (nominal set), A given set will be considered optimum if 
there is not another set, significantly different, with a lower cost 
function value, 

The definition of "significantly different" is developed in terms of 
the statistical lateral dispersion of individual flights from the specified 
trajectories. Galloway [37] states that the lateral dispersion of landing 
aircraft is approximately 750 m. at a distance of 16,000 m. from the 
airport. This corresponds to an angular dispersion of approximately 
3“. For takeoffs, where the paths are less tightly controlled, the 
lateral dispersion is given in terms of the standard deviation (a) of the 

f 

side-line distance (y) from the specified path. For straight-out take- 
offs, the standard deviation grows linearly with the distance 1 and is 
given by 

Oy 0.08 1 . (7-8) 

After a turn, the dispersion Is expressed as 

= 3.5 X 10^ ajz , 


(7-9) 
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where qJz (the standard deviation in altitude dispersion divided by the 
average altitude during the turn) is generally a constant; the value 
a_/2 ~ 0.2/ for variable weight transports, will be used here, 

Eqn. (7-8) may be used to estimate the angular dispersion on 
straight segments of a takeoff by assuming that all of the flights will lie 
within 20y of the specified paths (a Gaussian distribution would contain 
0Z% of the paths within 20y of the mean path). The outermost path, to 
either side of a straight segment, will then form an angle, Q^jigp/ with, 
the segment, given by 

®disp = = S'”' 

A set of trajectories will be considered "significantly different" 
from another set only if any path lies outside the dispersion limits of 
the corresponding path in the other set. As an example of how this 
may be used to restrict the amount of searching required in the optimi- 
zation process, the flight paths at the Lambert-St. Louis International 
Airport will be considered, Figure 7.7a shows the twelve nominal 
takeoff paths and three nominal landing paths. The shaded regions 
represent the predicted dispersion about each path. Note that the 
trajectories with turns possess more average dispersion, c.f. Eqn. 
(7-9); this is due to the extra drifting during a turn. 

Searching for the optimum trajectory set within the dispersion 
limits of this set is unprofitable since such an optimum would not by 
definition be significantly different from the nominal set. Of course, it 
may happen that the nominal set is the optimum, but regardless, any 
searching must occur "outside" the disperson boundaries. Different 
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starting points (representing significaivtly different trajectory sets) are 
used to perform the optimization a number of times. The lowest valued 
local minimum is then taken to be the global minimum. 

Choices for the starting point are limited by two factors: the 

lateral dispersion on these paths and the constraints (particularly the 
aircraft dynamic*' and passenger comfort restnctions). For example as 
shown in Figure 7.7b, there are only five choices for track #7 starting 
points. Moving track #7 farther to the east would cause it to overlap 
track #8, in which case the two paths would be replaced by one. Con- 
sideration of cases farther to the north involves the same result with 
track #3. Thus, the presence of many trajectories limits the required 
amount of global searching, an Important saving since the amount of 
computation needed for the optimization (local) increases with the num- 
ber of flight paths. At smaller airports, where the number of trajec- 
tories is small, more starting points will be allowed, but the local 
optimization will require less computation. 

In situations where the search has converged to a local nuninvum 
not significantly different from the starting point (but the existence of 
one is suspected), restarting the optimization after perturbing the 
corner points by some amount (-1,000 m.) may free the search from the 
region of the local minimum, Figure 7.8 shows an example of such a 
situation using a fictitious population distribution. Here, the short- 
range noise effects, as before, create a local minimum for the nominal 
landing path directly over the two lower population clusters. Repeated 
perturbation of the y-coordinates of the corner points eventually allows 
the optimization to find the global minimum. It is well outside of the 
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dispersion limits of the nominal path and the reduction in Nil' is 39%. 
The simplicity of the contrived problem provides for an easy decision 
about which way to perturb the corner points. In a real case, the 
direction may not be as obvious, 





CHAPTER VIII 
RESULTS 


Two computer programs for implementing the system described in 
Chapter VII have been developed, The first, MANIP2, is a revision of 
MANIP (see Ref. [26]), with the following refinements: 1) optimization 
of a flight path's vertical profile is included, 2) takeoff and landing 
trajectories may be optimized, simultaneously, 3) a larger variety of 
commercial aircraft noise data is included, and 4) fuel and time of flight 
have been included in the cost function (as described in Chapter VI). 
The second program, CYCLIC, has the same features as MANIP2, 
except that 1) it employs the equivalent aircraft concept for computing 
noise levels, 2) the optimization scheme is modified to reduce storage 
requirements for large airports (i.e., many trajectories), and 3) fuel 
and time of flight are not included in the cost. A more complete 
description of CYCLIC is given in the Cyclic Optimization section of this 
chapter. Both programs were written in FORTRAN IV for use on CDC 
Cyber 170 series machines. Appendix D contains the source code 
listings . 

The results of several investigations appear in the following sec- 
tions. Two airports were used in these studies: Phoenix Sky Harbor 

International and Lambert-St, Louis International, Their respective 
population distributions are depicted in Figures 8,1a and 8.1b. In each 
grid section in Figure 8.1a (Phoenix) the number shown is the popula- 
tion count for that sector; smaller cell sizes in the St. Louis data 
necessitate (for legibility) the use of population density shading in 
Figure 8.1b. The body of results is divided as follows: 1) takeoff 
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Phoenix (Sky Harbor International Airport) Population Distribution 





InternaUonal Airport) Population Distribution 
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trajectories, 2) a comparison of optimum landing paths using Nil, LWP, 
and HAPN as the annoyance criterion, 3) many takeoff and landing 
paths — cyclic optimization, 4) a study of the effects of including the 
fuel and time of flight costs in the total cost function, and 5) the role 
of the population distribution in determining the optimum paths. 

In all of the follow! ng- sections, the flight paths are represented 
by their ground tracks, with the flight path angles, Yc"Climb, and 
Y^-descent, indicated for each segment. 

Takeoff Trajectory 

Using the Lambert- St. Louis Airport community as a test case, a 
single trajectory was optimized^, with LWP as the annoyance criterion. 
Both the nominal and optimum paths as shown in Figure 8.2. The 
nominal path was approximated by a five segment trajectory, with a 5° 
flight path angle (y); the optimum path was restricted to three seg- 
ments and Y^^ax ^ Nearly all of the variation between the two paths 
is in the ground tracks, This phenomenon, seen throughout the re- 
suits that follow, occurs because 1) variations in the population density 
distribution occur only in the ground plane, and 2) the average dis- 
tances from the flight paths to the population centers change more 
rapidly due to horizontal variations in the flight paths than due to 
vertical variations (owing to the small angles between the ground and 
the trajectories). The exception to this is, of course, when a segment 


^ Aircraft-trajectory assignment data for St. Louis were provided by 
the NASA-Langley Research Center [38] . 









Figure 8.2 Single Takeoff Trajectory - Nominal and Optimum 


110 


near the runway (and closer to the ground) passes close to a popula- 
tion center. Then an Increase In the flight path angle y will produce a 
rapid decrease in the noise level at that center. 

A point approximately 90 km, east of the airport was chosen as the 
outermost endpoint. This allowed the position of corner point #2 to 
vary without being restricted by the geometric constraints. The com- 
parison of annoyance measures is as follows: nominal-LWP = 2571.1, 

Nil ” 0,161; optimum-LWP = 1267,9 (50.6% lower than nominal), Nil = 
0.228 (42% higher than nominal). This increase in Nil, while LWP is 
decreasing, is another recurring phenomenon. It is caused by the 
change in total population used to normalize Nil (i.e. people receiving 

> 55dB), 

Comparison of Annoyance Models 

In Chapter III, it was shown that for a simple population distribu- 
tion (2 concentrated clusters of people), and a straight-line flight path, 
the measures Nil and LWP indicated significantly different optimum 
positions for the path. Here, Nil, LWP, and HAPN will be compared 
tor general trajectories and a real population distribution (Phoenix Sky 
Harbor Airport community), This section also serves to illustrate the 
calculation of optimum landing paths. 

For landings on the runway indicated in Figure 8,3a there are two 
entry points (initial approach fixes). This figure shows the nominal 
ground tracks; the glide slopes are 3°. An estimate of the mixture of 
aircraft types operating at this airport was obtained from data in the 
Official Airline Guide [39]. Because the nominal paths are joined 
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together over a large portion of Uietr lengths, it Is assumed that sepa- 
ration for safety purposes is not a requirement for traffic control here; 

* 

therefore the separation constraint is bypassed. Note that there are 
three segments on each path, with some overlapping. 

With Nil as the annoyance measure, the optimum pair of trajector- 
ies obtained is ^hat shown *in Figure 8.3b, Values of Nil, LWP and 
HAPN for this pair are given in Table 8.1, along with comparisons with 
the corresponding nominal values. 

With either LWP or HAPN as the annoyance measure, the optimum 
pair is radically different, as seen in Figure 8.3c. The annoyance 
measures associated with these flight paths are given In Table 8.1. It 
is seen that the improvement in the annoyance based upon LWP or 
HAPN is greater than that of Nil. Further, Nil increases when LWP 
and HAPN decrease, and vice versa, another indication of the normali- 
zation problem in Nil, 

Fuel and Time of Flight Study 

As an example of how the system may be used to study associated 
problem, the cost of fuel and time of flight (Cp r£,Qp) has been incor'^ 
porated in the total cost function . In this investigation, the relation 
between Cp .pQp and Nil is examined. 

Using the modified cost function, Eqn. (6-13), the optimization was 
repeated for increasing values of Wp ,pQp (the fuel and time of flight 
weighting factor); the airport used was Phoenix Sky Harbor. Figure 
8.4a summarizes the results, showing the relation between the Nil, and 
the fuel and time of flight cost (cp <pQp). Accompanying it are Figures 








Figure 8.4a Nil vs. Fuel and Time of Flight Cost (c. 




8.4b through 8.4f; which show the optimum landing trajectories for each 
value of Wp rj,Qp used. 

As expected, the flight paths become shorter (straighter) as 
Cp ipQp increases. Trajectory #1 shortens rapidly (Wp ,j,Qp - 0.1) to a 
straight path , but trajectory #2 is hindered by two factors: a large 

cluster of people southeast ‘of the runways (see Figure 8.1a), and the 
final angle of approach restriction. For Wp .pQp < 3 shortening this 
path will result in Nil increases that are greater than any decrease in 
the fuel and time of flight cost; hence, there is a flat region in Figure 
8,4a for 0.1 < Wp .pQp < 3. Larger values of Wp .pQp increase the 
relative importance of Cp ,pQp, and the total cost is reduced at the 
expense of Nil. It is impossible, though, for trajectory #2 to become 
completely linear because of the final angle constraint, The effect of 
Cp rpQp is tlien saturated at Wp rj.Qp = 100. The extrapolation for very 
small Wp ,pQp is based upon the value obtained for Nil in the Com- 
parison of Annoyance Measures section (where Wp ,pQp - 0). 

Cyclic Optimization 

At large airports such as Lambert-St. Louis International, the 
large variety of aircraft that must be considered necessitates the use of 
the equivalent aircraft concept in order to keep the computation time 
manageable. Further, the number of trajectories being optimized (15 at 
St. Louis) will require a considerable increase in storage requirements, 
compared with those needed for Phoenix. If the storage required is in 
excess of the allocated resources, the dilemma may be remedied by 
optimizing the flight paths one at a time, cyclically. The power density 
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Figure 8.5a St. Louis Nominal Trajectory Set 




Figure 8.5b St. Louis Optimum Trajectory Set 
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contributions from the (temporarily) constant trajectories are stored and 
then summed with those of the path being optimized in order to deter- 
mine the value in each grid cell. As a result, the optimization 
storage requirements are for only one trajectory, but the annoyance 
measure is for the entire set of paths. 

It should be noted that this approach does not guarantee con- 
vergence to a minimum; however, the results for the St. Louis case arc 
excellent, Figure 8.5a shows the nominal paths at Lambert- St. Louis 
International (tracks 1-12: takeoff, 5®; track 13-15; landing , 3®). 

Each trajectory is approximated with 5 segments. Judging by the 
population distribution (Figure 8.1b), it was decided that paths 1,2,7, 
8,0,13,14, and 15 could give the most improvement in the annoyance 
(LWP). Using the program CYCLIC to cyclically optimize these paths 
(see Appendix D), and allowing a maximum of three segments per tra- 
jectory, the results obtained were: LWP decreased 26%, Nil increased 

1.5%, from the nominal vlaues. Figure 8.5b shows the improved tra- 
jectories. 

Population Distribution and Aircraft Mix 

In order to investigate the relative importance of the population 

distribution vs. the aircraft mix (i.e., mixture of different aircraft 

#■ 

types) on the optimum flight paths, a study was conducted using two 
different population distributions (one fictitious, one real), and two 
different aircraft mixes. 

Using the fictitious population distribution shown in Figure 8.6, a 
single optimum landing trajectory (3®) v;as computed, using the CYCLIC 
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program*® Two different mixtures of aircraft types were investigated; 
the degree of difference is evinced by the difference in their equivalent 
aircraft noise coefficients: 

Mix #1: Cl = 170.70 = 24.23 

Mix #2: Cl gq = 166,00 C 2 gq - 28.60. 

The optimum paths for the two mixes were nearly identical; no tv;o 
corresponding corner points differ by more than 36 m. 

With a real population distribution, though, a strong dependence 
upon the aircraft mix is evident. For the Phoenix airport, two mixes 
were examined: 

Mlx#l: = 169.47 = 24.37 

Mix #2; o^_ = 166.00 Cj gq = 28.60, 

Mix #1 corresponds to the true mix at Phoenix, which was used in 

Comparison of Annoyance Models section and for which the optimum 
results appear in Figure 8.3b (Nil) and 8.3c (LWP). Mix #2 is all 
DC-10 aircraft, with the equivalent amount of passenger seating of Mix 
#1. In the case of Mix #2, the optimum paths are significantly differ- 
ent, with either LWP or Nil as the annoyance criterion. Figure 8.7a 
and 8,7b depict the results. 

The consequence of this example is that the relative importance of 
the population distribution vs, the 'aircraft mix in determining the 

optimum trajectories cannot be stated in general. For a simple distri- 

bution of people, corridors for optimum trajectories that are 


2 Although only one trajectory, with three segments, was being opti- 
mized, CYCLIC was used so that the equivalent aircraft concept could 
be employed . 
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Figure 8.7b Optimum (LWP) Landing Trajectories at Phoenix Sky Harbor Airport (DC-10 mix) 



independent of aircraft mix may exist (such appears to be case in the 
fictitious distribution in this example). A real population distribution, 
though, is usually more dispersed, which results in more people (on the 
average) being closer to the flight paths. Because of the nonlinearities 
in the noise and annoyance models, changes in the aircraft mix result 
in a change in the annoyance of unpredictable magnitude. Therefore, 
new optimum flight paths must be computed when the mix changes.^ 


^ For a particular situation, however, experience may indicate that 
small changes in the mix do not necessitate recomputation of the 
trajectories. 


CHAPTER IX 


CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDIES 

It 

This study has developed a method to assess the ability to reduce, 
through operational controls, the annoyance caused by aircraft noise in 
a community. Test, cases, using a computer implementation of the 
method and actual airport situations, have shown the feasibility of using 
trajectory modifications in the solution to the problem of aircraft noise. 
The amount of reduction that can be gained depends upon the airport, 
aircraft mix, constraints, and model of annoyance being used. 

A comparison study of the annoyance measures Nil, LWP, and 
HAPN has shown that, in addition to resulting in different percentage 
improvements, the optimum trajectories that result from using each of 
these measures can be significantly different. Because of its mathe- 
matical definition. Nil has been shown to be an inferior indicator of the 
amount of change in annoyance for a given community. 

The method has a modular structure which allows for straight- 

s» 

forward modification of any of its components, e.g., the flight path 
model, the aircraft noise model, the optimization algorithm, etc. In 
addition, the method may be expanded to include related problems in 
the optimization. A particular example of this was demonstrated by 
including the "costs" of fuel and time of flight in the cost function. 
The modified system was then used to examine the influence of these 
factors (along with the annoyance) on the resulting optimum flight 
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Of particular Importance is the finding that the population distri- 
bution in a given community does not necessarily determine a unique set 
of optimum trajectories; the mixture of aircraft strongly influences the 
solution. This is unfortunate from a computational aspect, since the 
mixtures on the flight paths generally will vary over time, requiring 
that new optimurti sets be calculated. 

Two concepts that may have applications elsewhere were developed 
in this work, First, the complexity of computing the flight paths was 
reduced by assuming a parametric solution and explicitly including 
restrictions based upon analyses of the linear dynamic equations of 
motion of the aircraft. This technique may be useful in situations 
where the computing resources available prohibit the customary varia- 
tional approach to flight path optimization. Second, the equivalent 
aircraft concept not only reduces the time required for calculation of 
noise levels due to multiple sources, but also provides a means for 
predicting noise levels as a result of changes in the number of aircraft 
on any trajectory, 

The difficulty associated with searching over many local minima was 
overcome by making use of the statistical lateral dispersions in the 
trajectories. These dispersions limit the size of the alternative regions 
(in the horizontal plane) in which to search for significantly different 
optimum paths. 

Finally, this approach to the problem of aircraft noise has shown 
that annoyance can be reduced significantly via operational control. 
The method provides a tool for those who must decide what action to 
take in remedying the problem. It can also be used to optimally site 
an airport in a community to minimize future noise problems . 


Recommendations 

It is felt that further work should employ a population distribution 
model that does not artifically creat optimum pathways. An analytic 
representation of the distribution (e.g. spline functions, or low order 
polynomials) would likely aid in the operation of the method, without 
necessarily increasing comput tional requirements, 

Even with such a refinement, the difficulty of local minima in the ' 
cost function , caused by real, macroscopic variations in the population 
density will persist. Use of the statistical lateral dispersion in flight 
paths has been instrumental in restricting the global search and should 
be automated. 

The assumption of constant runway designations, i.e., takeoff or 
landing, which was made in this study, is not generally true. Future 
work should include the possibility of having to compute an alternate 
set of flight paths, perhaps because of a change in the wind direction 
during the day. The cumulative power density distributon from the 
previous optimum configuration would have to be retained and combined 
with that of the new set in order to include the total effect over the 
period of time in question. It is envisioned that air traffic controllers 
could then use this method to compute alternate optimum flight paths in 
real time. ^ 

Finally, community land use planning could employ this system by 
restricting the maximum noise levels in proposed sensitive areas 
(schools, hospitals, etc.). The cost function would be adjoined by 
penalties for violations of these restrictions. Industrial areas where 
high noise levels from aircraft would have little effect, would be 


weighted less heavily in the annoyance calculation. As the final step, 
the coordinates of the airport centroid would be included as unknown 
parameters to be determined (with some restrictions) in the optimization 
procedure. The result would be not only tl\e optimum set of trajector- 
ies, but also an optimally positioned airport. 


APPENDIX A 


SPECIAL CONSIDERATIONS FOR THE NOISE LEVEL CALCULATIONS 


Energy Addition 

Sound intensity (power/area) is usually expressed in the form 


u 


L = 10 . logj, i 


(A-1) 


where I is the intensity being characterized and Iq is a reference 

intensity (usually taken as the lowest audible intensity at 1,000 Hz, 
-1 fi 

approximately 10" watt/cm^). The units are called decibels (dB) and 
provide a convenient scale for the extremely wide range of intensities 
that may occur. 

For a given volume of space (in which the intensity of sound is 
constant), the acoustic energy received is proportional to the intensity 
and the time of exposure. More precisely, 

E = alAt (A-2) 

where At is an increment of time during which I does not change appre- 
ciably, E is the increment of energy received, and a is a constant that 
depends upon the volume of space being considered. The sound level 
may then be expressed as 


and 


L = 10 • log 


10 E, 


(A-3) 


E = Eq • 10^/^° (A-4) 

where Eq is the energy received from the reference sound. 

This provides a powerful tool for "adding" sound from different 
sources. Suppose that n different sources produce sound with levels L. 
at a particular location. The total acoustic energy per unit volume of 
space is 
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n n t . / 10 

= 1 E. = E„ I 10 ^ (A-5) 

lUi X u 

Eqn. (A-3) is used to calculate the intensity level in dB: 

n L./IO 

I-YOt = ^ (A-6) 

In Chapter III/ the measure uses the energy addition method to 
obtain a weighted average of noise levels occuring over a period of 
time. The use of energy addition is often referred to as adding on an 
"energy basis/' although it is actually power densities that are being 
added. 

Conversion of EPNL to 

The measure Effective Perceived Noise Level (EPNL) is based upon 
Perceived Noise Level (PNL) in the following way (see Ref. [5], pp. 
287-300.): 

EPNL = IPNL - 12 + (onset cocrections) + (impulse corrections) 

(A-7) 

where IPNL is the Integrated Perceived Noise Level 

16 PNL./ 10 

IPNL = 10 1 log^Q Z 10 ^ , (A-8) 

i=l 

The values of PNL are taken every 0.5 sec. for 8 sec., hence the 
upper limit of 16 in the summation. 

Neglecting the onset and impulse corrections (they are not signifi- 
cant for aircraft noise) gives 

16 PNL./IO 

EPNL ~ 10 • logjQ ( Z 10 ^ ) - 12 (A-9) 


. ■ili~ r i i i a i I III I I" 


) 


(A-10 


~ 10 • logjQ 


16 

Z 10 
i-1 


PNL^/10 


16 


The last expression is just the average PNL (on an evergy basis). 
So for a constant PNL, 

EPNt ~ PNL (A-ll 

Finally , for nonimpulsive sound 


giving 


PNL ~ L^ + 13 dB, 
L^ 2 EPNL - 13 dB. 


(A- 12 


(A- 13 




APPENDIX B 
THE LINE SEARCH 


In Chapter V, several of the optimization methods require a one- 
dimensional, or line, search in each iteration. The task is to find the 
value of « that fiinimizes the cost function t, i*e. 


minimize + ttdp 

where xj^ is the current position (at the k-th iteration) in E" and dj^ is 
a unit vector that defines the search direction, relative to 

A successful method for finding a, which has been employed in 

A 

this work, is that of assuming a polynomial form in a for f(Xj^ + 


the exact minimum of which is determined by the calculus of one vari- 
able. The form used here is cubic, i.e,, 


f (xj^ + adj^) = kg + kja + k20(^ + k^a^ , 

for which the minimum occurs at a*, such that 


3£ 

9c( 


a=a* 


= + 2k2«* + akgO^-^ ~ 0 , 


and 


= 2k„ + 6k-a* > 0 . 


(B-1) 


(B-2) 


(B-3) 


The last equation eliminates one of the two roots of Eqn. (B-2), so 

that a unique minimum is given once the k. are specified . This is 

accomplished by calculating f and at two values of a. Denoting 
9f 9f 

f(«o)' ^0' ^0' ^1' respectively, 
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the minimutn, c< 2 ' of the cubic function Is given by the following?^ 

(B-4) 


fft - £. 

"l ° *0 *1 ■ ^ ^ ' 


“2 = '"5 ■ 


“2 ■ “l ■ *“l ■ “0^ ^ 


^1' * “2 - “1 


1 . 


(B-5) 


(B-6) 


■1 *0 

Calculation of the directional derivative is normally performed 

T 

by projecting the gradient = y f(x) onto the direction d|^, i.e, 




(B-7) 


where (*) denotes the inner product in Since derivatives of f must 
be approximated numerically in this work, though, it is simpler to use 
the definition of derivative to obtain 


9a - Aa 


(B-8) 


where dj^ is a unit vector along dj^ and Aa is a small increment. 

Choosing values of a^ and a^ for Use in the cubic fit presents no 

problems if f is cubic: the fit will be exact regardless of what values 
of a are used. With a non-cubic f (as in this work), the value of 02 

A 

given by Eqns. (B-4) - (B-6), may not be the true minimum along dj^; 

however, as stated in Chapter V for the modified D-F-P method, 02 

need only be accurate enough to insure that >0. A logical choice 
for «Q is zero, since then fQ = f(Xj^), a value known from the previous 


^ From Luenberger [23 p. 142] . 
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iteration. If can be chosen such that It and Oq “bracket** a*, the 
true minimum, Is greatly increased, A clever way of achieving this^ (if 

ft 

f is not varying rapidly) is to make depend upon f^ and f^ in the 
following manner t 

=min l-afo/fg', (B-9) 

Where is a constant,' limiting value. The expression -2f0/fQ is 

seen in Figure B,1 to represent the distance from oq ~ 0 to the point of 
Intersection of the tangent to f at wq and the a axis. Taking twice that 
distance for gives a point on the other side of the true minimum a*. 
The resulting cubic fit is shown as the dotted curve, with 02 as its 
minimum. In cases where «2 sufficiently close to a* to give 

m 

g > 0, the cubic fit is repeated, with «q » U2, and computed as 
before. The convergence is quite rapid, generally requiring no more 
than two fits. 

A 

Finally, it should be noted that the use of d^ rather than d^ in 

• A 

Eqn. (B-8) means that M is the magnitude of the increment (Aadj^), 
which should be small to obtain accurate derivatives. Use of dj^, whose 
norm varies between iterations, would give unpredictable results. 


2 Fletcher and Powell [34] first proposed this scheme. 



APPENDIX C 


DETERMINATION OF TANGENT POINTS AND 
. CENTER OF CIRCLE AT A CORNER POINT 


In order to evaluate some of the constraint equations in Chapter 
IV, and for purposes of sampling an aircraft's position as described in 
Chapter VII, it^is necessary to know the cartesian coordinates of the 
tangent points on a trajectory (see Figure C.1). The value of is 
determined by the dynamic and passenger comfort constraints (Eqn. 
4-13), and the corner point coordinates are the unknowns that the 
optimization scheme determines. 


Ground Plane 

First, the value of a is given (by the Law of Cosines) asj 


a = cos"^ { 


a2 + J,2 » q2 

2ab ' 


(C-1) 


where 


The 'distance 1 from the corner point of interest 

R 

1 _ min 

^ "o? ' 
tanj 

so that the tangent points may then be expressed as 

1 




P1P2 


(Xj - x^) + X^ 


(C-2) 


(C-3a) 
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tl = = (n - * H 


(C-3b) 


P 1 P 2 


*t2 “ ^4" ^ 


(C-3c) 


yt2- = 


CVq - Yo) + Yo* 


(c-3d) 


^2^3 


where Xj, are the x and y coordinates of the i-th corner point. 
The linearity of the segments requires that 


■ ^2 


h ■ ^2 P,P 


(C-Aa) 


1‘2 


\2 ■ ^2 


^3 ■ ^2 P„P 


(C-4b) 


2 3 


from which the values of and z ^2 obtained (z. is the z 

coordinate of the i-th corner point). 

Although the coordinates of the tangent points in the ground plane 
and specify the circular arc, it is more convenient to express it 

in terms of R and the coordinates of the circle's center, (x„,y„). 
This point is just the place of intersection of the two radii shown in 
Figure C.l. The coordinates are given by 


b^-bi 
X = ^ — 

c - 1TI2 


(C-5a) 


= 


"*1 “ "’2 


(C-5b) 






APPENDIX D 


MANIP2 AND CYCLIC PROGRAM LISTINGS 
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CQSTi CDSTl # AAC CS # FCKIT# RES ULTKRCSCVP#ON ELTON# FPL y/TPNSPCSi 

FGPA05>GAFA#eFfCP#LINESCH#CHECK#CCPNEP#OYNAPIO#NII#NOISEt 


FAI 

FAI 

FAI 

FAI 

FAT 

FAI 

fat 

FAI 

fat 

FAI 

MAT 

MAI 

.MAX 


eCMMON nTPaJ#NMAP#NSEG#XF( 3#6)# YM(3»6 Jf ZM(3#6I#ARPAY(57t#9)#SP0SITFAT 
U 3#200#3)#X0(3J#XF(3)»YCm#YF (3) fZOO)#ZFO),NP05ITS(3)#RHC( f7EI#FM 
II $EG1# IPRCF#PAXn,PeNCRIT#AWEIT»TVICiA#ALWP# FATCFC# FATOf W MAI 

CCFFON /dynamic/ RADIUS#PWEIT1#PWE IT2#ZWEIT1#GMAMIN#GMAMAX#NSAFP MAI 
COFMON /THRESH/ ALMAXf MAXFL IT# TN6 I TI» TVE IT2#NPLANE FAI 

CCMFCN /CROSS/ X8EGIN» XFlNAli YD1S#CWEIT1#CUEIT2 FAI 

CCFFCN /print/ 5XCENTR(3#3) #SyCENTR(3r3»rSXT!(3#5)#SYTM3#5l# SZTl (MAT 
S3»5 )f SXT2(3#5) #SYT2(3»5)#SZT2(2#5> »SAN6LE(3#3 ) MAI 

CDFHCN /S.AFPLE/ 5AFPT# PE R IOC# SA ML THi NAC ( 3 02 ) # CNOI S ( 16# 2 > MAI 

DIMENSION XNQyOC#U# 0ELTAX(3C)# LABEL(P) 

OIFENSICN XPCRT( 3 I# YPORTO)# JN( 30 ) MAT 

DIFENSIQN NOCfl 3 C( 3 # 2 J# NDC 9 SAM( 3 # 2 )# NOClC( 3 # 2 )# N 707 SAM ( 3 # 2 ) # N 72 FAI 
I 0 ( 3 #n# N 72720 C( 3 # 2 )# N 72 7 S AM( 3 # ?. ) # N 737 I C 0 ( 3 # 2 ) # N 737 S AM ( 3 # 21 # N 7 MAT 
f 472 C 0 ( 3 # 2 I # NL 1011 ( 3 # 2 )# NA 30 C( 3 # 2 )# NBACI 11 C 3 # 2 )# KVC 10 ( 3 # 2 )# NCVFR! 
5990 < 3 # 2 )# NSSn 3 # 2 )# NAMAC( 32 l FRI 

EOUrVALENCe (NAC(n#NDC 830 (in# (NAC( 7 )#N 0 C 9 SAM(Xn# (NACa 3 )#N 0 ClMAI 
loan# (NAC( 19 »#N 7 C 7 SAF(in# (NAC( 25 T#N 72 C(l)i# ( MAC ( 31 ) #N 727200 ( IFAI 
"l)|# (NAC( 37 >#N 727 SAFa)J# ( NAC ( A 3 ) # N 737100 ( I ) ) # ( NAC ( A 9 ># N 737 SAM (IF A I 
m# (NAC(S 5 )#N 7 A 720 C (1 n# IKACI tl) #NL 1011 (in# ( NAC ( t 7 ) # NA 300 ( I ) nM AI 
I (NAC( 73 l#NeACiimi)# (NAC( 79 )#NVC 1 C (1 n # (NAC(B 5 )#NCV 99 Can# (NFAI 

iAC{9inNssTnn fat 

DATA ((CNOIS( r# J)# j»l#2)# I-l#16)/167.7AOP.32#)5A.2e#2A.65#15e.22#MAr 
f2e.PC#lA7. ll#23.fcO# 152.6P#21.95#151.A7#22.6t#l2fc.28# ie.7A,iet .7Af 2F#T 
IQ.*'2#159.31#26.9F# m.B9»2A.7B# lAE . 9A, A9 07F# A7# 37.20# I f 2. 72# 27FAT 
i.Cfc#UA.e 5 # 22 . 95 ,U 9 * 3 A# 27 .X 6 #lA 8 . 22 #lS. 06 / * FAI 

OATA NAFAC/»'DC-«-30 ”#” »»#»»DC-9 W/SAF*»#w enG "#«0CMAI 

1-10-10 '♦#” »'#»'7C7 W/SAH «#»»ENG »‘»”720 **#” MAI 

I ^#*'727-200 "#*» **#”727 W/SAF **#**EN6 "#»*737-FAI 


HC0/20”#”0 "#”737 W/SAF ”#»»FNG 

I ”#”LX011 ”#*♦ *V#**A20C 

n ”#« ”#”vc-ic *'#” 

I ”#”CCNCCROE S”#”ST ♦*/ 

DATA NAC/AeCACy 


**#”7A7-20P 

V pt* 

”#”CV-090 


#*,it 


MAI 


”#”BAC-X1MAI 
*>” F/T 

FAT 
FAI 


N AHE t 1ST / ACFtX / N0CF3C#N0C9SAM #N DC 1C#N7C7SAM# N72C# N7272CC#N727$AFFAI 
I#N7371C0# N737SAM#N7A7200#NL10U#NA3CC#NeAClU #NVCX0#NCV990# NS ST FAI 


1C 

20 

3C 

AO 

?C 

EC 

70 

PC 

OO 

XOO 

no 

12C 

130 

lAC 

ISO 

16C 

170 

XBC 

IOC 

200 

PIC 

220 

23C 

2A0 

260 
270 
2P0 
290 
300 
310 
320 
330 
3AC 
350 
360 
370 
3P0 
390 
AOC 
AlC 
A 2 C 
A3C 
A^C 
A50 
A60 
A 70 
AKO 
A9C 

see 

51C 

520 

530 

*A0 

550 

56C 

57C 

5?C 

500 

6fC 

610 

620 

62 C 

6AC 


OWGINAL 
OP POOP QUA 1 . 1 TV 




ooor>rkr»oraor>r»ooooor>or>orvr> oo o o o o o r» o n o r> o r» o r> o oo o r> r» r> o o 


. INPUT PAP^f'trePSi 

. lENTFY • «FIxeUH FCf? fUEO ENTRY PT5« • 

, • «VR»I4*^ rCR VARUBIE ENTRY PTS. • 

. ICNTINU • 0# INlTIil TP4J*S ARE STRAICMT tIKE5 • 

. ■ If IMTIAI T»AJ*S ARE USER INPUT 

• $AHP1 • TiaE IN SECCNOS BETWEEN NOISE lEVEl 5ANPLE5 

f PERIOD • TINE IN HRS. OVER WHICH THE SIHUtATIQN OCCURS . 

. ACSPO • A/C SPEED IN HETERS/SEC. 

. NTPAJ •• NUPaER CP PIKIHT TRAJECTORIES . 

. NPAP • NUMBER OE POPUIATICN POINTS CN MAP . 

. nsec • NUMBER OF llNf SEGMENTS CN EACH TRAJECTCRY . 

. MAJtn • MAXIMUM NO. OF ITERATIONS AILCWEO . 

• IF CMY PESUIT OF INITIAL CONDITION IS NEEDED# SET MARIT • 0 . 

. OEIH • PERTURBATION TMETERS) IN YiY CIRECTICNS AT CORNER . 

. POINTS FCP CRADIENT CALCULATION • 

. STCPCHO ■ STOP CRITERION FOR SUCCESSIVE CCST CHANGE 

. XCfYCf lOf KFftFflF • STARTING AND FINAL POINTS OF TRAJECTORIES . 

. XFCRTfYPORT • AIRPORT LCCATICN (METERS) 

. XMfVMfIM (IfJ) • CCCRD»S OF J-TH CORNER PT. ON I-TH TRAJ. 

• PVETTl • PENALTY iiEICHT ON OYN, CONSTP. CN TRAJ*S 

. PWEJT2 ■ PENAITY WEIGHT CN CONSTR. ON FINAL ANGIE . 

. RADIUS ■ MIN ALICWEC TURNING RADIUS (METERS) 

. FATCFW • FUEL AND TIME CF FLIGHT WEIGHT 

, XMlNfXMAXfYMlNf YMAX ■ BOUNDARIES OF AREA EXCLUDED . 

. FROM MAX. NOISE (THRESHOLD) CONSTR. . 

. 4LMAX • MAX. allowed A-LEVEL NOISE 

. RATIO • FRACtiCN CF FLIGHTS ALLOWED TO VIOLATE MAX. A-LEVEL . 

. TWEin • PENALTY WEIGHT ON MAX. NOISE AND NUMBER CONSTR. . 

. TWF1T2 ■ NOT USED NOW. BUT SOME VALUE MUST BE INPUT . 

. XBFGINf XFINAL ■ BOUNDARIES CN X-AXIS WHERE SEPARATING . 

. CCNSTR. IS IMPOSED 

. YDIS • MIN. SEPARATING DISTANCE (Y-AXIS) BETWEEN TRAJ’S* 

. cmrifCwEue ■ penalty weights on sep. and crcsscver 
. CCESTR’S. respectively 

. IPPCF • Cf hC PROFILE CPTIHI/ATICN . 

. If profile CPTIMUATION prCLRS 

. NOTE? IPPQF is ALWAYS SET • I IN THIS PRGGPAM 
. zwein • PENALTY WEICHT ON VEPTICAL SLOPE CONSTR. . 

. OMAMAX ■ MAX. ALLCWED VERTICAL ANCLE ON EACH SEG. • 

. GMAMIN • MIN. ALLCWED VERTICAL ANGLE CN EACH SEC. . 

. PENCPTT • HTN. SIGNIFICANT CYN. PENALTY LEVEL 
. AWEIT • SCALING CN MI 

! ""aCMIX ■ NAMELISTI ACTYPE (Ifl) • NO. OF ACTYPE ON TPAJ. I . 

, IN DAYTIME 

• ACTYPE (I#2) • NC. OF ACTYPE CN TRAJ. 1 . 

. AT NIGHT 


CALL OATE (DAT) 

ENCODE UOfOOrOf IDAT) PAT 
CALL REMARK (ICAT) 

CALL JPaRAMS (JN) 

ENCCOE (7f «5020f JOPNAME) JN(U 
call PEMAPK (JCONAME) 

WRITE (6f9U0) 

WRITE (6f903C) 

WRITE (6f9DA0) 4C6NAMEflOAT 
NPLANE ■ 0 
C 

c 


MAT FfC 
MAI FB'O 
MAI (^70 
MAT FBO 
MAT W 
MAI 7CC 
MAI TIC 
MAT 720 
MAI 730 
MAI 7AC 
MAI 750 
MAI TEC 
MAI 770 
MAI 7MC 
MAT 7<0 
MAI PCC 
MAI RIC 
MAI P20 
MM B3C 
MAI PAO 
MAT R50 
MAI 8EC 
MAI 070 
MAI 0FO 
MAI PRO 
MAI POO 
MAI PIC 
MAI 020 
MAI P?0 
MAT 04C 
MAT PfC 
MAI PRO 
MAI P7C 
MAI P0C 
MAI PPO 
MAT lOCC 
MAI ICJC 
MAT 1020 
MAI 101C 
MAI lOAC 
MAT lOfO 
MAT IQEO 
MM 1070 
MAI 1O0C 
MAT lOPC 
MAI UOO 
MAI IllO 
MAT UPC 
MAT 1130 
MAT U40 
MM 1150 
MAT U60 
MAI U7C 
MAT 11 8 C 
MM IIPC 
MAI 12CC 
MAT 1210 
MAI 1220 
MAI 1230 
MAI 1240 
MAT 1250 
MAT T2EC 
MAT 127C 
MAI l?BO 
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c , . >**I 1Z«C 

C . INPUT *U THE iNfOEHmON . P*I IJOO 

C . . P*r IJIO 

PM 1370 

C P4I 1330 

PfPO (7»<;120) (LPEEUnM^Ua) P«t 13«0 

PIAO {T»0050> lEHTPY «AI 13SC 

IP n£M»Y,NE.''V*PU''. AND. lENTpy. HE. "PIKED"! CAU ERPOP (A»"mEr,APAI 1360 

U ENTPV POINT TIPEv CHECK DATA"! HM 13TC 

C NAI 13PC 

C NAI 1300 

C . . . HAT 1*00 

C . ISEOl • 1 • piked ENIPy PIS.. 1ST SEG. INCtUOEC IN All CAICS . NAI 1*10 

C . ISEOl • 2 I VAPIAEU ENTPV PTS.k 1ST SEG. EKClUOEO . NAI 1*70 

C . .NAI 1*30 

C NAI 1**0 

C NAI 1**0 

ISEOl • I NAI 1*60 

IF IIENTPy.£0."VAPIA"| ISEOl • 2 NAI 1*70 

READ (7.*) ICNTINU.SANPT.PEPICO.ACSPO ‘ NAI 1*60 

READ (7»*l NTRAJiNNAP.NSEG.NAKIT NAI 1*00 

READ (7»*) DELH.STOPCHC NAI IJOO 

SAN.LTH • SAMPT*ACSPC NAI 1510 

OC 30 1 ■ IfNTRAJ NAI 1«2C 

READ (7.*) KCm.yolII.ZOd) NAI 1S3C 

READ (7,*) XF(I).yF(II,2F(I) NAI 15*0 

READ (7#*) XPORKDiYPORim NAI 1550 

C NAI 156t 

C NAI 1570 

C . .NAI 1*80 

C . CALCULATE INITIAL CORNER POINTS . NAI 15'0 

C . . NAI 1600 

C NAI 1610 

C NAI 1620 

XM(I.l) . XC( I) NAI 1630 

YM( 1,11 ■ YOI II NAT 16*0 

ZN(I,1) • ZOm NAI 1650 

XM(l,KSEG*ll ■ XF(II NAI 1660 

yN(I,NSEC-*l) • YFIII NAI 1670 

ZNII.NSEGtll • ZFU) NAI 1660 

XN(I,NSFG*2I ■ (XPORTt II.XPmi/2. NAI 1600 

yN(I,NEEG*2l • (rPCRnil.YPIIll /2. NAI 17CC 

ZN( I,NS£G*2I • ZFm/2. NAI 1710 

IF (ICNTINU.EO.il CO TO 20 , NAI 1720 

DC 10 J • 2, NSEC NAT 1730 

^ >M(I,J| • XCni«(FLOAT(J-l)/FLOAT(NS6GI)*(XF(I)-XOml NAI 17*0 

yn(i»ji • vnm*(FLrAT( j-ii/Fio*T(NSEGii*cYftii-yo(iii nat itjc 

ZN(I,J| . 2C(lM(FLO*T<J-lMFLO*T(NSECII*(ZF(II-ZOmi NAI 1760 

10 CCNTINUE NAI 1770 

GO TO 3C . NAT 1760 

20 READ (7,*1 (X6(I,JI,YN( I,J!,2N(I,JI,J.2,NSECI NAI 17*0 

30 CONTINUE NAI 1*00 

OC *C I • l.NNAP NAI 1610 

OC *C J • *.6 NAI 1*20 

ARRAY(I,J| • 0. NAI 1630 

*0 CONTINUE NAI 16*0 

READ (7,*l PUelTl.PVElTE, RADIUS, FATCFW NAI 16JC 

SEAO (7,* I XNIN,XNAX,YNIS,YNAX » NAI 1660 

READ (7,*l ALNAX, RATIO, TI.EIT1,IWE1T2 NAI 1670 

READ (7,6| XEEClN,XFINAL,Y01S,CkElTl,CW£lT2 NAI 1660 

READ (7,*) lPRCF,2WElTl,CNAniN,GMAnAX,P£NCRIT,*UEIT NAI 1*00 

READ (10, •! ((ARMY(I,JI,J.l,3),(*RPAy(l,JI,J*6,9l,RH0III,I.l,N6AP6AI 1<>00 

t| Nat lolO 

READ (7,ACNIXI NAI 1920 







ICBOF • 1 

.00 JC I • IthXm MI 

00 So L • 1. J2 fSI 

NPt*N£ • KPI.*N£+N*C(I»L> f*! 

COMINUfc t"*! 

00 60 I • 1»NH*P f#I 

**R*r(i»io> • 1. f**! 

IF (*Plh.l£.*«BAY(I» l(.ANO.FFP*Yn. II.Le.FPAX.FFO.mh.te.FFPAriPFt 
t l» 2 l.«N 0 .AFFAYn» 2 l •t.£«YnAX.AND.AFF«Y(I» JI.N£«O.I APPAYIIflO) ■ F«t 

» 0. MI 

I CONTINUE N*I 

HAXFUT • inX|NFLAKf***TIO*O.S) NAI 


PRINT INfOFHATICN INPUT 


WRITE (t»9i3C) (i*ean)»i»i»»» - hri 

WRITE (6»9160) NAXlT>NTR«J»NSEC«NnARrCElH*STCPCHC N«I 

WRITE (6»9C6C) S*NPT»P£R100MCSP0 NAT 

IF (IENTRY.EO.''FIXEt«) WRITE I6f9070) NAI 

IF <IENTRY,EO.«VARIA'*» WRITE I6»90E0» NAT 

IF (IPRCF.EO.U WRITE (6»9090I HAT 

IF (IPRCF.EC.O) WRITE <6»91C0» . NAT 

WRITE <it«l90l ALNAX>NAXFUT>XHIN>XNAX>YNIN^YHAX NAT 

WRITE (6»9160» XBEGIN»XF1NAI,TDIS NAI 

WRITE I6»9170l RA0ILS*PW£IT1>PW£IT?»ZWEIT1»TW£ITI»TW£IT2»CWEIT1*CWNAI 

fcIT2>GNANINK<HAnAX«PENCRIT.AWElT>FATaFW NAI 

00 90 I • X»NTRAJ NAI 

WRITE (6»9180J I> X0( 1 1 > YOI I ) > I0< 1 1 »XF< 1 1 «YF ( I) > IF (tl NAT 

NSEGl • KSEG«1 NAI 

DO 70 J . I.NSECl NAT 

WRITE <6.92001 J.Xnn.JI.YNn.JI.INfl. Jl NAT 

3 CONTINUE NAT 

WRITE (6.92101 NAT 

00 flO I • 1.31.2 NAT 

IF (NAC(I.U.NE.0.0R.NAC<I»L*1 » .NE.OI WRITE (6.922CI NANACUi.NAT 

1 NAMAC<L*n.NAC<I.U.NAC(I.t«l» NAI 

: CONTINUE NAI 

WRITE (6.91901 XPCPTdI.YPORTdl NAI 

3 CCKTINL'E NAI 

N ■ 2*(NSEG-1I»NTRAJ NAI 

NVAR • 2 NAT 

IF dPPOF.EO.l) N . 3*(NSEG-1»*NTPAJ NAI 

■ IF dPROF.EC.U NVAR • 3 NAI 

OC ICO I • l.NTRAJ NAI 

NSEGl • NSE6-1 Nil 

00 IOC J • 1. NSEGl NAI 

00 IOC K • l.NVAR NAI 

L • d-U*NVAR»NS66U(J-n*NVAR«K NAI 


IF IK.EO.n XNCWd.n 
IF (K.E0.2) XFOW(L.l) 
IF (K.E0.3) XNCWa.l) 
C!3*'TINUE 
OC 110 I • l.N 
OELTAXTH • pEtN 
CONTINUE 


XN( I. jn I 
YNd.J.ll 
2Nd. j*n 


STIRT CPTINIZATICN 




aVV* 


149 


C Z97C 

c . **#I ?!»c 

C*U ONtWTOW «l‘*»IT»STCPCM»N»)lH0W*0tlT»)l»3C»?CNCI»IH ?5«IC 

SfCP ?fcco 

c rn 26 K 

»CIC FCSH3T UlOl NT 26ZC 

9C2C fO»rM <*7I NT «« 

9030 fC(f|**T «//.•• NC *rBlfNT NOTSE IMClOO£D"»//» »**r 26*0 

9C*C fCFK*T ,»»/»lOX>".»'>UX»".**»y»lOX»''.'*»2X»67»2HM 269C 

»W, f IX# »•»*•»/» 10X» •."» 1 IX* "•"»/* If*»I 2660 

*0X>n. .#•.•••••••••*#/> T*XT 2670 

9050 F0R6*TH*9» , I'** 2660 

906C fC«H*T «6X#*'60F ION CXlCUtATlCN# SANLE TI»t • ".N.C#" SeC."»/#:XN*I 66 C 0 
»,«*/C LANDING PATTEFN FEPEATS rVfSt «#P5.I»" HPS*"#/#SX#*'A/C SPtfPNAI 270C 

« • »#f6.2#" Hn"ittn NAt rno 

9070 FCRNAT t/#*X#»riXtO ENTPY POINTS**#/! NAT 27^0 

906C FCRNAT «/#*X#**VARlA£LE ENTRY PCINTS**#/! N*I 2730 

9C9C FCPNAT j/#*X#**.. PRCflLES DPtlNIZEO .***#/) NAI 27A0 

9100 FOPNAT T/#AX#**#. CNIY CRONNO TRACKS AR£ 0PIINI2E0 ..**#/) NAI 2750 

9110 FCRNAT aHl#IOH****»*****#///l NAT 276C 

9120 FOBNAT IBAIOI NAI 2770 

9J3C FCRNAT nX#6A10l NAI 2760 

9UC FORNAT I///#5X#*'NAKINUN. ITERATICN SETt **# 12# /#SX#**NUNBER OF TRX"#''NAT 27»C 

fJECTCRIESi **#I2#/#5X#**NUHBER OF SECNENTS ON EACH TRAJECTORY! **#I1#NAI 2*00 
1/#5X»**NLNBER OF PCPUATION POINTS ON tHE NAPl **# IA#y»SX#**PERTliR6 TNAI 2P1C 
ARAJECIORIES IN X AND Y DIRECTIONS **#F10.5#** NETERS FQR CAICULATINGNAI 2P20 
1 GRA01ENTS**#/#9X#*'£XIT CRITERION FOR CNEHTONi N#IPE16.9#/» NAI 2630 

915C FCRr*T T5X#*<N0 6LCCK NAY RECEIVE **#F5.1#** 06 (OR HIGHERI A-LEVEL**#NAI 2«AC 
1** NOISE HCRE THAN **.I9.*< TINES A DAY**. //* 3X. **. .THIS CONSTRAINT ".**NII 2*50 
lEXClLOES BLOCKS WITH CENTROIDS INSIDE THE RANGES! **#F8.1#** TO **#1 ENAI 266C 
f.l#** (X) AND **#Fe.l#*' TO B.FB.I#** (Yl***//) NAI 2670 

OlfcC FOPNAT (6X#A0HWITHIN THE AREA OF X»CC0RCINATES EETWEENfFP #1# 5N ANCNAI 26*0 
» .F6.1,*1H# THE separating DISTANCE BETWEEN TR A JECT#5H0RI6S# / #5X# INAI 269C 
»*FSNCLLC BE AT L£AIT#F6. 1#7H NETERS#//! NAT 2900 

917C FCRNAT (/#eX#**OTHEP CONSTRAINT PARAHETERS**#/#*X»**NININUN RAOILS**#**NAI 2*IC 
% CF CURVATURE (IN X-Y PLANE PROJECTION! • **#1PE9.2> /#AX»**0YNANIC PNAI 2920 
tENALTY WEIGHTS! PWEITl • **#E9,2#" PWEIT2 • •*»£9.2#** IWEITl • **#NAI 2930 
*E9.2./#AX#*’HAXINUH NOISE LEVEL (THHESHOLOl **»*'PENALTY WEIGHTS! TWNAI 29*0 
lEITl • *'. t9.2»**TW£IT2 • *'#E9.2» /» *X.**S£P*RATION PENALTY WEIGHTS! NAT 2«50 
tCWEITI • **#£9.2#" CWEIT2 • "#E9.2» /#*X#"GHANIN « "#0PF5.1#" GNANANAT 2960 
»X • '*>F5.l#**PENCR IT • "#IPE9,2#" AWEIT • "#E9.2# 9X# "EATOFW • **#E9,NAI 297C 
»2#//#5X#"INF0RNATICr< CF EACH INITIAL TRAJECTCRY ! "#/ / ! NAI 2980 

9160 FCPNAT (10X#15HFL1GHTPATH NC! • II# /#12X#13HINITIAL X#Y# 2# IX # 23HCOCNAI 29«0 
tPOINATES IN NETERS! # 3(FB . I #3X! # /# 12X#11HF IN At X#Y# 2# 1X#23HCQOROINNAI 30CC 
TATES IN NETERS! # 3(FE . I# 3X !# /# 1 2X # 1 AHIN IT I At C0RNER#18H POINT POSINAI 3010 
..TTinNS! #/#IAX#10HCCRNER NC.#6X# IHX#9X#1HY#9X# 1H2# /! NAI 3020 

9190 FOPNAT (I0X#"90NWAY LOCATION! " #f * . 1 # I OX# F9 .1 # / / ! NAI 3030 

9200 FOPNAT (lBX#Il#tX.Fe.l»3X#F£.l#3X#FE.l! NAT 3C*f 

9210 FCRNAT (15X#"A/C NIX"# T*5# "NO. OF A /C"»//»T*0# "OAYT INE"# 3X# "NIGHT TNAI 3050 
TINE"#//! NAI 3060 

9220 FCRNAT (T20»2A1 0# T*2# 13# T52 » 1 3/ ! NAI 3070 

END NAI 3CPC 


SUBROUTINE COST ( IGRAC. TCTAL# ANIl # PNAlTY#CLOSE. THRESH! CCS 10 

CONNON NTFAJ#NHAP#FSEC*XN(3#6l#YHI3*6!#2H(3»« !#ARRAY(!76#9!#SP0SITfCS 2C 

T(3#2CO#3l»X0(3l#XF(3l»YO(3l#YF( 31 » 2C( 31 # IF( 3 1 #NPOS ITS( 3) #RHC( 5761 »rr$ 30 
TISFGl.lPRCF#NAXIT#PfNCFIT#AW£IT»rVICLA#ALWP#FATOFC»FATOFW CCS 90 

CCNNON /CCRNER/ ANGIE(3#3!#P0SIT(3.200.3!#XCENTP(3#3I»YC£NTR(3#3!#C(’S 50 

tZCeNTR(3#3).»Tl(3.5l»YTl(3#5!»2Tl(3»£!#XT2(3»5!»YT2(3#5l#2T2(3»5l#CCS 60 

TNPCSir(3l,0IS(1 1 CCS 70 

CCNNCN /PRINT/ SXCtMP(3»3!»SYCE*!TR(3»31#SXTI(3# 5!# SYTI ( 3 » 5 l#S 2 n(rcS 80 
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ONE 

ECO 


OIHENSION 500(30.30) 

ONE 

Etc 


0|j^5NSXON 0(30,1} 

ONE 

220 


n ■ 0 

ONE 

230 


CALL COSTl (O.N.rNCW.XNCU.ANll.PNALTY. CLOSE. THFESH.NOEC) 

ONE 

240 


IF (PNAITY.GE.PENCRIT) GO TO 10 

ONE 

250 


CALL HONIT ( 1 T. FNOW. N.FNCW.ANII. PNALTT. CLOSE. THRESH. NOEC ) 

ONF 

260 


!f (HAXXT.COrO) STOP 

ONE 

270 


GC TO * 

ONE 

200 

xo 

bPITE (6.0020) PNAITT.ANII 

CNF 

290 


CALL CCPNER (XNCW.)..NDEC) 

ONE 

300 


IF (H4XIT.E0.0) STOP 

ONE 

no 



ONE 

320 



R»i 0^E 

330 



. ONE 

340 


STEP, 11 SET S • IDENTITY HATRIX ANC CALCULATE GRADIENT 0 

. ONE 

35C 



t ONE 

3EC 



■I. ONE 

370 



ONE 

300 

2C 

DC 30 I • l.N 

ONE 

390 


DO 30 J • l.N 

ONE 

400 


S(l.J) • 0. 

ONF 

4X0 


IF (I.EO.J) S(t.J) « 1. 

ONE 

420 

3C 

CCNTINUE 

one 

430 


CALL F6RAC (N. FNOb. XNCb.GNOW. DELTA X.NCFC ) 

ONE 

440 



ONE 

450 



1.. ONE 

460 



. ONE 

470 


sr^p 2» SET 0 • -SG 

. one 

400 



. ONE 

490 

0^»^5 51C 

40 

CALI HPIT TN#N,1#SaGNQW,0#N0EC#N0EC#1I 

ONF 

920 


OQ 5C I ■ l#N 

ONE 

930 


ona> ■ -o(ia) 

ONE 

5 40 

50 

CO^TINU£ 

ONF 

990 

OKe 5ftC 


t 

STEP 31 LINE SEARCH ALCNG 0 TO FINE AF* THAT SATISFIES P0>0 . OHA 

. ONF 

ONE 

V ONF 

K • 0 CNF 

K • K»1 ONE 

CALL LINFSCH «K »N» FNOV» XN0H»D»AFA2» TTEHP* FSM ALL. ANtI» PNAL TT»C LOSE »ONE 
ITFPESH.NCEO ONF 

CALL FGP40 (N.FSHALl.ATENP.CN'EAT.DELTAF.NOECI ON'F 

00 7C I • l.N ONF 

ou.n • GNEXTn.ii-GNcvn.n one 

P<hll • AFA2»0(I,1» ONE 

CCNTiNUE ONF 

CALL TPNSPOS (N06C.1.P.PTI ONF 

CALL MPLY (l.N.l.PT.O.PO.l.MDEC.U ONF 

IF (PO(l.l) .GT.O.I GO TO 00 * ONF 

IF (K.CT.3I PEIUPN , ONE 

FNOW • F5MALL ONE 

DC EC 1 • I,K ONF 

»N0N<I.1I . TTFNPd.ll ONE 

CONTINUE ONE 

GC TC AO ONE 


'll 1 ‘tfl'itti* ^fr 




oor>r>Oir>or><^ 


c 


c 


c 

c 

c 

c 

c 

c 

c 


90 DC ICC I • l#H 

. KNcwn#u • )(T£Hp<i»n 
100 CCNTIKUC 
IT • n>i 

PRCNT • ResnPSP^iU-PNOKWPNOW) 

IP (PRCNT. Ce.STCPCHO GO TO 120 
wRm (6r9C30> PPCIT#STOPCHG 

IP (PNALTy.c£.PCNCprn go to uo 

CAU NCMT nT^XPOV#N,PSHAU#4NnfPPAlTyKlOSEfTHP6SH#NDPC) 
R£TUPN 

UC ►RITE IGf^OlOl PNAUY^ 

CAU CQRNeP (XN0k#N#N0ECI 
fltETlRN 

120 IP (IT. IT. MAT IT) GO TO lAC 
MPITE'(£»90A0) 

IP (PNAITY.CE.PENCPIT) GO TC 13C 

CAU MONIT (IT#TNOW#N^PSMAUpAMI#PfvALTY#CtC5e#THPeSHfN0EC) 


ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 

CNF 

ONE 

ONE 

ONE 

ONE 

ONE 

ONE 


RETURN " ONE 

I3C WRITE (6#9010) PNAITY ONE 

CALL CORNER (TNCWf N.NOEC ) ONE 

return one 

UO IP (PNALTY.GE.PeNCRlT) GO TO 150 ONE 

CALL MCMT ( I T^ TNON#N» P5NALL# AMI . PNAL TY# CLOSE #THPESHrNDEC> ONE 

GO TO UO ONE 

i5C WRITE (e»9020) PNALTYrANII ONE 

CALL CCRNER (RNCW #N#NDEC I ONE 

16C CCNTINUI CNE 

ONE 

. . ONE 

. .ONE 

. STEP XF »IT»» IS INTEGER MULTIPLE CE N GC TC STEP U . ONE 

. IE NOT# update S • ONf 

• IS ”IT»' INTEGER MULTIPLE OF N? . ONE 

. . ONE 

ONE 

ONE 

IP ((FL0AT(IT)/PLCAT(N)).NE.PL0AT(IT/N)) go to 17C ONE 

PNCU • FSMaLL ONF 

GC TO ?0 ONE 

ONE 

• ONE 

. . ONE 

. UPDATE MATRIX $1 CC TO STEP 2 • ONE 

. . ONE 

. ONE 

ONE 

i7C CALL MPLY ( N#N. US # C, SO#NOEC# NOEC# 1) ONE 

CALL TPNSPOS (NDEC# I #C, QT ) ONE 

CAU MPLY (N#UN#SC#OT#5C0#N0eC#l#N0Ea ONE 

CAU MPLY (N.N^N#SCCrS#SCCS#NOEC#NOEC#NOEC) ONE 

Call MPLY (1#N#N#CT#S#CS.1#N0EC#N0EC) onp 

CALL MPLY U#N^ l#CS#C.CSC»l#NCECf 1) CNF 

CALL MPLY (N#l#N#P#PT#PP#NaeC#l#NOEC) CNF 

00 UO I « 1#N ONP 

DC IPC J • l#N ONE 

S( I# J) • (S(I#J)-SQOS( J, J)/CSC( Ull )*(POa#l )/OSQ(U U )^PP( I# JI/ONC 

% PO(l#l) ONE 

180 CCNTINUE one 

FNCW • PSMALL OFF 

DC UO X • UN CEP 

GNCWU.l) ■ ONEXTU#!) ONE 

19C CCNTIME ONE 




PIO 

E 2 C 

R30 

RAO 

E50 

RED 

R70 

0PO 

090 

900 

91C 

«20 

9|0 

9A0 

9fC 

9ec 

970 

960 

990 

1000 

loic 

1020 

1030 

lOAC 

1050 

lOfcO 

107C 

1060 

1090 

IIQQ 

mo 

1 I 2 C 

U3C 

IIAO 

U5C 

lUO 

U7C 

U6C 

1190 

1200 

1210 

1220 

1230 

12A0 

1250 

12 EC 

1270 

126 C 

l?9C 

ncc 

T3IC 
1320 
MIC 
13A0 
1350 
II EC 
1170 
1 3P0 
1 39 C 
lAOC 
TA 1C 
lA?r 
1A3C 
1 AAC 


CO TO 

C . 0H 14^0 

SOIQ <Mi40HU FESyiTl DYNAMIC CONSTI»AINTS VtOlATTQN^ /# XOY ^CHPeKCKf 1470 

lAlTY •flPElei*)) Cf.t |APC 

«?C20 fOPMAT HOKTTI OYNAMIC CONSTPAlNTS VIOIATICN# A# lC)(»«?MPINAC^f 14«>0 

HTY •rlPU6.4>/rlOFf*'ANn • . 0.# THEN ANII HA5*»# »*NOOKF 1?00 

IT g|SN COMfUTdD IN THIS ITIPATICKi QUI TO lAPCE PNAITY»») ONE 1«10 

<502C FGPPAT (2)tf-PCt4MVr CHANCE IN SUCCESSIVE COSTS »*rlPElC.|f USS*‘QN^ 1S?0 

tKAN STOP CfltEPlCK ♦♦flPnOiSI OTI 1550 

9040 FOPHAT (2Y»29HfAYirUH iTEFAnON SET FEACHEOl ONE V?4C 

END CNF 15?0 


C 

C * 

C * 

C . 

C , 

c * 

c 

OIMFNSIGN 4aOECrHCEC)i MHOECfNOECIi CaPEC#NDFC) 
DC 1C r • IH 
DO Id J • V>K 
cnf>l) - 0 . 

00 to K • S^H 

ca#j) • till jHAtb KMa(K,j) 

K CmiNUfc 
FEtUFN 
END 


SUftFCUTlKE MPIY < S#C#lDEC>MOECi NOtC) 
CAiCUUTE MATRIX NUlTIFUCAnCN C • AC 


MFL 10 
HFl 2 C 
. HPl SO 

. NFl 4C 

• MFl 50 

• MPl 60 

• PFl 70 

PC 

MPl 90 
MPl toe 

MPl no 
MPl UO 
MPL UC 
MPl i4C 
MPL tfC 
MPl tfC 
MPl 170 


C 

C 

C 

C 

c 

c 

r 


SUaRQUTINE TRNSPOS IM,N#A,|$J 



. TPANSPOSF OF matrix A 1$ R£TURN£0 |N MATPlX & 




D|Me^S^CN AIM|N)| 
^00 It I • UM 
DC U J - WN 
QlJiH • AUfJI 
Id CCFTIM'E 
PtniFN 
bNO 


TPN to 
TFN ?C 
TRN 30 
TFN 40 
TFN »C 
TFN 60 
TFN 7C 
TFN PO 
TFN «C 
TFN ICC 
TPM 110 
TFN 120 
TFN 130 
TFN lAC 
TFN 15C 


SiPFCUT IN£ FG«AD <K#F*XrCfOaTAX,NCEC) TCP IC 

t FFF 70 

r ..M -•••«•••**«•««• FCF 30 

C . • FCP AO 

C . CAICUATI ORAOUM CT CDSt F VltH RESPECT TO UNKNOWN X • FCp fO 

C f • FTP 60 

C ..... FCP 70 

c POP PC 

DIMENSiCN XlNOFC^lh OiNOfcCMUf OF ITAX tNDFCf U FOR «C 


00 10 I • l$H 
, CU>1) • 0* 

|HOO(If3).£OtO) CC rc 10 
MI#U • KU* U^OItTAXm 

cm C05TI a#NifF#)i#INIUPN^LTVKlOSE#THF£5H^NO£C> 
OUfU • <PP-Fl/0£ITAKU) 

K(I»AI • Kn#U-DELTAMn 
CCNTINUI 
«erM#N 


fC» 120 
FOP 130 
rci^ 140 
FGP 150 

rcp uo 

pop 170 
FGP ICC 
FOP l^C 


SU6FCUT1NE CAFA (Nf F#OAFA^G#Xi»Of NOFC) 


CALCULAYF OPAOIfNT CF F With RFSPfCT TP AFA 


oiHfcNsioN xeNO£c#u# x>itit>tCp\), XTcrpnoui 

00 10 I • 1#K' 

nEhp<r#i) « xn#ntpAFA^on#u 

10 CONTINUE 

CAU COSTl (WN,FF#XT£MPMNnfP4AlTY^CLCSE#THPESH> 

G • (FF-FMOAFA 

PETIFN 

END 


SUBPCUTINF FRROP JK, MESSAGE) 
DIMENSION MeS$AGE(12)# ITEXTUZ) 
DATA 1TEXT/12*« 

00 1C I • If K 
10 ITEXtm • MESSAGEtn 

WRITE (6f9010) meXTUif IulfK) 
STfP 
C 

901C FORMAT (//flXfMEPPCP J **#1ZA10#//) 
END 


IPO 10 
IPO 20 
IPO 30 
IPO AC 
IPC 50 
IPO 60 
IPC 70 
IPO CO 

I eo 90 

IPO IOC 



SfPPCUTlNE UNESCH <K # N# FNO V# KNCWf Df AFA2, XTEMP# FSMA LL# ANI I # PN AUY# UN 
1CtCSEfTHPE3HfN06C) UN 

OIMENSION 0(NOEC»n# XT£ MPIfO EC f I ) f XNqw(NOEC>l) UN 


CORK FIT 6 Y INITIAUY lETTING AFA • 0 AND OAFA • 5 


DNCRM • 0. LIN lie 

OC 10 I • IfN UN 12C 

DNCRM • ONORH^CUf 1)**2 UN 13 O 

CONTINUE I IN l AC 

ONCRM ■ SCRTIONCRn UN IfO 

IF (ONORM.EO.O, ) CAU CHECK (1 ^ N# f NOW, XNOW f Of Al AOf FAFAQ# GAFAOf AFA It IN UC 


UN 130 
UN lAC 




llFAF*^•GAfAl#Ul|tU^U^iNCCC) 


UN 
UN 
. UH 

. N0RKAU2I 0 • UN 

, . UN 

UN 

UN 

Do to I • i#N UN 

on/n • onuMCNCPN un 

2C CONTINUE tIN 

OAFA »n5. HN 

AFAO • Of 

FAFAO • FNOW UN 

KK ■ 1 UN 

CALL OAFA (N#FAFAO^OAFA»OAFAO,yNOVfO»NDfC» UN 

IF (GAFAO.FOtO) CALL CHFCK U#N#FN0W»FN0U#0# AF>0# FAFAO# OAfAOf Af AUUN 
«FAFAl#GArAUUl#UUlH»a#NpCC) UN 

IF f-FAFAO/CAFAO.lF.Of ) CAU ERROR (6|« -FAFAC/CAFAC .l£. Off TRAJI IN 
t» 5FT PROBAFLY NFAR NIMHIN FT.**) UN 

AF*l • ANINmOOOif-FAFAC/GAFAOl UN 

30 00 AO I * IfK UN 

YTENPUfl) • >NCWMf JHAFAMOUfU UN 

AC CONTINUE UN 

CAU CnSTi (IfNfFAfAlfXTEFPf ANnpPMlTYfClOSEiTHRFSHfNOEC) UN 

IF IKK.EO.AI CAU 6RRCR (7f»'PR0BtE« IN LINESCH. CANU FIND LOWER UN 
tCCST AFTER 3 ClfilC FITS»’) UN 

CAU GAFA (NfFAFAlf DAFAfGAFAlfKTlNPfOfNOEC) UN 

IF tAeSUAFAO-AFAUFAFAD.LTf pOOn CAU CHECK <2# N»FNCW> XNOWfOf Af ALIN 
lOfFAFAOfCAFAO#AFAlfFAFAUCAFAlfUl#UeifU2#NOeC) LIN 

U1 • GAfAO^0AFAl-3.*(FAFAO-FAFAl) /f AFAC-AFAl) UN 

UU ■ Ul»*2-GAFA0<GAFAl UN 

IF (UUl.LE.O.) 02 • Of UN 

Ir aOl.GT.Cf) 02 • SORTdUn UN 

AO • AKINUAFACf AFAU I IN 

A I • AMAKUAFA0fAFAl» UN 

GO • GAFAO LIN 

FO ■ FAFAO LIN 

r,I • GAFA). LIN 

FI ■ FAFAl UN 

IF lAFAO.UfAFAn GO TO 50 LIN 

GO • OAFAI LIN 

FC • FAFAl I IN 

Gl ■ GAfAO tIN 

FI • FAFAO LIN 

50 IF <Gl-G0#2.AUP.EO.0.l CALL CHECK (3# Nf FNQWrXNOWr OrAO# FO# GO# Alf Fl f L IK 
""5Gl »Ul#UUl#U2iNDEC I UN 

AFA2 • Al-Ml-AOl A(Gi^U2-Un/(Gl-GO*2.^U2 ) UN 

DO 6C I • l#N UN 

XTENPdfl) • XNOWd# lMAFA2*0n#n UN 

EO CONTINIE UN 

CAU COSU (OfN#FSNAU#XTEMP#ANn#PNALTy#CLCSe#THRESH#NOECI UN 

WRITE U#9ClCi AO# Alf AFA2#ONOPh#FC#Fl#FSNALL#GO#Gl UN 

WRITE LIN 

IF dSNAU.GE.FKOW) GO TO 7C LIN 

return LIN 

70 AFAl • AFAl/2. UN 

KK P KK*1 • UN 

GC TC 30 I IN 

LIN 

GOJO FCRFAT d5X#“A0 • IPE13.6#5X#»<A1 • »♦# E V3 f fr # 5X #» AF A2 • # 6 1 3 . 6# 5 XLIN 

IffCNORfi • «#F13.6# /#15X#*‘F0 - w#E13 .fe#5X# «Fl - *»# El 3 . 6 # 5X # »*F$ 7 A LL LIN 
1- tV#£l3.F#/#15X»‘*G0 • »*#En.6#5X#»»Gl ■ LIN 

902C FCRHaT (5X#»<0 • ♦V#lGax#iPE‘;.?f?X n lin 


JOC 

?00 

21C 

220 

230 

RAC 

2«C 

260 

27C 

2RC 

2«C 

300 

no 
120 
330 
3A0 
35C 
36C 
370 
IRC 
390 
AOO 
A 10 
A20 
A3C 
4AC 
AfO 
4FC 
47G 
410 
49c 
?0C 
510 
520 
53C 
54C 
?5P 
56C 
570 
BRC 

frOO 

61C 

F20 

630 

640 

65C 

660 

ETC 

66C 

600 

7Q0 

710 

72C 

730 

7AC 

750 

760 

770 

7A0 

7PC 

ACC 



ENO tU' tIO 

SUgECUTINE CHECK I ICHECK* N» FNOW>KNOIi«0#«C>FO«CC»« l»Fl»ei>Ul>UtI>U2rHE 1C 


IfNDECI 

CHF 

20 

OIKEKStCH XNCUtSCEC*!)* 0(N0EC>1) 

CHE 

30 

hPITE (6»10l FNOk#N«n>XKCH(I«ll»0(l4U>I>l4NI 

CHf 

AC 

IF nCHECK.EC.l) STCP 

CHE 

EC 

WRITE I6.2Q) AO,FC.CO»*l»Fl»r,l 

CHP 

EC 

IF (ICHECK.E0.2) STCP, 

CHE 

7C 

WRITE I E. 30» Ul.inJl.U2 

CHE 

PO 

STCP 

CHE 

OO 

C 

CHF 

ICC 

10 FCPFAT IZOX.RWHTHIS IS SUIRCOTTFE CHECK WHICH GIVES RU 

THE>35H INCHE 

no 

IPQfNATKH IK SUBROUTINE UNESCH..//.30K.7HFNCW • .1PE16 

37KfiCHE 

12C 

IHXNCW.17Y.1H0. / /»• (25X.I2.3X.lPElE.<).AX»lPeiE.F. / 1 1 

CHE 

130 

20 FCRMAT I//.20X.NA0 • «. 1PE1E.9. 5X.-F0 ■ ".E1F.9»SX."C0 

• ^#Eie.9f/CHE 

JAO 

»2CX,«41 • «.EU.9»SX»''F1 • «.EU. 9. 5X»"G1 • ".E16.9I 

CHE 

15C 

3C (//#20X,«HU1 * # IE 5)r# lUl • f lPei6*9#5X^ 5HU2 ■ #lPEieCFF 

uo 

I.P) 

CHE 

17C 

END 

CHE 

leo 



SU3RCUTINI CORNER IX.N.HDECl 

CCP 

1C 


COKKCN NTRJU.NH*P.NSEC.XK<3.6|>YK(3.6l»ZHI3.e>.«RRRY<976.9t.SP0SITCrR 

2C 


f (3.200.3)»YCI3I.XFOI.YOI3I.YF<3).2C(3I»ZFI3».NROSITSI3».RHOI576I.CCR 

3C 


SISESI.IPROF.HAXIT.PENCRIT.AWEIT.TVIOLR. ALUR.FRTOFC.FRTOFU 

crp 

40 


OIHENSIQN XINDEC.l) 

CCP 

50 


NVIR • 2 

CCP 

FO 


IF (IRROF.EO.l) NVAR • 3 

CCP 

70 


OC 1C 1 • l.NTRAJ 

COR 

eo 


NSfOl • NSEG-1 

CCP 

OO 


00 10 J • l.NSEGl 

CCP 

IOC 


00 lO K ■ l.NVAR 

CPP 

nc 


L • (I-U»NVRRAK5EGl*(J-ll*NVAR.K 

CCP 

12C 


IF (K.EO.l) XKII.J+ll . KU..U 

CCP 

130 


IF IK.E0.2) YHII.J.ll • XIL.n 

CCP 

14C 


IF (K.E0.3I mi.jAii • xa.n 

CCP 

X5C 

IC 

CONTINUE 

CCP 

IfC 


0 0 3C I • l.NTRAJ 

CCP 

17C 


WRITE (6.9CICI I.XHII.ll.YNII.U.ZKlI.n 

CCP 

IPO 


DC 2C J ■ 2. NSEC 

COP 

190 


WRITE I6.9020J XNII.Jl.YHn.JI.ZrM. J» 

CCP 

200 

2C 

COMINLE 

CCP 

21C 


WRITE (6.O020) XH(I.NSEG*lI.YH(I.KSEG*l».ZKM.NSEC*l) 

CCP 

22C 

3C 

CONTINUE 

CCP 

230 


RETURN 

CCP 

240 

C 


CCP 

250 

9C1C 

FCRHAI I12X,"FLIGHI PATH KOI "» 11 »/. I*X. "CORNER PTS."./.UX." 

(»%F7CCP 

?6C 


l.C. F7.0>",".F7.C."I") 

CCP 

2 70 

9020 

FORMAT IIRX,"(",F7.0,«,«,F7.0»".".F7.0."I«T 

CCP 

2P0 


END 

COP 

2 CO 


sueroutine oynario IPNAITYI 

DYH 

1C 


CCPflCH MRlJ»NK4P»NSEG.)rH{3.6>»VM(3»tl»2M3»tli.*«PAY(576»<J)»SFOSnCyH ?C 

»(a.2CC»3»/»C(3l «>f O)»Y0(3l »YF(3»»2Ct3).lF<3»»NPOSITS(3»»PHC(gr6l»tlYN ?C 

»t SYfil»lPFOF»H*Y IT,p£SCPII>AViEIT>rVIClA,AtWP*FATCFC» FATQFli PYN *0 


coHhCN /oYNific/ RACtus.PkemjPwtnz.rfceiTweNiwiH.cHifAXiNsif'p oy>» 
. CCffCN /CORNER/ *NCIE « 3r3>» POSI T( 3» E0C> II , RCENTRI 1*3 >ftCE MR ( |» 3> ,PYM 
IEC£NTRII,3l»3ri(3r5IMTU3»tl»3UI3,5l»)(T2(3»5l#n2(3»5l»2T2(J.5l»PYN 


iNPOSITOlfCIEIII . PVN 

REAl LENGTH P»N 

l£NOTh(R»8.C»0»E»PI • SORT! l*-ei*«2»(C-0l**ER«E-F|*P2) OVN 

PI • ATRNa.I*Rf RTN 

C OYN 

C PVN 

C . . PVN 

C . FIPST CALCVLATE PARANETERS AT CORNER POINTS . DTM 

C . " , . OYN 

C NYN 

C OYN 

00 6C I • l.NIRAJ OYN 

00 10 J . 2.NSE6 OYN 

YT2U.JI • C. OYN 

YTUt.JI • YTZn.J) PYN 

KT2II/J) • YT1(I.JI OYN 

PTUl.JI • KTid.JI CYM 

2TUI.JI • 0. OYN 

ZTZd.J) • 0. PYN 

ANGLEII.J-U • KTUI.JI OYN 

YCENTR(I»J-1I • ANGLE<I»J»ll OYN 

YCENTRn»J-M • YCENTRI I. J-ll PYN 

10 CONTINUE PYN 

OISII) • 0. OYN 

TTHI.ll • XNd.ll OYN 

XTEd.n • <Nd.ll OYN 

YTld.ll • YNd.ll OYN 

YT2(I.H • YNd.l) OYN 

ZTUI.I) • JNd.l) OYN 

2T2(l»n ■ ZNd.ll PYN 

ISEG2 • ISEOUl PYN 

00 AO J • ISEC2.NSEG CYN 

A . LENGTH (<Nd.J-ll.XNd.J).YNd.J-l I, YNd.JI.ZNd.a-D.ZNd.rYN 

i Jl) OYN 

e • lENGTHlXNd. J).XNd.J»ll.YN(I.J).YHn>J»l).2Nd.JI.ZNd.J*CYN 

t 1)1 OYN 

C • LEN6TH(<N(I»J-l)»>Nd»JYll.Y«d.J-U»YNM» J*I).2Nd. J-n.ZrYN 

t Nd.Jfll) OYN 

C PYN 

C PYK 

C . . PYN 

C . 00 THE INC SEGNENIS PORN A STRAIGHT LINE? . OYN 

C u . OYN 

C CYN 

C 0 YN 

IF (AeSIA.p-Cl.CT.l.l GO TC 20 PYN 

SEG • LENGTH(<T2d» J-i)*XHn.JI.YT2d» J-1 1. YHd. 31.212 d.J-n.PYN 
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ZTUI#J) • ZH(l#J)^ALCNG 2 *fZHn>J-l)- 2 «CI#J)MAXY OYN 
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f l#ZTUI#J)> TYN 

DISN 2 !* AltNG 2 PCOS( 0 ETA/ 2 i) DYN 

OISC 2 • PE 0 ltSZSlN( 8 ETA/ 2 .> CYN 
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AC CONTINUE OYN 
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OYN 
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CXICULATE VERTICAL SCOPE CF EACH SEGMENT AND . DYN 
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i TlMAMAXl<C.#TANOCGrAMlN»-VSLDPf )) 4*2 DYN 

50 CONTINUE OYN 

AXY • LENGTHC XHCI #NSEG ) # XN Cl# NSEGf 1 ) . YM C I ^NSEG) # YM Cl #NS EG ♦ 1 ) # C. # PVN 
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PNAtTY • PNALTYtPVEIT 14 CANAXl CO. # C 4 C 0 . ♦ ALONGl-A X Y ) ) ) • * 2 1 PUE I T 2 * COYN 
« AhAXICO.# (?. 3 fC 5 -ANGn )442 OYN 

DISCI ) ■ DISC nnCNGTHCXT 2 CI#NSEG),XHCr#NSEG 4 l)#YT 2 n#NSEG)#YMCI 0 YN 
I #NSEG^U#ZT 2 C IiNSEG ) # ZMI I #NSEG# 1 ) ) OYN 

FATOFC ■ F 4 T 0 FC#f ATOFW*OISt n DYN 

6 C CONTINUE OYN 

RETURN DYN 
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HA 
HA 
HA 
HA 

HA 

gYTfPNAL ACV»Y HA 

INTIOEft OAO HA 

tOClCAl HAGlCBfAgWTPAj^IgNAm HA 

CCHHCN NTPAJ,HHAP>N$eG#>H(U^13)>tHM5#U)WHU5#13WSPCSn(U#250HA 
l#3N)f0n5J#)'F(I5)fYCfl5hYPaMf30U5)fZA < 15) i NPCSlTSUf )» AHa(973 )HA 
$(il5EOli IPACF/MA)Ur# PeNCPlTiAVgXt# TVICiA|AtWP^P60PU#NTACPT^THPl.0N(HA 
lS73)iaPO<15)rlTGf ICPTrITfl€PClCBiNEVtfA4#APPAY(571ilO) HA 

CCfHCN /GYNAHIC/ PADIUSi PWB ITl# Pwf lU/Zvei Ti#6HAH|NU5)#GHAHAA(l5 )HA 
5fKSAHP HA 

CCHHCN /THRBSH/ AirAY#MAXFUT#TWimiTWen2rNPlANK HA 

CCHHOH /CPOSS/ XCEClKiKf IKAUVOlSfCWIITliCUElTZ HA 

CC^HCN /PPINT7 S>C£MPn5»12) >$ YCENtP(15#12)^$XTm5#l?)f SYTU15#1HA 
i2)#SZnM5#l2)»SXt2(l5#U)rSYT2U5fl2)#5ZT2n5#12)#$ANGUfl5n2) HA 
CCHHCK /SAHPU/ 5AHPT#PERlOD^$AHLTH^NACU5r32)#CNClS(U#2) HA 

CCHHdN /INAett/ ZENABtt HA 

ClHtHSltN KHOwriC#!)# OElTAXOCIf lAPEtrA) HA 

CI-ENSICN XPGATaj)# YPOPm5)f JN(30) HA 

OJHEI^SICN KOCl30a5>2)» NCC9$ AH < 15 » 2) ^ N0C1C(15#2)^ N7Q7SAH(15^ 2)#HA 
I H72GU5#2)# N7272CCfl5r2)# K7E75AH ( ) 5# 2 ) ^ N737100( I5r 2 N7375AH(HA 
il5i2)# N7<r7200n5fE)# NUOn(15#2)/ NA3Q0<15^2)# N0ACUXU5#2)f NVHA 
ACICU5#2)» NCV990f 15*2)# HS$m5#2)# NAHACrSZ) HA 

ftOuiVAUXCE rNAGm^HDCB3G(l))j ( N AC (3 U# N0C9SAH ( 1) )# OACUDiNPCHA 
llO(l))» (hAC(9U#N7C7SAHa) )r IHACnZI) »K720iU I# |NAC (151 ) #H7272QHA 
iO(n )> rKACaAUiN727$AHm )# (HAC(2U)#N73710Cm)» (KAC(2Al)iN72HA 
A7SAH(D), (NAC(27U#f7A72CO(U)# (NiC( 301 ) aHUOI M 1) ) # (NACOaUiHHA 
iA300 tll)> |NAC( 3fel)»NeACUl(n ) # (HAC(391) #NVC10MI)i ( NAC( A 21 lf HCHA 
iV990MI)# iHACIASn^HSStd)) HA 

dATA NAHAC/»’CC-e-3G h,«oc«c vySAH*«»« ENG *‘#*«OCHA 

i-ic-)o *•#« »v*ac7 wsAH **pt*^m ha 

1 ••|»'727-|00 ^/$AH «i»*ENG «,«'737-h| 


ncc/?o»»#«o 


**#«737 V/SAH »«>«fNG 


H/H7A7-200 


HA 


it 


290 

300 

310 

320 

330 

?AO 

350 

3fc0 

37C 

3HC 

390 

4CC 

910 

A2C 

A3C 

AAC 

A5C 

AfC 

A7C 

ARC 

AOC 

5CC 

510 

52C 

530 

5AC 

550 

500 

57C 

fPf 

«0C 

0OC 

010 

620 

630 

6AC 


c 

c • 
c • 
c • 
c . 
c . 
c . 
c . 

C i 

t • 

c . 

€ i 

c « 
c • 
c • 
c ♦ 
c « 
c • 

C I 

c • 
c • 
c • 
c • 

c . 
c . 

C f 

c . 
c • 
c • 
c • 
c . 
c . 
c # 
c • 
c • 
c * 
c « 

C .» 

c > 
c » 
c • 

c • 
c • 

C f 

c . 
c • 
c . 
c « 
c • 
c # 
c * 
c . 
c • 
c . 
c . 
c • 
c . 


• «#nion •»#h>ko 

II 1*11 

I ^tnoHCom 5^#«sr **/ Pk\ 

okn Nic/ifO*o/ ra>t 

nmiur nmifi Nccnc#Noc^$i>‘#Nocxo#NTOT$>^^#Nnc#N7e?20o#NTi7$i?»***i 

f#N7371OO»AT37S4^i#N7%7eOO>NUOU#MICO##<0ilCm#NVeia#NCV9^O#NS5T 1*17 

Hmiuj (QPCiPf cPO#Mxna#ziNmi nn 

i^n 

Pi 


*n 

• 7*11 
. PM 

• PM 
. PM 

• 7*A! 
» PM 

• 7»Ar 

. Mr 

• PM 
. PM 

• PM 
» PM 

• 7* AT 
. Tt A f 
. PM 

• T^AT 

• T^AI 

• fAI 
f 7* A T 

• HAT 
. 7*AI 

• 7<Ar 

• 7< AT 

• TiAT 

• T^AT 
. T-A l 

• 7>'AI 

• PM 

• PM 

• PM 

• PM 

• 7* AT 

• 7»AT 

• 7* A I 
. 7*AI 

• PM 

• PM 

• FAT 
. 7* A T 
i 7*At 

• 7* #T 
f ?*AT 
. rij 

• fAT 
. PM 

• PM 

• 7* At 
. 7*At 

• 7*AI 

• hat 
. HAT 
. HAT 


INAiTT AAAAHITIASi ^ 

OAOIA • NAT^tUSTi CAOill 

HA Kite 

aNASil 


iiT^Tpy 

IC7<T17<U 

SAHPT • 
PfRICO 
MtPO • 
WTPA4 
NHAP 


TIAJ* TIDt klTH CAOlHAi I 
TOTAl NO# {JA filCIAt mAAtlOT^S 
.lAtr.# 7-C0dA03 CA CCAHti ATS# 
mOkfO TO VAAV 
» tAALSe.# I-COCAO$ MHO CONSTANT 
f HHKgrw PCR FlKeO INTAY ATS. 

• hvaIIah FCA VAAime intay ats, 

• Ca INiriAt TRAJS AAI STIAIOHT LlNtS 

• W INITIAH TRA^S MB USCA INPUT 

TiHr IN SCCCKOS eiTVISN NOISE UVEt SAHAtlS 
- TIME IN HAS* CViR WHICH THE SIHUIATICN QCCUAS 
A/C SAEEC IH Pittn/SH* 

NUHEER Cr AU6HT TRAJICTOAIES 
NUH0EA OF PCAUfcATTCN POINTS CN HAP 


NSIO •* NUH0EP OP IINE SieTlfNtS CN EACH TRAAECTdAy 
HAKIT • HAKIHUH no. CP ITIRATICNS AUOHfO 
IP ONLY RiSULT DR INIUAl CCNCjnCN 13 NEICIO* HAKlt • 0 
CCLN • PEPTUfiP/ATICN IHeTERS) IN K*Y OIRrCTIONS M CCRNEA 
PCINTS PCP CMDIENT CAtCULAtICN 

STDPCHC • STOP CPnCPION POP SUCCESSIVE COST CHANCE 
KCVYCiZC#KP^YP#lf • STARTING AND PINAL POINTS OP TAAJECTOAIES 
kPCPTiYPCPT ^ AIPPOPT IdCATICN fHETERS) 

XPMPMP itpM • CCCPDS OP J-TH COPNPP PT. CN 1-TH TAAJ. 

PWEITI f PENALTY WflOHT ON OYN. CONSTP. ON TPAJ5 

PWEIT2 • penalty WEItHT ON CONSTP* CN PINAL ANCLE 
PAOIUS ■ HIM ALLCWEC TUPNINC RADIUS 77'ETEPS) 
KHIN^KHAKfYHlN^yrAy • EOUNOAftlES CP AfEA EKCLLCEO 

fPOP 7*4 K* NOISE (THPE5H0L07 CONSTR. 
ALHAK • HAK. AUCWEC AHEVEL NOISE 

^ATir • PPACnOK OP PLIGHTS ALLCWEO TO VIOLATE HAK. A-LEVEL 
TVPiu • penalty weioht on hax. noise and NUHREP CONSTR. 

UaiTE • NOT USED NCV* CUT SOHE VALLE HUST EE INPUT 

KEICINiKPINAL • feOUPCARIIS CN K-AKIS HHERE SCPAPATINC 
CONSTP. IS IHPOSEO 

YOlS • PIN. 5BPAKATIHC DISTANCE IT-AKIS) 6ETW6EN TPAJS* 
C>EITl#CwEITe » PENALTY WEICHTS ON SEP. ANO CPOSSCVPP 
CCNSTPS. RESPECnVELY 
IPPOP • 0# NO PPCPUE CPTITtlZATION 

Ip PPOPIU CPTIHIZaTION OCCIPS 
IPPOP IS ALkAVS SET • I IN THIS PRCCPAH 

• PENALTY kElC-HT ON VERTICAL SLOPE CONSTR. 

. HAX. AltCKCO VEPTICAL ANCLE CN EACH SEC. 

• HIN* ALLCkEC VEPTICAL ANGLE CN EACH SEG. 

* HIN, SIGNtPICANT OYN. PENALTY LEVEL 
SCALINO CN MI 


NOTE* 

Zk?m 

OHAHAK 

GHAHIN 

PENCRIT 

AkEIt • 

ACM IK ■ 


nahelistj actype riMJ 


I 


NO. OF actype on TRAJ. 

IN OAYTIHE 

actype llfZr ■ NO. CF ACTYPE ON TRAJ. I 
AT NIGHT 

GHAHIN a I ■ HIN. PLIGHT PATH ANGLE (DEG I ON 
TRAJ. 1 

CHAHAK m ■ HAK. PLIGHT PATH ANCLE (CEG» ON 
TPAJ. I 


EfO 
Rec 
67 C 
EPC 
6RC 
700 
71C 
7?C 
730 
7AC 
7E0 
760 
77C 
70C 
7RC 
POO 
eio 
P20 
«7C 
6A0 
850 
PEC 
R7C 
PPC 
P«0 
ROO 
<710 
770 
RIO 
OAO 
R50 
R60 
0 70 
RPC 
ROC 
KCO 
1010 
1070 
1070 
10*C 

10 SO 
1C6C 
1C7C 
lOPO 
1000 
noo 
nio 

M2C 
II 30 
IIAC 

11 EC 
iiec 
117C 
VIPC 
11 CO 
IPCO 
1310 
l??0 
1?30 
Ut-c 
13*0 
1260 
1370 
13FC 
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c 

c 

c 


• » 

# PAI 

’ " ^ . .. T r 

C^Li UCCV* f«CVM>7T»»C) 

riP 1 

PAI 

CAU OAII fOATI 

PAT 

EHcoee (vo»9qio»io*n o*r 

PAT 

CAU AlPAAP MOAT) 

PAI 

C*u JMRAfS (JN) 

PAI 

encoci (^(qoEO/JcifoKeJ jn(i» 

PAI 

CALL KIPA6K (JCfeNAPf) 

PAI 

KKITE 

PAI 

PPITfc (6# *030) 

PAI 

mAIU (6»*0*0) 4CINAPI/1DA7 

PAI 

NPii^t • 0 

PAI 


PAI 




. PAI 

9 IHAUT all the INFCAPAflON 

# PA! 


• PAI 




PAI 

tfiC (7»«150» (U6Hm#l«i»8> 

PAT 

*fAC (TiOROEI>l 

PAI 

»E*D t7»0CS0» lENTPY 

PAI 

IF (!IHIFr,NE,<'V#RI*t'.*FO.IENT«ir.N6,HFJFeo«) 

CAU £6606 <*»P|iUOAPAl 


H CNTRY PCtM TYPb CHICK CnAt*) >1! 

Pkl 


• ••••«* i • i •*••***»*•***** <^ *••••* 9 •••»*••*• P « * • « • 9 f p> « f • 9 P i • i • • « i i « » •» 


* ISfeCl • I * fnu IKTRY PT5.I 1ST 5IG. INCLVOID IN ALl CiUS • 

9 I$fCI • 2 i V#MMU IN1PY Pn., 1ST $10. IKClUOlO 


itm • i 

IP n5NTPYp|Ci'*VAPlA»*) iseoi • 2 

«|AC (7iH ICNTlNU^SAJ^Pt^PeUqOMCSPD 

PIAO (?#•) NTAA4#NrprPT#KTPCNST#NHAP^N5fC>NAYIT 

P IAO (7#* I DllHrSTCPCHC 

SA.HTH • $*nPT*ACSPC 

OD 30 I • I9M9A4 

«r4D (?/•) YCiCPOun>YocoPorrnizorcPfHin 
P2A0 (7r*I YHCP0(m#YMCA0(Un2M0P0nn 
riAo r7r*i xpoPncpom)#YPCPticRcm> 



• CALCULATE IMTIAL CCPNIP POINTS 


it-(nPDMra) • xcicRDmi 
YP(o«o(n#ii • YO(OPtm)) 
iniopcm#!) ■ ic*LPD(tn 

ifHi0PC<IbN$|O4l) 9 XF(QPOm) 

Yt(CR 0 (I)#N$ECil) • Yf(C»0<in 

IH«CP0m^ASFG4l) • ZFfOPCm) 9 

YK(0PO(I)iNSE0^?) • LYPOPUCPOni KF <0PC 1 U )) /2. 
Yrt(OROtn#NSHC4 2) • {YPCPTtOPOrilMYPIdPtHl)) >/?p 
zr(CPDm#N$ic4pj p zFLOPom)/2. 

IF UCNTlNU.iO.l) OC TO 20 
DO 10 J • 2 /NSEC 

XHUfJ) * tci IHiFLCAK J-U /FLOAUNSFCnA(XF( n-YOCin 
YMII# J) • YC( l)4(PLCATU-U/FtOATINSPc n»<YF rn-YD( n ) 


pn 
PAT 
HA I 
Pat 
PAT 
PAT 
PAI 
PA I 
PAI 
PAI 
PH 
PAI 
PAI 
PAT 
PAI 
PAI 
PAI 
PAI 

PAT 

PAI 

PAI 

PAT 

PAI 

PAI 

PAT 

PAI 

PAI 

PAT 

PAT 

PAI 

PAI 

PAI 

PAI 

PAl 

PAI 

PAI 

PAT 


u«o 

HOC 
I A 10 


mo 

1110 

IIAC 


1110 


ilfO 

1170 

mo 

HOC 

1*00 

1*10 

1*10 

1*10 

1**0 


1*fC 

1*10 

1*70 

1*10 

l*flO 

1«0C 

151C 

IS2C 

1530 

15*0 

1550 

1 A ^ A 

1570 

15FC 

15*0 

ucc 

UlO 

lft?c 

1630 

16*0 

1650 

1660 

1670 

mo 


16*C 
1700 
1710 
17?0 
173 C 
17*0 
1750 
1760 
1770 
17P0 
1700 
IPCC 
IPIO 
1620 
1630 
16*0 
1650 
1660 
1670 
15 60 
1**C 
1000 
1«10 
1020 




m 
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c 

c 

c 

c 

c 

c 

c 


. • 70m*»FlO*T»J-i)/FLO#T(NSI:C))»(ZFn)-ZOnn fM 

coNTiNue 

r,c 1C 30 t‘n 

»f/t0 «7»*» »XI^«0^0(l)»J>*VH(0lt0n>»J»»Z»*«0fDn»»J>*J*2»NSEGI fll 

CONtINUe »*4I 

00 *0 I • l»NN*P P*I 

CO 40 J ■ 4>6 

*pf*y»i»j) • 0. >*M 

CCKTlKOe fH 

nto I7**> PWfiin»PW£IT2»P40IUS PII 

ff.no (7»*i )iniN*xr*7»¥»'iN»YH*x nn 

PE4C (7,*t 4LI‘4X»»iTrD>TkeiTJxTfcEITZ P*I 

Ptxo «7/*> XBE6IK»XfIN*l»yOIS»CWEIIl»CW6tT2 f*t 

READ t7»*> lP«OF>ZbEin»PENCRtT¥iWEir f*PI 

REAC I7>*) ((CN0IS(CP0<n»Jt>ail>2l>I*l>NTRPJ) HRI 

READ U0>*> ((APRAYdf 3 )»(APPAV( !*<))« >I>6>P)»RP0(I)« 1*1 tN’KAPPAI 


READ (7«ACHIXI PAI 

READ <7»*) <6MAHlNICP0»lM»I»l»f<TR4JI - PAI 

READ (T»*( <6MxHAX(C»0nn»l»l»NTPA4l »AI 

IPPCF • I HAT 

DC 90 t • IfNTRAJ MAI 

CO 90 L ■ 1>32 MAI 

NPLANE • NPLANE»NAC(UU MAI 

) CONTINUE MAI 

DC EC I • UNMAP MAI 

ARRAYiUlO) • 1. MAI 

IF I*MIN,lE.AFPAY(I.l(.AND.ARPAY(I»l).lE.XMAX.ANOtYMIN.LE.AFRAY(MAI 
1 If2I.AND.APPAYU*2).LE.¥MAX.ANO.APRAY(U3).NE.O.) ARRAYIUlO) • MAI 

1 0. MAI 

3 CONTINUE MAI 

MAXFLIT • IFIX«NPLAN£*PATID*0.9I MAI 

MAI 

MAI 

. MAI 

PRINT INFORMATION INPUT . MAI 

. MAI 

MAI 

MAI 

WRITE (t»9I6C) iLAeEL(I).I*UP| MAI 

WRITE »6»9I40) <CPOm»Ul»NTPAJ»,MAXITG MAI 

WRITE U.917C) MAXIT.NTRAJ.N$EG>NMAP.OELH«STQPCHG MAI 

WRITE <6.9180» NTRAJ.NTROPT^NTPCNST MAI 

WRITE I6.9060J SAMPT.PEPIOO.ACSPO MAI 

IF IIENTPY.EO.t'FIXEC") WRITE (6,9070» MAI 

IF (IENTRY.E0."VAR1A») WRITE (6.9C8CI MAI 

IF (IPROF.E'O.l) WRITE (E.9090I MAI 

IF MPROF.EO.O) WRITE (t»9l00l MAI 

IF (2ENABLEI WRITE (E.9IIC) MAI 

IF ( .NOT.ZENABLE) WRITE (E.912C) MAI 

WRITE (6.9I90I AIMAX.MAXFL I T» XM IN . XMAX, YMIN. YMAX MAI 

WRITE (6.9200) XBEGIN.XFINAL.YOIS MAI 

WRIIE (E.92ICI RAOIUS'.PWEITUPWEITZ.ZWEITUTWEITUTWEITZ.CWEITlfCWMAI 
IEIT2.PENCR IT.AWEII MAI 

WRITE (6.922C» MAI 

WRITE (6.9230) (ORD( I ) >GMAH IN (ORO ( I) ) i GMAMAX ( CRD ( I ) ). I • I. NTRA J ) MAT 
DC 90 I • I.NTRAJ MAI 

WRITE (6.9240) ORO ( I ). XO(CRD( I) ) > YO(ORO ( I ) ). 20 ( CRD ( I ) ) . XF (CRD (I )M AT 
> ).YF(CRD(I)).ZF(CRD(I)) MAI 

t'SEGl • NSEG.l MAI 

NSEC-2 • NSEG*2 MAI 

CO 70 J • 1.NSEG2 MAI 

WRITE (6,9260 J. XM)ORO ( I ). J) . YM( CRD ( I ) , J) , 2M(CP 0 (I ), J ) MAI 

1 CONTINUE MAI 



r> r> o ri o o o 


ao 

<?c 


wRITit ( 6 f 9270 r 
DO 60 L • lp^X^t 

IF <K4Cn^U.NE.0.06*NAC(I#L>U.NE.O> WRITE (^#9260 

I NAHACU^Df KAC(I#U#NACU»I^U 

CONTINUE 

WRITE ( 6 f 9250 ) FFCRT (ORDI 1 Hf YPORT(ORO( I M 
CONTINUE 
N • Z*(NSEC-n 
NVAF • 2 
IF nPRCF.EO.U 
IF U>RQF.E 0 . 1 J 
I « 1 

NSEGl • NSEG -1 
DO 100 J » i#N$EGl 
OC too K • IfNVAR 
t . n -1 >^NVAR*N$EGl*IJ-n^NVAR^K 


N • 39 <NSEC-n 
NVAR*« 3 


100 

UO 


IF TK.EO.l) 

IF (K.E0.2) 

IF (t<.E0t3) 
CCNTIFUE 
00 no I ■ IfN 
DElTAXdl ■ DEIH 
CONTINUE 


xNawTun 

XNDVU#n 

xNcMun 


KpicPoaiFJn) 
YP(CR 0 m#J 4 U 
ZNIORDI I)#JfU 


START OPTIMTZATICN 


OC 120 K ■ IfNHAP 
TPPIOMKI • It 
120 CCNTINUe 
NOEC • 30 
ZEPCaCB • .TRUE. 

CALL PNGPOBN {NAX ITGfSTCPCHGfNf XNOWfOEITAX^NCEC ) 
STOP 


9GIC PCPPAT <A10) 

902C FCPPAT lAT) 

Q03C FQPMAT (//#« NC ANBUNT NCISE INClUOEb»S / / ) 

FCFHAT (lhlp9X>M /#! GXf» .«f UK # ” f / f lOX f »». «f 2X t A 7 f 2NA I 

iK#«.'S/FlO>F«.»*FnKF«. »$ /f ICK# A IOfIK# ”.»*#/ UKf »* t ?»f / f IN AI 


NAT 

MAT 

NANACIUfRAT 

NAT 

NAT 

Nat 

NAI 

NAT 

NAT 

NAT 

NAI 

NAT 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

• NAI 

• NAT 

. NAI 

NAT 

NAT 

MAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 

NAI 


9 LJSC 

90 EC 


fOKi ••.n*F/) 

FCPNAI tA5> 


NAI 

NAI 


FOFHAT (SXf^FQP IDN CAICUIATICNf SAtPlE TINE • ** fF4 .0 f ** S EC . »♦# / f 5 XNAI 


1f'U/C LAhDING PATTEFN PEPEATS EVERY ^fFSIfN 


9070 
9C8C 
9Q9C 
91 CC 
91 1 C 
9I2C 
91 3C 
914C 

91 5C 
016 C 
917C 


HRS.”f/f 5 Kf»»A/C SPCEPNAI 
NAI 
NAI 
NAI 
NAI 


1 • »*fFF.?f« N/S»V///) 

FCFNAI (/FAKF’^FiyeO 6KTFY POINTS”# /I 
FCPNAt (/F9KF”VAPlA6tE ENTRY PCINTS”#/) 

FOFHAT (AfAKf”.. PPCFHES QPTINIZEO ..”f/ 1 
FCRNAT t/FAKF”F. ONLY GFGUNO TRACKS ARE OPTINIZEO *.”#/) NAI 

FORNAT (/fAKf”.. Z COORDINATES OF CCFNER PTS. ALLOWED TO VARY”f//INAT 
FCPNAT (/fAXf'S. Z CCOFOIFATES OF CORNER PTS. HELD CONSTANT”#//) NAT 
FCFNAT (IHXf 1 Oh ••♦♦♦♦ •*♦♦#///) NAI 

FCPPAT (/ f5Xf”IMT]AL CROER QF TRAJSl ”# IM 12 f 2X ) f /# 5X# ”NQ. OF GLCNAI 


SHAL ITERATIONS 
FORNAT IPAIO) 
FCPNAT IIXfPAVO) 


n?) 


NAI 

NAI 

NAI 


FORNAT I///fSXf”NAXINUN ITFRATION SETi ”f T 2 f /f f X f ”NUNB ER OF TRA”f”NAI 
tJECTORlESl ”f I 2 f/f 5 Xf”NUNBER OF SFCNENTS CN EACH TRAJECTORYi ”f 11 fNAI 
S/^1Xf”NUMBER of RGPLLATIQN points of THE NApi ”# lA# /f 5 Xf”PEPTUPB TNAI 
iPAJtCTCRIFS IN X ANC Y OIPtCTlOFS ”fF10,5f” NETERS FOR CALCULATINGNAI 
% GFADIENT5”f/f£Xf”6XIT CRITERION FOR CNEWTONi ”fIPE 16 . 9 f/I NAT 


2970 

256C 

299C 

?eoo 
2610 
762C 
263C 
26 A 0 
2650 
2660 
267C 
26P0 
2690 
2700 
2710 
2720 
2730 
27A0 
275 C 
2760 
2770 
27P0 
2790 
2«00 
2P10 
2P20 
2«^0 
2«A0 
2850 
2860 
2870 
?88C 
2890 
2900 
2910 
2020 
2930 
29A0 
295C 
2960 
2070 
2cpO 
2090 
300C 
301C 
3020 
3030 
30AC 
3050 
3060 
3070 
36P0 
3090 
3100 
311C 
312C 
313C 
31A0 
3150 
3160 
31 7C 
21 6 C 
3190 
320C 


r><-> r> o o o o o 


9160 FCRMAT (7^5X#«TCT*L ^UheER CF TFAJ6CTCRIFS • 2 ^ /, «)T# »’NiniEF CF 

ITPAJECTCRUS TO BE rPTIHIZeO • 1 2# 5X^ »»NUPeFR CF TRAJECTCRUS FKRI 

IHO CONSTANT • m 

919C FCRHAT (5X»»’N0 BtOCK MAY PBCFIVE «»Ftai« OB (DR HIOHEP) A-LEVH«>MRT 
»»* NOISE HCR6 THAN Tl^ES A 04 Y»»# A/f 3X# **• . TH IS CONSTRAINT w,i*rAI 

IBXCLUOES BLOCKS VITH CENTROIDS INSIDE THE RANCESI TO »»#FFNAI 


lai” (X) AND w#Fe.l»« Tt (YI»«W7) 


,HAI 


9200 FORNAT (5X740HW1THIN THE AREA CF X-CCCPOINATES BETWEENae. ANOHAI 
• #F6.1#VIH# THE separating DISTANCE BETWEEN TRAJECT#5HORIES#/#5X# IHAI 
lanSHDULO 6£ AT LEAST# F6.l#7H HETEPS#//) HAI 

o?lC FCPHAT (/#eXf mother. CONSTRAINT PARArETERS"#/# AX#**NININUH RAOIlS»«f **NAI 
I CF CURVATURE (IN X-Y PLANE PROJECTION) ■ «# 1RE9.2# /# AX^WDYNAriC PHAI 
SENALTY WEIGHTS* PWflTl • «»E9.2#»* PWEIT2 • «iE9.2#»» IWEITl • »*#NAT 
IE9.2#/#AX#»*NAXIHUH NOISE LEVEL (THRESHOLD) *«#»’PENALTY WEIGHTSI TVNAI 
$EIT1 • •’#E9a#«TWEITZ • »»#B9.2# /# AX#»SEPARATICN PENALTY WEIGHTSt HAT 

scwem i ♦*#E9,2a‘ cwenz ■ «#e9,2#/,ax#»» pencpit • «#ipe«. 2 #»» awfnai 

»IT • *^#E9.2#//) NAI 

9220 FOPRAT < /# 6X#”VEPTICAL SLOPE CDNSTHAINTS«# /# 8X#^TPA J N0.»*#2X# ••GHANNAI 
fIN CNANAX«#/n FAI 

9230 FCPNAT (1 IX # 1 2# 5X# F5 . 1 # 3X f F5. X » HAT 

92AC FCPHAT (J OX aSHFL IGHTP ATH NOt # 12# / # 12X# 1 3HIN ITI AL X# Y# Z# IX# 23HCQCH At 
SPDINATES IN RETERS* #3 ( f 9a #3X ) #/ # 12X # UHF INAL X# Y# Z# lX#23HC0CP0IKrAT 
SATES IN NETEPS* # 3 ( F9 • 1# 3X) # /# 1 2X # I AMIN IT 1 A L CCPNER#ieH POINT POSIHAT 
STIONS* #^# 1AX#10HCCPNEP N0.#6X# IHX# 9X#IHY#9X, IHZ# /) HAT 

025C FCPNAT (10X#»*RUNLAY LOCAT lOKl »*#F 8 . I# lOX# FB« 1 #// ) MAI 

9260 FORMAT ( 1 EX# 12# fcX# F9, 1# 3X# F9.1# 3X # F9 . 1 ) MAI 

9?70 FORMAT (15X#»»A/C MI X«' # TA5#«NC. CF A/C«# //# TA0# WDAYTIME"#3X#«NIGHTTMAI 
SIME»^ // ) MAT 

92PC FCPMAT (T2C#2 AJO# T A2# 1 3# TS2# 13/ ) MAI 

END MAT 


321C 

3220 

3230 

32A0 

3250 

3260 

327C 

3280 

3290 

3300 

33X0 

3320 

3330 

33AC 

3350 

336C 

3370 

3380 

3390 

3A00 

3A10 

3A20 

3430 

3440 

3450 

3460 

3470 

3400 

349C 

3900 


SUBROUTINE COST ( IGMO# TOTAL# AN I! # PNALTY# CLOSE # THRESH ) COS 

COMMON NTRAJ#NMAP,KSEG#XM(15#13 )# YM (15# 13 1 1 ZM( 1 5# 13 )# SPOS ITn5#25CCCS 
1# 2)#XD(15)# XF(1 5)#YCa5)#YFa5)#ZO(l5 )#ZFa5)iNPOSlTSU5) #PH0(573)CCS 
S# ISE61# IPRCF# MAXIT#PENCRIT# AWE IT# TV IOLA#ALWP# PEOPLE #NTP OPT# TMPLONT CCS 
S573) #OPD( 15)# ITG# ICPT# IT#ZERGL0B#NEWTPAJ#ARRAY(573#10) COS 

INTEGER QRO CCS 

LOGICAL ZEPGLae#NEWTRAJ CCS 

CGMMCN /CORNER/ ANGLE ( 15 # 12 ) # POST T ( 1 S # 250# 3 ) # XCENTR ( 1 5# 12 ) # YCENTR (C CS 

51 5# 12) #ZCENTRa 5#l2)#xn< 15# 1 2) # YT 1 (1 5# 12 ) # ZTl (1 5 # 1 2 ) # XT2 (1 5# 12 ) # YCPS 

ST2I 15#l2)#ZT2(15#12)#NPCSn(15)#0IS(15) CCS 

COMMON /PRINT/ SXCENTR (15 #1 2) #SYC EN TR (15# 12 ) # SXTl ( 15#12 )# S YTX (15# ICOS 

52 )# $2T1 U5#12)#SXT2U5#12)# 5YT2 (1 5 # 12 ) #SZT2 ( 15# 12 )# 5ANGLE (1 5# 12 ) CCS 

TCTAL-0. CCS 

AMI. TOTAL CCS 

CtCSE-ANII COS 

THRESH-CLOSE CCS 

PNALTYMMPESH CCS 


CALL DYNAPID (PNALTY) CDS 

CCS 

CCS 

. CCS 

SUBTP. NOISE EFFECTS OF TPAJ. BEING OPT. . C0< 

DLFING PPEV. CALL TO COST • •CCS 

. COS 

CCS 

CCS 

IF (.NnT.ZERGLCB) CALL ADOR (-1) CCS 

IF (PNALTY.lt. PENCPIT) GO TO 10 CCS 

TOTAL • PNALTY CCS 

IF (MAXIT.NE.O) FETURN COS 


10 

20 

30 

4C 

50 

60 

70 

80 

90 

IOC 

110 

12C 

13C 

160 

150 

160 

170 

IRO 

ICC 

2CC 

210 

22c 

23C 

2*0 

250 

260 

27C 

2RC 

29C 

3CC 


172 


10. 

cm NOlStL (THPeSH) 

cos 

31C 


CALI Nil <ANin 

cos 

3.20 


IF INTRAJ.ec.l) GO TO 20 

cos 

330 


CALL CRC5CVR ( NPCS IT# POSITfClCSE} 

cos 

3A0 

^0 

TCTAl • ANlI^THReSHfClOSEAPNALTY 

CCS 

350 


IF nOPAO.EO.n PFTIPN 

CCS 

3EC 


or 5C M • UNTRAJ 

COS 

370 


1 • DRonn 

CCS 

3F0 


smruu • xTui#u 

CCS 

390 


sKT2(i#n.« xT2n#n 

CCS 

AOC 


SYTKWn ■ YUn#L> 

CCS 

AlC 


SYT2(I#1> • YT2n#V) 

CCS 

A20 


szTin#n ■ zTin#i) 

COS 

A30 


S2T2U#l) - tlZiUX) 

CCS 

4AC 


DO 30 J ■ 2#NSEG 

CCS 

A5C 


SXCENTR(I#J-U • XCkNTPn#J-U 

CCS 

A6C 


SYCENTRMiJ-n • YCENTRn#J-U 

cos 

470 


SXTIUfJ) ■ >Tia#J) 

CCS 

A8C 


SYT1(I#J) ■ YT1U#J) 

CCS 

A90 


SXT2( I# J> • XT2(I#J) 

COS 

SCO 


SYT2ri#J) ■ Yr2U#J» 

COS 

510 


S2TK 1# J) - ZTin#4) 

CCS 

520 


SZTZMrJ) • 2T2U#J) 

COS 

530 


SANGlEfliJ-ll • ANGLEdfJ^n 

cos 

540 

30 

CONTINUE 

CDS 

550 


Npcsns a ) • NPcsmx) 

COS 

560 


IPCST « NPOSITS(I) 

COS 

570 


DC AO J ■ l#IPQST 

CDS 

5flO 


£POSITn#J#U • PCSITUiJ#!) 

cos 

59C 


SPC5rT(I#J#2) ■ P0SIT(l#J#2) 

cos 

6CC 


5PCSIT(I#J#3) • PCSITa#J#3) 

CCS 

610 


ccntinue 

cos 

620 

50 

continue 

COS 

630 


ftgrUFN 

CCS 

640 


end 

cos 

650 


SUPRPUTINE COSTl ( KHAD»N,F*X*ANIIii PN*LTY»ClCSE»THPeSH»NOtC 1 

CPI 

IC 

CC^PCN NTRAJ,NfAF,NSEGrXH(l5»13lfyP(i£>13l»ZP(l5»l3».SPtSlTa5» 

250CP1 

20 

1. 3)»XQOS»*XF(HI>YC(15I.YFaS(»ZO(15l*ZFn5)*NPQSITS(I5I.PHO(573>COI 

30 

Sr ISEGlf tPROFrMAXITrPENCPITr AWE! TrTViei A* ALWPr PEOPLE rSTRCPTrYnPlONICCl 

40 

1573 liOPO ns )#1TG# KPT/ IT# ZERGLOB# NEVTPAJ^ ARRAY (5 7 3# 10) 

COl 

50 

^ INTEGER CPO 

CDl 

60 

logical zcrglcb#sevtraj 

CCl 

7C 

OIMENSICN MNOECil) 

CCl 

PC 

NVAP • Z 

COl 

90 

IF (iPPOF.ECa) NVAR . 3 

COl 

IOC 

r • 1 

COl 

110 

NSEGl • NSEG-l 

CCl 

120 

OC IC J • l#NSEGl 

CCl 

130 

DC 10 K • I# NVAP 

cri 

140 

L • (Hl)RNVAR*NSKlf( J-l )*NVARtK 

cri 

150 

IF (K.EQ.l) XN(CRD(n#J^l) • X(lfl) 

cri 

160 

IF (R.EO.Z) YN<aRD(n#Jtl) • X(L#1) 

CCl 

170 

IF (K.F0.3) ZNlOFDi n# Jd ) • X(LaI) 

CCl 

18C 

10 CONTINUE 

CCl 

ICO 

CALL COST nGRAC#F#AMI#PNALTY#CLGSEfTNPE5H) 

cn 

200 

RETURN 

COl 

210 

END 

CCl 

220 



FUNCTION AACOS (AfB#C) 

AAC 

1C 


X ■ (A**2^09*2-C^*2)/(2.»A9B) 

AAC 

20 


IF (ABSIXI.IT.1.1) GC TO 10 

AAC 

30 


WRITE (6f9ClC) XfAfEfC 

AAC 

90 


STOP 

AAC 

50 

1C 

IP (X.GT.l.) X • 1. 

AAC 

60 


IP (X*LT.-l.) X ■ -1. 

AAC 

70 


AACOS - ACOS(X) 

AAC 

PO 


RETURN 

AAC 

90 


’ ^ 

AAC 

IOC 

9010 

FOKfiAT <29H TitOUBLE IN AiCOS> ttirtpC ■ >4(1PE16.9> 3)0) 

AAC 

no 


END 

AAC 

120 


SUBR0UTIN6 dCNlT F#4MI»PNAL TYiCLCSE »THRESH^NDEC) 

CCHHON NTFAJ#N?*4P#NSEG#)iHa5#l3)#Yf'( 15^ 131 >2rU5^ 13 )#SPOSn(l5»25Cf*CN 
»#3)iXC(15.)#XEa5)#YOU5>#YFn5) rZC(l5> # ZP (15) ^NPDSITS ( 15) #RH0<573)P0N 
»# ISEG1» XPPOP#HAXIT#PENCFIT#AWEIT#TVIQU#AtWP#PEOPtE#t^TPCPT#THPLON (HaN 
I573)#CR0( 15)»ITG#ICPT»IT»ZEPGICB^NEWTRAJ#AR®4Y(573#10) HCN 

INTEGER ORO * )^CN 

tCGICAt 2ERGLa8#NEWTRAJ ^CN 

COMMON /PRINT/ SXCE N TR ( 15^ 1 2 ) # S YC ENTP <15r 12 ) # SXT 1 (15^ 1 2 )# SYTUl 5 # I MCN 
S2)#S2Tl(ie# 12)>SXT2(l5#X2)#SYT2a5# 12 »#SZ T2( 15# 12 SINGLE ( 15# 12) PTN 


1C 


MCN 

MCN 

MCN 

MON 

MON 

MON 

MCN 

MCN 

MCN 

MCN 

MCN 

MCN 

MCN 

MCN 

MON 

MCN 

MCN 

MCN 

MCN 


OIMENStCN X(NOeC>l) 

PI « ATAN(1.)^A. 

NVAR • 2 

IF (IPPCF.EO.U NVAP • 3 
I ■ I 

NS661 • NSEG-l 
DO 10 J • 1/NSEGl 
OC 10 K • liNVAP 

L • (I-l)*NVAR'*NSfGl4(J-l)*NViPtK 
IP (K.EO.l) XMlOROm# J^l) • >(L#1) 

IF (K.E0.2) YM(CRDn)# J^l) • X(trl) 

IP (K.E0.3) ZM(CRC(D# J^l ) • >U#1) 

CCNTINOE 

WRITE (6#901C) lT#P#ANn#PNAlTY#Cia5E#THRESH#TVrOLA# ALWP 
MAX - I 

IF (ZERGLCB) MAX • NTPAJ 
DO 3C I ■ If MAX 

WRITE (6f 9020 CRD ( DfXMlCRO ( I) # 1 ) f YM ( CRD ( I) # 1 )# ZM( OR D( I )# 1) 

CO 20 J ■ 2#NSE6 

WRITE (6f9C3C) XM ( CRO ( I )f J )# YM ( CPO (I ) #J )# ZM( ORO < I) f J ) fSXCENTR ( MCN 
CRDd ) f J-l) fSYCEFTF(CRO(Uf J-l) fSXTUORDd )# J ) f SYT I ( GPO d )f J ) # MCN 
5ZT1(0R0( Df J)f SXT2(CR0(nf J) f S YT2 (ORO ( I) f J ) f SZT2 ( ORO ( I) # J ) # SAMCN 
NGLE(CROd)f J-1) AieCWPI MCN 

CONTINUE MCN 

WRITE (6V9CVC) XM(CRO(l)f NSEG^l »f YM(CRO( I)f NSIGtllf ZMIORDd )fNSEMCN 
^♦l) MCK 

NPCSITI . NPOSITS(CPCdn MCN 

CONTINUE MCN 

WRITE (0f9C5O) (IFIXARFAYtKf 3) )f ARRAY(Kf A)fK-IfNMAp) MCN 

ZEFGLC0 • .FALSE. MON 

RETUPN MCN 

MCN 

9C1C FORMAT (lHlf9Xf IIHITERATICM f I 2# / f 12 Xf 12HTCTAL GGSTt . IPEl 6 .C, / , if fN 
12Xf 17HANNCYANCE (NII)J f IPElt .9 ./# 1 2Xf 31HPENAL TY ON CYNRMIC CCNSTRMCN 
!AINT» f lPei6.9,/,l2X,3<*NPENALTY ON SEPARATING rCNSTRAlNTl flPEU.PMCN 
If /f l2Xf ZFHPENALTY ON THRESHCLD NOISE* # IP El t . «f / f 12X f pTHP ESP VIOIAMCN 
ITICN • ”f tl6.9f /fi2Xf>‘LWP • ♦»fElE.9f/) MCN 

9020 FORMAT ('.2Xf»PLICHT PATH NOi « f T2 # /f T20f »‘CdPN6 P PTS . ”f TAOf »»CENTER MCN 


20 


3C 


1C 
20 
30 
90 
50 
EO 
70 
PC 
90 
100 
110 
120 
130 
190 
15C 
lEO 
170 
IMO 
190 
20C 
21 C 
22C 
23C 
290 
250 
260 
27C 
2PC 
290 
ICO 
31C 
320 
330 
39C 
350 
360 
37C 
3PC 
3 90 
9CC 
91C 
A20 
93C 
99 C 
95C 
9E0 
970 


174 




• OF CJFCU»’#T7^i»»TANGFNnAL PTS. »*# T llCf »'AKGLfc (OECiEE PPCN fEC 

i7.Cr*> »#F7*C#»N«#f7.C#»«)«) ^^0 

903C FC«K3T n4^^»M»*>F7,Cr«#»*FE7.0fP»^rF7.0i»’>»»»>’(^»F7.0#»'#«rF7.0#.w>«#2hCN SCO 

M••^•#f7.0#»»#»^F7.C#*♦p»»F7.0#»'MM» 3X^FE«n HCW SIC 

9C4C FORMAT a4l^^«<«^P7.0^»«#M#F7.0#H,M,^7,C,«)«> PfK f2C 

9C5C FC«rA7 <57f lH^$)t#F9.5) 530 

^NO 940 


SUEFCUTINE PESULT C X#N^Ff ANH,PMlTYf CLO$E ^THFESHf NOEC ) Pf$ 

COPMON NTFAJfNP4F#NSEG#PH<15#l3J»YMU5fl3>#2H<I5fl3l»5Pa5im5»250PfS 

• O)r>0(15)iPFa5)#YC(15)#yFn5»#ZaM5)#ZFa5)>NP0SIT5IU>#RH0(573 >PE$ 

If lSECl# IFP0F#MAXn^FENCPU#UEITrTVJ0lA#4tWPfPeCPlE#NTPCPTf THftON(PES 
1573WCRD(15)#ITG^ICFTf IT#ZEPGLC0#KEWTPAJ»4PP4y<573#lC) fIS 

INTEGEP QPO 

tOGlCIl 7ERGLCft#hEWTRAJ 

CCJ-PON /PftlNTV $XCeKTFa5#t2)rSyCENTR(l5#12)#SXTUl5rl2^syTia5APE$ 
l2)#$ZTlUf»l2) » 5X12(15/12 )#SYT2 (15 » 12)# SZTZaSaZJfSANGU U5» 12 J PH 
OIPENSION K(NDEC#1) RE5 

PI • AT/NM.)**. PH 

NV4P • 2 , PES 

IF (IPPOF.EC.l) NVAP ■ 3 PES 

I . 1 RES 

NSEGl ■ NSEG-l PES 

00 10 J • 1#NSEG1 PES 

DQ 10 H ■ If NVAP PE5 

I • n-n4NVAP*N$ECl*(J-n*NVAP4K PES 

IF (K.EC.U XP(CFC( n# Jtl) ■ Mtfl) PES 

IF (K.E0.2I YP(CPC(n#J^U • K(L#1) PES 

IF (K,EC# 3 ) ZP(0FD(n#J4n • X(Lfl) PES 

10 CONTINUE PfS 

WRITE (6#9010) ITfFf ANIIfPNALTYfCLQSE# THPESHf TVIOUfALWP PES 

OC 3C I • IfNTPAJ PES 

ViPITE (6f902C> OFr(I)fXfi(CRO(Uf l)f yM(aR0(nfl)fZH(OPD(X)fn PES 
OC 2C J > 2fNSEG PES 

WRITE (6# 0030) yP(OPO( U# J) # YP( OPOCI ) f J ) f ZP( OPO (I ) # J) # SYC ENTR (PES 
i CPOd lf J-U fSYCENTF(OFO(I)f J-1) fSKTKOROU )f J ) # SYTl ( ORO (1 ) r J ) # P ES 
SZTK CPO( I)f J l#SXT2(CP0(n» J)f SYT2(aPD( I) #J ) fSZT2 (0P0( 11# J )» SAP ES 
NGl.£(OPO( I) f J-1)^160*/PI PES 

CONTINUE PES 

WRITE ( 6# 9040) XP ( CRO ( D# NS EG^ U # YK (QPD ( t ) #NSEG41 ) # ZP( QRO ( I )# NS EPFS 
Gfl) PES 

WPITE (6#9050) PES 

NPOSITI ■ NPCSITSUPO(I)) RES 

VPITE (fc#9C60) (tSPDSmOPO(I)f J#K)fK-l#3)#J*l,NP0SiTI) PES 

CCNTINUE PFS 

PETUPN PES 

: PES 

9010 FCFPAT nHi#9Xi llHITtPATICNI # 1 2# /# 12X# 12HT0TAL COS T I # IP El 6. 9# / # 1 P E5 
I2X#17HANN0YANCE (NID^ # IPE16#9#/ # 12X# 31HP6NALTY CN DYNAMIC CCNSTPPES 
fAiNTi f IPE16.9> /,12X#34HPENALTY ON SEPAPATING CONSTPAINT» /1PE1E.9PES 
t. /ilZXf 28HPENALTY CN THPFSHCLO NOISE* » 1 P E Ifc • 9, / r 12X , »THP6 SH VlCtAPES 
ST ION • «> El6r9» /, UX/HWP • '*#616. 9#/) PES 

9C2C FCPPAT U2Xf»»FUGHT PATH NO* «# 12 W # T20 # «CORNE P PTS . »'# T40# »’CBNTER PES 

• OF CrPCLE»’#T74,MTANGENTlAl P TS . *’# T1 1C /*ANGU (DEGREE )«#/# 14X# »M ’^iFPES 

7.0f*Vf »# F7.0# f )»M PFS 

90 3C FO^PAT (14X#<M»»# F7.0 # «# ♦’# F7. Of »*#« I F7.0^« )”r« (‘VF7. Of **#»'# F7*0# ♦*)»»# 2PFS 


20 


30 


i( »•( »»f F7. Of M^u,F7.0/^«fF7* Of ”)•')# 3XfFF*l) 

90 40 FOR PAT (14Xf 'M“fF7fcCf ”f »’f F7.0f»»f*SF7.Cf «)»') 
9050 FCFPAT <//#12Kf3thFlIGHT TMJECTOPY X> Y# Z, 
$4X»lHYfl4X,lHZf / ) 

9C eC FOFHAT ( 1 4Xf 1 PE I C . 3i t X# El 0. 3# 5Xf E 1 C • 2 f 5X ) 


PES 
PFS 

IN HETEPSf /f 19Xf IHXf IPES 
PES 
PFS 


10 

20 

30 

40 

50 

$C 

TO 

PO 

90 

100 

no 

12C 
130 
14C 
150 
IfC 
170 
IPO 
190 
200 
210 
220 
230 
240 
250 
260 
270 
2P0 
290 
300 
31C 
320 
330 
340 
35C 
36C 
370 
360 
3*0 
400 
410 
42C 
‘3C 
440 
450 
460 
4 70 
4PC 
49C 
500 
510 
52C 
53C 





o o o o r> o o r>o o o o r» rji o o o o o o or> o o o r»r» o o 


175 


SNO 


SiaRObTINC CFOSCVIi 


T>IS IS A OUCNY (CUnN£» REPLACING THE ORIGINAL CROSOVR# 
WHICH DOS NCT WORK fC8 TRAJS THAT OCNT SHARE A COPHON 
SeCHENT OF THE K-ARIS 


RETURN 

EKO 


SUBROUTINE ONEWTDN (N»XNOW,DELT*X* STORCHG,ND£C I 
CQHHQN NrRAJ.NNAR.NSE6»XHaJ»13l>YNU5»13 (.IHtl! 


»573).0R0a5l#ITG»ICPT.IT»Z£RGL0B.NeWTRAJ»ARRAY{573»10> 


THIS OPIIHIZATICN EHPLOYS SELF-SCALING# RESTARTING. 
OLASl-NEVTON METHOD. 

REFtRENCEi O.G. LUENEERGER INTRO. TC LINEAR AND NONLINEAR 

PROGRAHHINGI P.20A. SEC.R.5 

hAXITf HAXIMUH NLHEER OF ITERATIONS ALLOkEO 

STCPCHGl STEP IF PERCENTAGE CHANGE IN SUCCESSIVE COSTS IS 

LESS THAN THIS VALUE 

N« CIMENSICN OF THE UNKNOWN X 

XNCWi PRESENT OR INITIAL VALUE OF UNKNOWN X 


OIMENSICN XNOWINCEC.U. DELTAX(NOEC) 

DIKENSICN GNCWnC.n. GNEXTI30.1). P(3C.l). 0(30.11. P 0 (I, 1 ) 
OlHENSICN OSO(l.l). PP(30.30». S00SO0.30). S(30.30» 

^ 0IPEN5ICN XTEMPOO.n. Pm.301. QS(1.30). Son.3C», om.30> 
OIMENSICN SCC(30.30J 
OIMeNSICN 0(30,11 
IT • 0 

Call COSTI (O.N.FNCW.XNCW.ANn.PNALTY. close. thresh.noeci 
IF (PNALTY.CE.PENCRII) go TC 10 

call HQN'IT (XNCW.N.FNCU.AMl.PNALTY.CLCSE.rHPeSH.NOEO 
IF ( MAX IT. £0.0 I RETURN 
GC TC ?0 

10 WRITE (6.9C3Ct PNALTY.AMI 
CALL CORNER (XNCW.N.NOECI 
IF (MAXIT.EO.O) PETLRN 


STEP II SET S • ICENTITY MATRIX AND CALCULATE GRADIENT 6 




540 


cto 

10 


CPC 

2C 

t • • 

fPO 

30 


CFO 



CPO 

?c 


CPV 

60 


CPO 

70 

* 

CFP 

NO 

> • • 

fpn 

90 


CPV 

100 


CPO 

no 


CPO 

120 


OM 

10 

-25CON6 

2C 

m 

)ONf 

30 

.ON(ONE 

40 


ONF 

50 


ONF 

60 

IS t 

ONF 

7C 

• 

o^e 

8C 

■ 

QH 

<i0 

• 

ONE 

IOC 

• 

ONF 

no 


ONF 

120 

• 

OF)F 

130 

■ 

ONF 

no 

• 

ONF 

150 

• 

CNF 

160 

• 

ONF 

170 

• 

ONF 

180 

‘ • * 

CNF 

1«C 


OFF 

2CC 


CNF 

21C 


ONF 

220 


ONF 

230 


ONF 

240 


ONF 

250 


ONF 

260 


OFF 

270 


ONF 

2F0 


ONF 

2«0 


ONF 

300 


ONF 

310 


ONF 

320 


CFF 

330 


ONF 

340 


ONF 

350 


ONF 

360 

• » 

CNF 

370 


ONF 

I8C 

• 

ONF 

3<50 

• 

ONF 

4CC 

«. • 

ONP 

410 


ONF 

42C 


ONF 

430 


2C OC 3C I • l.N 
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t CC 3C J * l#K 

sn^j> • 0. 

IF n.fo.j) sa#j) • i« 

CONTINbE 

CAU Fa»40 (Nr FKCU^KKOk#GNCb#OetTAXfHrECi 


STEP Zt SET 0 i -SG 


AC C4U HPIY (NrNf WSfGNCWrOfNOeCrNOECrn 
DO 50 I • 1#N 
OUrU • -OU/1) 

5t CCNTINl’E 


STEP 3t MOST IMPCPTANTI 

UNt SEAPCH AUNG D TC FIND AFA THAT SATISFIES PO>C 
IF COST FUNCTION IS ANALYTIC# SLIGHT MODIFICATICN IS 
NELDEO (SIMFllFICATIdNMOFTIONAU 


CNF 

K • C ONE 

60 ^ • K41 one 

CALL LINFSCH (K#N#FNQk#XN0W#0#AFA2# XTEMP^FSMALL# AKn#FNALTY#CLOSE#ONE 
ITHfiFSKiNDEC)# PETUPNS 1220) QN| 

IF nT^l.LT.MAKIT) GO TO 70 ONE 

n • tin ONE 

WPITE (6# 9050) ONE 

IF (PNALTY*GE*PENCRIT) CC TC 15C ONE 

IF (FSMALL.LEtFNOW) CALL NONIT (XNCW# N# fS MALL# ANM#PNALTY#CLOS|#THONE 
IPE$H#N0EC) ONE 

FETUFN ONE 

70 CALI FGRaO (N#FSMALL#XTEtP#CNEXT#OELTAX#NOeC) ONE 

00 PC I » X#N ONE 

on#n • cNEXTn#i)-GNrw(i#i) one 

P(I#n ■ AFA2*0(1#1) CNE 

EO CCNTINUe ONE 

CALL TPNSPOS (NOEC# 1# P#PT) ONE 

CALL HPLY Jl#N#l#PTfC#PO#I#NDEC#n ONE 

^ IF (POUfD.CTiO.) GC TO UC ONE 

IF (KtLE, 1) GO TO ONE 

VPITE (6# 9010) ONE 

PEttPN ONE 

9C FNCk • FSMALL ONE 

CO KC I • 1#N ONE 

XNQW(I»l) • XTEMPU#!) ONE 

100 CONTINUE ONE 

CC TC 60 ONE 

UC CO 120 1 ■ I#N ONE 

XNCWdfU ■ XTEMPU#!) ONE 

12C CCNTINUE ONE 

IT « nn ONE 

PPCNT > A8SnFSMAU-FNDV)/FN0W) ONE 

IF (PRCNT.GE.STCPCHG) GC TO lAO ONF 

WRITE I6#904,C) PRCMrSTCPCMG ONE 

IF (PNALTY.GE.PfcNCFlT) CC TO 13C ONE 

CALL MONIT t XNOW# N# FSMALL# A M I# PNAL T Y# CLO SE #THP6 SH ^ NOEC ) ONE 

ONE 

PETUPN ONE 


PETUPN 


rn r» r> t*> o rv o o o n r% r» o o o 
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13C VRiie ie#9oao> PfAirt 

CAU CQRNffl UNOk»N#NDIC) 

AeTUPM 

14C IF lIT.Lr.f^AFm GO TO UO 
yffire (6»9050» 

IF (PNALTY.CF.F6NCFIT) CO TO 150 

CALL nOhU <XN0^#N»FS>^AU#AKII#PNAlTY|Cl05E#THPeSM#N0CO 
RFTUF^ 

150 WFITI U#9020) P^AiTt 

CAU CONNER (XNOWfH^NDlO 
PiTUPN 

160 IF (PNAlTY.Ce«Pl6CRIT> CO TC I7C 

CAU ♦^ONIT (XNOW#N#FSHAU»ANII#PNAlTY#ClOSE#T^<PeSF>NOtCI 
GC TO IBO 

170 WFITE (A#9C30> PKALTY#ANM 
CALL CORNfP rxK0W,N#N0|C) 

IPO CONTINUE . 


* * 


1 • • • • 




STE? 

41 

If "IT" IS INTEOEK MUlTIflE CF N 00 TO STEP II 




IF 

NOT> UPDfTE S 



IS 

HITH 

INTEGER multiple CP HT 


• • < 


» • • • • 




IF |rFLQATriTI/FLCAT<N)).Ne.FLOAT(IT/NI) GO TO 190 
FNOW • F5NAU 
GC TC 20 


UPOATE HATRIX S| GC TO StEP 2 


19C CAU NPLY <N»N# l»SfCrS0#N0EC>N0EC> U 
CALL TftNSPCS (NCEC # lrO#OT ) 

CALL PPLY tN,l>N#SC#OTfSOO#NOEC>l#NOEC7 
CALL PPLY <N#NfN#5C0f S#S0CS#N0EC#N0ECfN0EO 
CALL NPLY (lFN#N,0T#S^0S#UN0EC#N0Ea 
call fPLY U#N#l#0$^0j OSQ»1#NOFC? U 
call HPLY (N»l»N#P»PT#PP»NOECUrND€C) 

200 I • UN 

• 1»N 

( 5 1 If J ) - S C 0 S ( Ip J » / 0 SO ( 1 • 1 ») • < P 0 (1 pX » / 0 S 0 < I p 1 ) 


CO 20C J 
SMf4) • 

% PCU#1) 
eOC CONTINUE 
PKVf • FSPALL 
DC 210 I • 1#N 

GNCW<IU) • GNEXTTWn 
?IQ CCNTINUE 
GO TC 4C 

?ZQ WMTE <6#9C6C> CFOil) 

peruPN 

"9010 FCPFAT |/,lOXf«AFTLP 2 CALLS TP LINESCHt CANT fINO TPAJ. 

SOkFP COST. RETURN TC PK*^FdeN.«) 

9C2C FGPrtAT <5X>A0H|N RESULT! OVNAMIC CONSTRAINTS ViaLATiaNf/ 
lALtY •»1PEU«9» 

90 30 FCLFAT (5X#39HIN NCNIT! OYNANlC CQNSTFAIMS VI0LAT1CN»/^ 
SLTY •fiP£U»9f/#lCXf*»ANir • «»ei6-9r'’MF • C«# THEN ANl 
IT PEEN COPPUTEO IN THIS ITEFATICN# OLE TO LARGE PNALTY**) 


m 
ONI 
CNF 

m 

ONE 
ONE 
ONE 
ONI 
CNF 
ONE 
ONE 
ONE 
ONE 
ONE 
ONF 
ONE 
ONE 
ONE 
CNF 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 

••••••••• ONE 

• ONE 

• ONE 

• ONE 

••••••••• ONE 

ONE 

ONE 

ONF 

ONE 

ONE 

ONE 

CNF 

ONE 

ONE 

ONE 

MPPM# J1/0NE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
ONE 
CNF 

WITH’*#*’ LONE 
ONE 

#10X# OHPENCNF 
ONF 

1CX#OHPENAONE 

I HAS»*f«N0ONE 
ONE 


lOPC 

1090 

IICO 

1110 

U20 

1150 

1140 

1150 

1160 

1170 

llEO 

1190 

1200 

121 c 

1220 

1230 

1240 

125C 

126C 

1270 

12PC 

1290 

13CP 

me 

1320 
1330 
13^0 
1350 
1360 
1370 
1360 
1390 
14 CO 
lAlC 
1A20 
1430 
1440 
1450 
1460 
1470 
146C 
1490 
1500 
1510 
1520 
1530 
154C 
1550 
IfEC 
1570 
1560 
15<50 
1600 
1610 
1620 
1630 
16AC 
1650 
166.0 
1670 
16P0 
1690 
1700 
1710 
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9C4C PCPfAT CHIK6E IN SUCCfSSZVf COSTS *•#1^110.3# jlCSS"ONe 


if« TNAN STOP CPlTiPION **,jp||0i3l 0N$ 

90SC PCfJ^AT U)T^29HHAUNtH mftATICN SIT PlACMeOJ * ONf 

9060 rCPP*T if 9^ •.••• OPTIMOnCN PAILS fCP TPAJ» ON THI5"#»* CQH 

UOPAl ITkiATIONH#//) CM 

IND ONI 


SUePCUriNt NPLY U^K#N4*#lgiCnoeC#»CECfN05C> Hfl 

C NPL 

C NPL 

C 9 . NPl 

C . CALCUlATt HAT«n NUtIPUCAlION C • 46 # NPl 

C . . NPt 

C 9^9 9 9 9 9 9 . 9 9 9 9 9.9 9 ( 9 . 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 fPl 

C NPl 

OlPiNSICN ML06C9>^CPC)9 0 (NOSC# NOPC C(tO€C#NOeC) PPL 

DC 1C 1 • 1#L PPL 

OC 10 4 • WN PPL 

cruj> • 0 . PPL 

OG 10 K . IrH PPL 

C(l#j) • cu# PPL 

1C CGMIM/L PPL 

^0TIPN PPL 

m PPL 


SLBPOITINE TPNSPCS {HiNfA/P* TPM 

C TPM 

9 . . 9 . 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 . 9 99 9 9 9 9 . 9 9 TPN 

C . 9 TPM 

C . TP^N'POSB OF rtATFiy A IS RE TURMEO IN PAIR I X 0 . TPN 

C 9 9 TPN 

C 9 9 9 . 9 ......... TPM 

C TPN 

OinENSlON A(p9N)f B<K#P) TPN 

OC 1C 1 ■ 1#M TPN 

OC 10 J • l.N TPN 

a(jfi) ■ Au#*n TPM 

10 CCNTINUE TPN 

•vPeUPN TPN 

EFC TPM 


SVEPOUTlNt FGRAO ( N. f . X. G>OEL TA P#NOEC ) FGP 

C FGP 

C 9 9 ••• . 9 999 . •.9. 9 •••••.. 99. 9.9 ••.• 9 9 99 9 9 9 9 9. 9. 9 9 . 9 9 FCP 

C 9 9 FCP 

C . CALCULATE GRADIENT CF COST F WITH RESPECT TO UNKNCWN X . FGP 

C 9 9 f fP 

C 9. 9 •••9 .99 ••••••.•• 9 9 9. 9. 9 9 . 9 9 9 9 . 9 9 . 9 •»9. 9 9 . 9 9 9 . 9 9 9 •.••.•• 99 •••• 9 99 9 . FCP 

C FGP 

COMMON /ENAEIE/ lENABU FCP 

LCGICAL ZfeNA 0 lE FCP 

OIMENSICN XCNDECM^ G(NOECrl)f D EL TA X (NCEC 9 1 ) FCP 

00 1C I • l.N FTP 

cn.l) • 0 . FG» 

IF ( .NOT.ZENABLE.ANO.MOOn^ 31 .EO.OJ GO TO 10 FC» 


17?0 

1730 

1740 

I7SC 

1760 

1770 


to 

?C 

30 

40 

SO 

60 

70 

#0 

«0 

ICO 

MC 

170 

13C 

140 

ISO 

160 

170 


1C 

2 C 

30 

4C 

?C 

60 

70 

eo 

PC 

toe 

1)0 

120 

130 

lAC 

150 


IC 

20 

30 

AC 

50 

60 

70 

PC 

4C 

100 

no 

120 

130 

lAO 


• XM>ll«OELTAXItl 

CtU CC5U n»N*EFrK»*MI>PMLTY«ClCSE>THXESH>NOEC) 
Cd^ll • tfF-fl/CHT*XII) 

X(I»1) • XU»l»-OElT*xm 
10 CCNTIhUE 
PETOFN 
END 


SUnPCUTl»«E C*F* IN#FfO«N^GfX*D»NCEC> 


CXICUXTE CFXOIEHl CF F fclTH FESPECT TO *Fi 


OIHEKSICF X(NOEC>n> 0(NCEC»1)» XTEFPI30>1> 

00 10 t • »>N 

XTEPP(I»U • X(l»lM0*F**0n»n 
10 CONTINUE 

C*U, COSTl fl»N»fF»XTENP,PNlI,PN*lTY»CLOSE»THPESH»NCECl 
C ■ (FF-F)/0*F* 

PETUPN 

ENO 


SIPPCUTINE ERRCR (K»P|I5«CE> 
OIHENSICN HEJSX6E(1?>» IIEXTll?) 
OiTi ITEXT/IZ*** »V 

00 10 1 » 1,K 

IlEXTdl • flESSFCEU) 

WRITE (fc/90101 UTEXTin.l.l»K> 
STEP 


10 


9cie 


FCfPAT 

END 


t//dX>"ePPCFt «>12U0d/l 


ICtCSE. THRESH. NOECJ. RETURNS IKABORTI 
OIHENSIDN CtNCEC.ll. XTEHPI NQEC. U . XNDWCNOEC.l) 


CteiC FIT BY INniFLLY LEITING AF* • 0 ANO OAF* • 5. 


ONDRR • 0. 

CO 10 I • l.N 

ONORH • DNCRH*C(I>1JR*2 
10 CCMIHUE 

CNURN • SCRTIONCRNI 

IF (CNCRr.EO.O.) CALI ChfCK ( l.N.FNCW.XNOV.D.AFAO.FAFAO.f.AFAO.AFA 
1. FAFAl.GAFAl.Ll.LUl.UZ.NDEC) 


ret 

lU 

ret 

m 

F6» 

lit 

rcn 

IH 

rc» 

IGG 

ffip 

ICC 

m 

210 

eif 

to 


20 

1 OAr 

30 

1 Clf 

*C 

1 CAr 

30 

> GIF 

oc 

» GIF 

70 

GIF 

FO 

GIF 

«0 

GIF 

IOC 

GIF 

no 

GIF 

120 

GIF 

130 

GIF 

no 

GIF 

150 

GIF 

uo 

EIF 

10 

FF« 

20 

FFF 

30 

FM 

10 

£FF 

50 

EF» 

60 

m 

7C 

FFF 

90 

FFF 

90 

f»F 

100 

#UN 

1C 

ilN 

20 

UN 

30 

UN 

1C 

UH 

5C 

UN 

60 

LIN 

7C 

UN 

fO 

UN 

CO 

1 IN 

ICC 

UN 

lie 

UN 

120 

UN 

130 

UN 

no 

UN 

150 

lUN 

16C 

UN 

170 

UN 

leo 

UN 

190 




c 

c 


180 


• 

« 

NO^Hiirze D * 

• 

IIN 

UN 

UN 

UN 

ICC 

tit 

zzt 

z$o 

• • ♦ 


UM 

140 


DO 20 1 • l#N 

UN 

190 


CM»l) ■ 0(I»l)/0NCl!H 

UN 

700 

?c 

CCMINU6 

UN 

770 

OAfA • 5. 

UN 

7F0 


AFAO • Of 

UN 

7^0 


fkfto • FNOW 

UN 

100 


KK • 1 

UN 

no 


CAU CAFA |N>FAFA0>CAFA7CAFA0>FN0h>0*NDEC) 

UN 

370 


If fOAFAO.fOtO) CAU CHECK U#N#PNCJW#XNOViO# AFAO#rAKAO#CAKAC|AFAiU IN 

330 


iFAFAi»CAFAliUl»lU>l2ANDEC)A FETURKS (901 

UN 

140 


AFAl ■ AMJNmOCC»»-F*FAO/GAFAO) 

UN 

390 


IF (-FAFAC/OAFAO.tr.O) AFAl • AMAR U-ICCO. »-F4FA0FeAFA0 > 

UN 

no 

30 

OC AC 1 • l»N 

UN 

3T0 


xTiKp(un • )(KOwn#iMAFAifpavn 

UN 

3FC 

MC 

CONTINUE 

UN 

340 


CAU cosn a»N,FAFAl»XTEHPiANn#FNALTY,CLOSE»THRESN»NOEC» 

UN 

4C0 


IF (KK.LI.A) 60 TO JO 

UN 

410 


wPin (6»90A0) 

UN 

4?C 


neruKK 

UN 

430 

5C 

CALL OAFA (N.FAFAl,CAFA»6AFAl»)fT6FF»0>N0EC» 

UN 

44C 


If IAESKAFAC-AFAII/AFAI) .LT..001I CAU CHECF (2»N, FNOW»MNOW»C» AFAL IK 

490 


»0#KAMO#CAFAOiAMUrAFAl#CAfAUUUUUl#UZ#NOfC)r KETUKNS 

UN 

400 


U1 • CAF40A0AFAI«|,»!FAFAC=PAFAI)/(AFA0-AFAU 

UN 

170 


UUl • WI**2 -GAFaC*CAFA1 

UN 

470 


If aUl.LE.O.) IZ • 0« 

UN 

440 


IF (UUl. GT. 0.1 U2 • JOFTTUUl) 

UN 

9 CO 


AO • ANINKAFAO.AFAI) 

UN 

510 


A1 • AHAFKAFAC.AFAI) 

UN 

970 


GC • GAKAO 

UN 

930 


FO • FAFAO 

UN 

9*0 


GX • GAFAl 

UN 

990 


FX • MMl 

UN 

960 


If UFAO.lT.AFAX) GC TO (iC 

1 IN 

970 


GO • OAFAl 

UN 

950 


FC • FAFAl 

UN 

900 


G1 • CAFAO 

UN 

6CC 


FI . FAFAO 

UN 

fXO 

fcC 

IF IGl-G0«2.AU2.F6.C.t CAU CHECK <3 >K«FNOU|XKOU.O»AO«FO<GO >A l.Fl 

#UN 

670 


16l.Ul»Ul»U2»ND£C)» PeiLRNS (90 » 

UN 

630 


AFA2 • Al-(Al-AO)*(0UL2-U)/(0I-6CA2.»U2) 

UN 

640 


UN 
1 IN 

690 

660 


• 

UN 

670 


RESTRICT HAXIKlt AFA2 TO l.EA . 

UN 

670 


* 

UN 

UN 

UN 

640 
70 C 
no 




AFA2 • AMNl(AFA2.l.tAI 

UN 

770 


IF (AFAe.GT.O.) AFa 2 • APiNU AFAZ « 1 .£9 ) 

UN 

730 


IF (AFAZ.IT.C. ) AFA2 . AMAA l{«FAZ.-l. E9) 

UN 

740 


PC 70 I p 

UN 

790 


XIEHPd.ll ■ XKCWn.l)AAFA2ADn,U 

UN 

760 

?C 

CONTINUE 

UN 

77 C 


WRITE Ifc.9010) AC» Al. AFA2»0N0PP»FC*FI»GC#G1 

UN 

790 


WRITE (t.oOlOl (DM.l) .I*liN) 

UN 

740 


CAU COST! 1 0. N, FSHAU, XIEHP. ANTI . PNALTY. CLOSE. THPESH.NOECI 

UN 

500 


WRITE (t» 90201 FSHALL 

UN 

910 


IF tf SHALL. GE.FNCM GO TC SC 

UN 

970 


RETURN 

UN 

H30 


to AMI •» ArAl/ 2 » 9 AC 

KK • KKtl UN 990 

OC ?0 10 UK i90 

PttUPH KAICf^f UN «T0 

C UN «9C 

OOlC FCPNAT fUMMO • WrU^m*6#SfIrNAl • MIU*6^9)f#HArA? • *f#IU«0i9XUN l«0 

IfWONCIN • N|6l3.e#MHX#*’FO • N,|i3,Or5)t#Nt| • ••if l|,(6?r/# |?X^«00 •UN 900 

i ^•»U3*6|^X^»•0l • UN 930 

9020 fOFNAT a5X#«FS9AU • UK OfO 

9010 fCXNIT (tXiHO • *»#IOtlX#l9l9t2#2XH UN «30 

9090 f09NAT W/UOXit'UNISCH CANT UNO tONff COST AfTlI 1 CUfcIC fITS»*l UN 9AC 

I NO * UN 950 


5189CUTIM CHICK nCHICK#N# TN0N#KN0W,0rA0|f0iC0>Ai liOUUJ #UUl#U2CHf 


^•NDfC), RtTUFNS (KAB0RT1 

CHf 

OINfNUCK XNCWINOICUl# OINCeCil) 

CHf 

wpm uiior FNow»NMirKNcwM#iboriuiii"UH) 

CHf 

IF UCHlCK.IO.il RITURN KA0ORT 

CHf 

WRITf Ui20) A0^F0»CC#A1#F1#G1 

CHI 

IF UCHICK.IO.ai RETURN KABORT 

CHf 

W«m U#30) OliUUl#U2 

CNf 

RITURN KABCRT 

CHf 


CHf 


10 rOKNAt j2CK#99HTHrs 1$ SUSACUnM CHICK WHICH CIVtS AU TH«|15H INCHf 
IKCRMTICN IN SUAOUlNe UNi5CH,U /i10X/7Hf NOW * #lAf 16*9^ /># 37X# AfNf 
|MKN0Vii7yUH0i//>-(25KU2>3XilKll6.M9XiUEie.9i/V> CHf 

20 fCKNAT (//»20X»»^A0 • •Vlf»U6*9#3X#»^F0 • »♦# lll.9i 3X|NC0 ■ ••#m,9WCHf 
I2CX|»AI • H,616t9»5X#K|s| • «, IH.9#5K#»»U • ♦*#U6*9> CHf 

3C fOftHAT (//rfOXrIhUl • xm6.9#5)i#AHtUl • UKU0.9#5I^3MU2 • ^IPIUCHf 
3*9> CHf 

£NC CHf 


10 

fC 

ac 

AO 

30 

10 

7C 

fO 

00 


ICO 

110 

120 

130 


lAO 


13C 

uo 

170 

IfC 


StaPCUTlNE CD9NIA (X^N^NOfCI CCK 

CCNNCH KTRA^| NKAP|N«fGiXfm5U3 )#yim3# V3>#2Nn5# 13 SPOS IT 1 15V23CCCP 
1#3I #lcn3)iPH13bVCU5)^yF(15»#7OUn#ZP n5)#NP0SITSa5>#PH0(573)CCP 
IUC60lilPPCF#NAyiT#P£NCPlT#AWIITi!VIOUrAlWP#PfOPLf#NTtCPT>TKPtON<CCf 
l371h.OPDU3)rtT0flCPrin#2|PGt08iNeWTPAJ#APPAy<t73f 101 COP 

INTtOfR DRO CCR 

tCOICAl ZeRGlOBiMhTPAJ CCR 

*^OIK|NSION X(NOlC#n CCR 

HVAP • 2 COR 

if nPPCPilC.U KVAP ■ 3 CCR 

nc 1C I « l#NT»CPt CCR 

NSfGl ^ NSIG-l COR 

OC 10 1 • IfHSlGl CCR 

CO 10 K ■ liNVAR CCR 

I • a-M^KVAR*NSIGn (J-U*NVAR4K Cpp 

IP (K.io.u xN<cPD(n#j>n • x(t#n ccp 

IF (K,|0.2> YN*DFrU)#aM) ■ xarl) CCR 

IF (K. 10.31 ZK(CPO(n>M) • X(L#U CCR 

10 COMINUf , CCP 

OC 10 I • l^NTRAJ CCP 

• Pita U#9010) WXK(CRD(n>lj#YN(CR0riMl)U«(0»CM1»l) CCR 

OC 20 4 • 2rN5IG CCR 

VRITF U>902C1 lK(OPO<I)#JUYN(CRC(l)f J)f2N(aP0(I)f4) CCR 

20 CONTI NUC CCR 

wRlTt i 6#9C2C) XHfOPCt l)#NSlGKll#yN(CR0(l bNS6G^l)/2N(0R0n )iNSfCCR 
I G*U CCR 


IC 

20 

30 

AO 

50 

AO 

70 

BO 

90 

100 

110 

12C 

130 

1*0 

150 

UO 

170 

IfC 

190 

200 

210 

220 

230 

2AC 

2fO 

2tO 


182 


30“ CCMIKVf ^7C 

cot tn 

C m ftO 

9C1C mj**? fI2K#HFt|CHT Hlki NOI »,U#/lUK>«CCtMP I^TS.**# /# («>F0CCt 100 

cct no 

0010 fPHiT CPR IXQ 

iU CC» 130 


oYNiito m^un nw 

COHUCN NtMa#Hf*i?>K5fO#XHii5n3lrYl^(lf#13)#2Pri5a3)^$tOSnnij2100m 
l#3by0tl5l»)(r(l5)fYCCl§J#VP<l5>#roa5)p2r(in#NPQ5fT5a5>#PNq(5n)PYN 
•#!SfCl#IMaf#HiKIT#Frl^CPIt#*Vtl1jTVIOt*#liWP|f»IOPUE|WTtOtTrTHPtON(OYM 
l$7lbOf^OlllKr{TOllOr»T#11#Z£fOlOB#N£liTtJka#AII9|iYll73#10l OYH 

IhT6(5|t OFO OYK 

idoieii OYK 

CWCH / 0 YKAH 1 C 7 tACIUS#PW||Tl#?WEn 2 #ZWIItl»<SH^HIHn 5 bOHlf»A)f(lf > 0 YN 

|>||€ 4 I<P (jYH 

COH^'CN /COtNf^/ 4NCafcCl5a2)#PCnm5#230#3)#XCIHTta5#ie)#YCfHU<0YN 
nj>l?J#2C|NTt(Ua2J#XTUl3iU)#Yrm3#lzama5^UblT2a§#UWY0YN 
ftzniiizj r2T2(i j^ui^Nposrr (15) #015(151 oyn 

UNOTH OYN 

tENOTN(^?|#C^Pll#F) • $0ttnA-ei*t24(C-0)#n4(e-F)t*2) DYN 

PJ • AtANCl •)»<»• OYK 

OYK 


Pli-$t CAUUtATI PAPmmS At CdPNCP POINTS 


HAX • NTACPT 
IK aeftctcft) HAx • MflAj 
INCLUDE ALL TftAJS KCP THE 
00 EO II • I#NAX 
I • CAP (ID 
00 10 J • Z^NSEC 


ZEPOTH 0LC8AL ItEPATtON 


10 


YTZ(DJ) 
YTKIf J) 
XT2U#J) 
XTUDJ) 
ZTl IIiJ) 
nZIDJ ) 
ANGLEa#J- 


YCENTR a# J-1) 
KCEHTPn# 4 *-D 
CCNTiKUl 


0 # 

YT2(I#J) 
YUn# J) 
XT2(1#J) 

0 « 

c. 

n • XTDIiJ) 


ANT‘tE(l#J-n 

YCENTPUi4-D 


015(1) 

KTl(Dl) • XH(1#1) 

XT 2 (I#D • XM(l#ll 
YTDDl) • YHDrD 
YT2(l#l) - YN(Dl) 

ZUdrl) • ZH(|#1I 
ZT2CD1) • ZH(DD 
ISEC2 *p ISEClAl 
00 AO 4 • nECZfKSK 

A • lEKCTH( XHn#4-D#XH( I#4)#YHn#J-l)#Yrn#J )#0##0f ) 
e • L ekoth ( )h(i# j ) d n d 4 > ii# yh ( d 4 i # yh u # ja i » # c. # o. ) 

C • LEnGTHDH (D4^1)DN(D J^Di YHT r#4-l ) #YN ( I^4D )# ) 


OYN 

OYN 

OYN 

OYN 

OYN 

OYN 

oyn 

OYN 

OYN 

OYN 

0 >N 

OYN 

OYN 

OYN 

OYN 

OYN 

CYN 

OYN 

OYN 

OYN 

OYN 

OYN 

OYN 

DYN 

OYN 

CYN 

OYN 

OYN 

OYN 

OYN 

PYN 

OYN 

OYN 

OYN 

CYN 

DYN 

OYN 


10 

20 

1C 

4C 

SO 

00 

70 

PC 

•0 

ioo 

DO 

120 

ISC 

HO 

150 

100 

170 

no 

i«o 

200 

210 

220 

23C 

240 

250 

240 

270 

2 P 0 

200 

ICC 

sio 

320 

130 

340 

150 

360 

170 

160 

390 

4CC 

410 

420 

430 

440 

450 

A60 

470 

4P0 

400 

500 

510 

520 

530 
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C . » 00 THE TWO SEOWfhTS FCRh A STRAICHT LINE? . OVN 5A0 

C • •OVN ESO 

C 9«0 

C OTW 9T0 

IE <*BS(*a6-C(.6T.1.I go 10 20 OVN J#o 

SEC ■ L£NGTH(XT2(I.J*l).XMU*J>.YT2(l*J*l>.Y(Ml.J>.ZT2tl.J-n.0VN !4C 

i ZHU.JM OYH ECO 

YTUI.J) * YNI.JI OYN tic 

YT2(I»J> • XMd.Jj OYN 620 

YTia.il • YNa.il OYN 6)0 

nzil.jl • YNd.J) OYN 640 

2Tia.il • ZfMI.JI OYN 650 

2T2(I.il • 2Ml.i| OVN 660 

ZCENTEa.i-1) ■ C. OVN 6T0 

YCEhTRa.J-ll • 0. CYN 660 

YCENTEM.J-ll ■ YCENTSI l.i-ll OYN 690 

ALONG 1 "0. OYN TOC 

ANCLEa.j-ii ■ 0 . OYN no 

CO TO 30 * OYN 720 

20 A»Y • lENGTH(XI'a.J-ll.)(Na»Jl»V«a.J-ll.YNa.JI.O..O.I OYN 730 

BAY ■ LENGTHaNa.JI.AHII.iAU.YNa.JItYI'a. J41I.0..0.) OYN 740 

CAY • LENGTHtXNll. J-ll.ANa.J.l).Y«a.J-ll.YNa.J*ll>0..0.> OYN 790 

BETA • AACCSTAAY.BAY.CAYI OYH 76C 

ALCNCZ . R40IUS/TANTBETA/2. I OYN 77C 

PNALTY • PNAlTY»PUElTl*aANAXHC..(ALCN62-AAYIII»*2«ANAAll0..<0YN 760 

» A10NG2-BAYIII»*2 OVN 7«0 

IP (i.NE.ISEGZ) PNALTY* PNALTY»PWE IT 1*UNAXU0.. UlONGlAALCNGOYN iOO 

» Z-AAYII)»*Z OYN 6|0 

C OYN *20 

C OVN 830 

C . .OYN 640 

C . CALCLIATC TANGENTIAL POINTS. CINTEP OF RADIUS R AT J*TH CORNER . OYN 690 

C . , OVN *60 

C OYN 670 

C OYN 86C 

ATUI.JI • AHa.i|Y4L0NGZ*(XM(I.J-ll-A«a.JII/*XY OYN 600 

YTia.il ■ YWa.J )tALDNG2»(YM( I»i-1 j-Y«a.i) l/AXY OYN 900 

ATZn.J) . ANa.aitAL0NG2*(A«a»i*ll-A«a.JIIABAY OYN 910 

YTZa.JI ■ YNa.J|.AL0NG2A(YNn.J.J.)-YNa.J)l/BAY OYN 920 

2 Tld.il ■ ZN<I>il*ALCNGZ»( 2 Ma.i-ll-ZN(l.il)/AXY OYN 030 

2T2a.il ■ ZN(I.iHAL0NG2*(ZNn.i.ll-2na.ill/BXY OYN 940 

SEG • LENGTNIXT2a.i-ll.ATia.il.YT2a. i'll. YTia.il.2T2a.i>10VN 970 

t I.ZTlII.ill OYN 9(0 

0ISN2 • ALCNG2»C0S (BETA 12.1 OYN 970 

^ 0ISC2 • RAOIUS/SINIBETA/2.1 OVN 960 

A • (XTia.il.ATZa. ill/2. OYN 990 

Y • (YTia.iHYTZa.ill /2. OVN 1000 

ACENTRa.i-ll ■ AMa.ilYOISCZ* (A-Ana.il|/I)ISN2 OYN 1010 

YCPNTRa.i-ll . YMa.il»0ISC24(Y-YNa.ill/0ISN2 OYN 1020 

ALCNGl . ALCNGZ OYN 1030 

C . LENGTHlXna.il.ATZ(I,J),YTia»il.YT2a.il.O..O.) OYN 10*0 

ANCLEd.i-ll • AACOSIRAOIUS.RAOILS.DI OYN 1070 

SEG • SEG.ANGLEa.i-llYRAOIUS PVN 106C 

30 oisai • oisa)*stG dyn ioto 

*0 CONTINUE OYN 1060 

CO SO i • TSeGZ.NStG CYN 1090 

C » OYN HOC 

C OYN lllC 

C . .OVN 1120 

C . CALCILATE VEPTICAL SLOPE Op EACH SEGNENT AND . OYN 113C 

C . ACO penalty POR E/CH segment with VERTICAL SLOPE GREATER . OYN 11*0 

C . THAN GNAHAA OR LESS THAN GHAHIN DEGREES . OYN 1190 

C . . OVN 116C 

C OYN 1170 




104 


c 


• DYN n»e 

>n • lENSTtMXCn* J4lW)INM»J»»VH|I»J*U»Yf'n»J>»0.»0,» OYK 11»C 

VSICP6 • OYN uoe 

(■NiitY • f>N*LiY«ni6iti*<*f*iino.#YSiore-T*Noic('*PA)t<n>H»»?Yj:pyN i?io 
» VEin*(*«*YMO.»T*NO(«H*«IN({ M-VStOPE)>**2 PYN U20 

5c coNiiKte CYN inc 

A»Y ■ LeNOTHrXhll»NSEOI»Kmi,HSl 6 Al IaYHU,NSECII»VH« I»NStG*l (» 0 ,,PYN 1**0 
I 0.1 OYN mo 

SXY ■ l£NCTH(XP(I*NSeGflWXN(I»N$EGA2WYH J»NSEC*l>»yN<l»NSEG*2)DYN l*fcC 

» ,c..o.i PYN nro 

CXY ’• LtNCTM(XHn»N«fOI»XN(t,NSEC* 2 »»YNH»NSEG|»VNU»NSI 6 * 4 >. 0 .»CYN U #0 

I 0.1 > 2*0 

ANC . A*COS«AXY>EXY»CXY» OYN MCC 

PNALTY • PN 4 UY*EkElTl*«»N*XUO.,«IXY**lON 6 >.>XXYn»** 2 *PWEI| 2 *(*OYN 1310 
I f'«Xl( 0 ..(Z.J! 05 -»*NGM >**2 OYN 1320 

OISUl • 0 IS(ll*l.fNCTHtXT 2 U»NSCG».XHtI,NSEGYl»»YT 2 n»NS£ 6 l»YNMDYN 1330 
1 .NSfOYll»n 2 (I»NSEG). 2 NU»NSEC*ill OYN 134 C 

tc CCMlNUe OYN 1350 

»En,«N PYN 1360 

END OYN 13 T 0 


SLPXOUTINt Nit l*MI) ><11 

CCPfON NTMJ»NhAP,NStG»Xtm5.13)#YNtl5»131 »2N(15fl3 )»SfOSITtl5»25CNtl 
»»3».*0tl5»AXE(151*YCI15l»YF«15).2On51.2f n£».NPOSITS«15),PHOI5T31NII 
t, lStOl» rPROE»M*XIT*P6NC)*lt>*V£IT»TVIOLA.*lUP.PeCPl£.NT»CPT.TNPlON(Nir 
»«731,0R0( 15l.nc.lrpr. ir./EtlOlCB.NEUTPAJ.ABPAYtS'JS.lOl N!I 

INTEGER DRO Nil 

LCflCH IIROLaB.NtWTPAJ Nil 

PEOPIE • 0. Nil 

ALWP ■ C. Nit 

00 10 K • l.NPAR Nil 

IP tARRAY(K.3M.T.l. ) GO TO 1C Nil 

ALCN • THPtONIK) Nil 

If tAlON.Lr.55.1 CC TO 10 Nil 

W ■ 3. 365-6*10, ♦♦I0.103»A10NI Nil 

N ■ V/( 0 » 2 » 10 .*»l 0 . 03 *AtDNI*l.* 3 E-**lC.**<C.Ce*ALDNM Nil 

PEOPLE • PEOPLE *AM<AY(K,3I Nil 

ARSAYIK.5I • APPAY(K» 3 I*W Nil 

ALLP • ALkP.ARPAYIKf 51 Nil 

10 CCNTlNUt Nil 

ANIl . AU,P/P50PU*AWEIT Nil 

PETUEN Nil 

ENP Ml 


10 

20 

30 

AO 

50 

6 C 

7 C 

PO 

00 

IOC 

110 

1 ?C 

130 

lAC 

150 

160 

ITC 

lec 
1 00 
200 
210 
? 2 C 


SUARGL'llNE NOISEl (THRESH) NCI 10 

CQHMOH NTRAJ.NHAP,NSE 0 .XH( 15 . 131 .YN(l 6 .U». 2 P(l 5 . 13 )»SP 0 $lT« 15 r 25 CNrr 2 C 
A. 2 i.'iO( 15 ).XF( 15 ).YC( 15 ).YF<I 5 I.IO( 15 ). 2 F( 15 IfNPasiTSIl.‘).RHO(; 73 )NPI 30 

I, ISECl.IPROFiNAXIT.PENCRlT.AUEl 1 . TVlClA.AlHP.PECPtE.NTRCPT.ThPLON (NPI Ao 
»573 l,aPC( 15 ).lTr.iICPT.n. 2 eRCl 0 B»NeNTRAJ,ARPAY( 573 . 10 ) NCI 50 

INUGEP OPO NTI 60 

LOGICAL IERCLCB.NEWTRAJ NCt 70 

CCCPCN /CCRNER/ ANGl E ( IS. 12 ). PCSl T < 1 5 . 250 , 3 ).XCENTP ( 15.12 I. YCENTR (NCI PO 
AU.iai.TCENTP ( 15 . 12 ),Xn( 15 . 121 .YTl( 15 , 12 I,m{ 15 . 12 ),XTZU 5 . 12 ), YNC 1 90 
»T 2 ( 15 » 12 ),IT 2 ( 15 , 12 ),NP 0 S 1 T( 151 , 0 IS(I 5 ) nOT IOC 

CCNCCN /OYNAHIC/ RAriljS»PkEITl,PHEn 2 »Zk£lTl.CMAPlNM 5 )»CNAPAX( 15 INn 110 
A.NSaHP NCI 120 

CCPCCN /THRESH/ ALCAX.rAXELtT.tVETTl,TkEIT 2 ,NPLANe NPI 130 

CO.HPtN /SANPIE/ SA 6 PT,PtRICiO»SAHlTH,NAC( 15 . 32 ).CNOIS(l 6 , 2 ) NPI 1*0 


X05 


c 

c 

c 

c 

c 

c 

c 


.OI^^NSIOK SfeOlTH<9) 

»tn lENcrn 

cacpo(4^9»c^D) • (A/e»*<c*d)fo 
IENCTH(A/B#C#D#PfP> » 50Rf{(A*B)4*B4(c-0>^^2P(e-f 
PI ■ ArPNii*iP4* 

MAX • I 


IKCUOB AU TRiaS FCP THB MIROTH GLOPAl mRATlON 


IF atPClCB) MAX • MRAJ 
00 KK • WHaX 

I • OPO^XM) 

xTan^NseGAiMXMnf NSEon \ 
xn(UK$IG4n»XT?(l#NS6C*U 
YT2( WN5EG4n»YMn#KSEG4n 
Ynn#N$EC4U-YT2n#NSFG4n 
iTeUf NSEGf U-Zhn #NSEG4n 
nm#NSECHr*ZTEUfNSBCM5 
OLZPM . izon HZF(I) WDIS (I ) 

po$iTu#i#u • xMUjiseou 
P0$ITn»l/2l • YMUiISIGl ) 

PCSIT(l#l#3) • ZM(I#1SEGU 

NSEOl • NSFO-1 

DO 10 J • ISeGX^NSEGl 

SEGLTHUI • tFMGTHniT2n#J>#XTUt>J4l 

I Ui j).zTi<UJnn 

APCUH(J) • RADlLS*ANGUn» J) 

10 CONTINUE 

StClTHiNSEC) • LeNGTH(XT2n#NS6G)#XMtI# 
I NSfGni^ZTZI WNSlGIf ZPO> NSEGn » ) 

II « I 

DC 70 4 • I$EGl#NSEGX 
2C IF IStCLTHNI.iT.SAMlTH) GO TO 30 

II ■ II ♦! 

IF m.0T.Z5C) VP ITE (6^9010) I#JiII 
IF (M.GT*Z50I stop 
PDSIT( l#n>l) ■ CCrROISAMLTHrSECUHU 
I fD) 

POSlTMilIrZ) • CDCPdISAMLTHf SECUHU 
I ^2)) 

^ P0SlTniIt»3) • PCS1TU#II-1»3)-0LZPM 

IF MPRGF.ec.l) PQSITa»n#3I • CCCPO 
% ♦lbPGSn<I#lI-lf3JT 

SeClTH(J) • SEGL1HU)-SAMLTH 
GO TO 20 

30 IF rSEriLTMU)>APClTH(JKGE«SAMtTH) GO 

PlEC^ • samlth-afcuhuj-segltfu) 

n • lui 

IF m.GT.2501 WPITF (6#9CI0) IiJfM 
IF m.GT.aSC) STOP 

PCSITdflWU • CCPF0<P1ECE#FEGITH( J4 
POSIT ( Idl^ZI • CDCPOI PIECE •SEGlTHUf 
POSlTU d da) • FOSlKlf n^lf 3A-01.ZPM 
IF (IPFOF.EC.U PQSmi#n#3i • COQPO 
1 J42)#ZT2nt J41J) 

SEGUMdtl) • SEOtTN(>ltlji-PlfCE 
CO TO 70 

AC PIECE • SAHlTH-SFGLtMI J) 

GAMMA « flECE/PAOMS 
ARCUHUI ■ APCITHUI-PIECF 


NPI 

NCI 

NCI 

NCI 

NCI 

NCI 

NCI 

NCI 

NOI 

NCI 

Nfl 

NCI 

NOT 

NCI 

NCI 

NCI 

NCI 

NPI 

NCI 

NOI 

NCT 

NPI 

NOI 

NCI 

NCI 

NCI 

NOT 

NOT 


)fYT2(t^4I,YTI riiJ4| bZTZNOT 

NPI 

NOI 

NCI 

NSEG4U# YT2f If NSEGI# YMU^NOI 

NOI 

NCI 

NOT 

NPI 

NOI 

NPT 

NOI 

i.xni If jfUf pcsiTdf ii-iNor 

NCT 

)rm Uf 441 ifPosrmfii-iNoi 

NOT 

•SAMLTH NCT 

(SAMLTHf SEGLTH<4)»Zmii 4NPT 

NPT 

NOI 

NO! 

TO AC NCT 

NCI 

NOT 

NCI 

NOT 

n#XTI(If 442)f XT2(lf J41) )NPI 

UfYTin»J42)fYT2Mf J4n )NP| 

•SAMITH NOT 

(PIECEf SEGlTHUtUimn »N01 

NCI 

NPI 

NCI 

NOI 

NCT 

NPI 


170 
lAC 
l«0 
200 
21 C 
220 
230 
2*0 
250 
200 
270 
220 
290 
300 
310 
32C 
33C 
3*0 
350 
3frO 
370 
300 
390 
AOO 
AlO 
A2C 
430 
*40 
450 
4ff 
470 
400 
400 
50C 
510 
520 
53C 
5*0 
55 C 
550 
570 
5 00 
50p 
600 
610 
620 
630 
640 
650 
66C 
670 
60C 
69C 
700 
710 
72C 
730 
740 
750 
76C 
T7C 
7A0 


186 



M • ll*l NOI 
IF m.CT.2$CI VeiTE IE«<iO)C) hJtU HOI 
If (II.GT.25CI STOP NOI 
UfPI . 4r*S2MyTin»J*I>-YCENTMI.J)l#imn.Jfll-KCCNTi<I* JiNOI 
M HOI 
IPtUS • 1 HCI 
At ■ ALfAt*AKCLt(l>J)/A. NCI 
K22 ■ XCENTP(I*JKPA01US*CCS(AAI HOI 
T2Z • yC£HTPM»JMP*OIUS»5IN|A*l NCI 
AOIS • S0RTIIX22-XHIIfJ*in**Z*|yZZ-I«n»J*I!l**2> NOI 

effis • scRiiixTii-i# jfi»-*Nn»j*in**2#(Ynn» j*i>-YN(i» j*iM**Noi 

21 NCI 
IF IBOIS.LT.AOISI IPltS • -1 HCI 
AA • AlFAl*IPUS»ftAP«A NCI 
POSITM.INn • xaMSII.J|*l<AOIU$»COS(**l NCI 
P0SII(I»II<2) • YCENTPI If J)ARAOIUS*SIHIAAI NCI 
FOSim»II» J» • POSITIIrII-l» 3l-OUPN*SAPlTH NCI 
IF MPPOF.EO.U PCSITM.IW3I • COOPO IP lECE, APCLTHI Jl ♦PIECI»2 INC! 
2II>J*U»ZTU WJ*1I) NOT 
IF lAPClIHI JI.U.SAHITMI 00 TO tC NOI 
AA • AAtIPll,S»SAPlTH/RAOIOS NOI 
APCITHIJI • APCUHIJI-SAHLTH NCI 
II • IIAI . HCI 
IF III,GT.2S0t kKITE 16>90t0) IfJttl NCI 
IF m.CT.ZiOI STOP NCI 
PDSITIIrtWlI • XCENTPiIHI«P*OIUS*CaSIAA) NCI 
PC$im»IIf2l • YCENTPI I,JI*RAOIUS*SINIAAI NCI 
POSITI I»II»3( • POSIT(I»II-1»3I-OL2PN*SANLTH NOI 


IF (IPPOF.SC.n P0SIT(I»IW3> • COORO(SAnTH»ARCLTHUMSANLTH»Nri 
2T2(I»JNl)»PCSlI(WII=l»3l> HOI 
GO TO 50 NCI 
PIECE • SAPIIH-ARCLTHIJI NCI 
II • II»1 NCI 
If (II.6T.2iCI VRITE (O^IOIOI NCI 
IF III. 01. 2501 STOP NCI 
POSITI I. II.2) • CCCP0<PIECE/SGCLrH(Jtn.m(I»JP2I.YT2(I«JPl) )NCI 
POSITII.II.n • CCCBO(PlECE.S£OLTH(J«l).XHn.JA2)>XT2(IiJ*lMNOI 


P05IIIl»n.3l • POSITn»H-l»3l-0LZPN«SAHTH NOT 

IF IIPBOF.EC.ll P0SIT(l»II»3l • COORD I PIECE» SEOITHIJP I > ,211 1 1 »NOI 
I J*21. 212(1, Jtin NCI 

SEG(.TH(J«l) • SEGITN(J»II-PIECE NCI 

CONTINUE NCI 

IF (SEOLTHINSECI .IT.SANLTHt GO TC 90 Hfl 

It • II ♦! NCI 

IF (II. GI. 2501 WRITE (fc,9CI0» I,J,II NCI 

IF (II. 01. 250) STOP NCI 

PCSIT( I,II,ll • CCCPO(SANUH,SEOl IH(NSeG) ,XTI( I,NSEGNII,P0SIT(I,NCI 
I INI, 1)1 NCI 

POSITI 1,11,21 • CCCRO(SA(lUH,SEGLTH(NSEG|,m(I,NSEG,l),PCSTT(I,NCI 
I 11-1,2)) HCI 

M0SirtI,M,3) ■ P0SIT(I,II-1,3)-CIZPN»SAPLT(( NCI 

IF (IPOOF.EQ.I) PCSMII, 11,3) • CCCRD (S ANITH, SEGITHI NSEC) , Ml ( I ,NCI 
I NS£C*1),PDSIT(I, 11-1,3)) NOI 

SEGITHINSEO) • StCLlHINSEO-SAHTH NCI 

CC TO EO NOI 

Nposmi) » n NOT 

IF 111. GT. 250 CAU. ERROR (8, "NO. Of TPAJ. SAHpLES .GT. 250, NCI 
I Se-OIHENSICN AFPAYS posit ANO SPOSIT") NCI 

rCNTINUE NCI 

NSANP . JFIX(PERlCD*3frOO/SANPT) NCI 

I viola • 0. NCI 

DO 110 K • 1,NNAP NCI 

AMEN .0. NCI 

NCI 



c*t. eNfJicy cohtMB. by AttaiiNT nqi$e» 

CUVRENTLY *HStENT LtVEl • 0« OB 


ii!f;AY(K»«) ■ Ms#r‘r*io.**iAr*BK-yio.) 

CCNTtNUE 


IF THIS IS THE IfHCTH GLOB*l ITERBTION# IHCIHOE 
NOISE EFFECTS CF *U TS/kJ$ IN TFPLON 
EISE» C*IC. EFFECT OF TF*J. BEING OPT. AND 
STCFE IN *PB*YI*»*) 


IF (.nct.zeshcb.op.hff.eo.i) c-c to ItO NOI 

CC 150 F • l.NHFP NCI 

8LKVI0L • 0, NCI 

IF liRR4Y<K,3),lT.l) CO TO 1*0 NCI 

00 130 II • 2»F»> NOT 

I • OROMII NCI 

NPCsiTi • Nposnm nci 

ENEFCY • 0. NCI 

CO 120 J • 1. NPCSITI NOI 

fUNC-E • SOFT (IPOS IT < U J>1 1»*aP*Y(F » 1 1 )*»2* (POSIT (I» J» 2V-FFF»N0l 
t YIK,2n**2*(PCSIT(I»J»3H**2I Nfl 

U • CN0Isn.il-CNCIS(If2)PH0610(P*NeE) NCI 

ENEFCY • £NEBCY»10.»*(*L/10.» NOI 

12C CCMINUc NCI 

TNPLCN(KI • ENEFCY»TNPtON(K) NCI 

130 CONTINUE NCI 

1*C TPPLDNIKI • 10.»FlCGIOnNPlON(K|/(PEPIOO«3600./S*PPTI) NOT 

eiKVICL • iN*Xl(C.»BlKVial-*lN*X»P**FUn NCI 

IF (lPR*Y(K.10I.t0.1.) TVICU • TVlCLFfBLKVICL NOT 

15C CONTINUE NCr 

Itc OC leo K • I.NMlP NCI 

IF <*PF*Y(K.3).IT,1) GO TO IfO NOT 

1 • QFCdl NCT 

NPCSITI • NPOSIIII) NCI 

tNERGY • 0. NOI 

OC 17C J • WNPOSITI NT! 

F*NCF • SOPT(IPCSIT(Id»U“*RP*Y(K>l) »*»2*(P0SIT(I. J, 2 I-*FP*Y(NCI 
$ R.2M •*2*<PCSlin» J.3II ♦•21 NOI 

AL • CNOIS(l»lI-'CNOISn»2)»ALOG10(RANGE) N(?! 

ENERGY • ENEFGY*10.*A(AL/IC.I NOI 

170 CONTINUE NCI 

APPAY(K.A) ■ IC4*AIOG10(ENERGY/(PEPIOO*3NOO./SANPTM NCI 

leC CONTINUE NTI 

NC* 

NCI 


AC(' THE NCISE FPCr TFAJ. EEING CPT. TC THE OTHER NOISE 


IF (zebglc0.cR..not.neiiTR*j) call ACOR (II 
THRESH • IVEII1»TV1CLA**2 
NEVTRAJ • .FALSE. 

RETURN 



90X0 FC»f«*T I"! • J • II • 

ENP 


NOT 

NOt 


U’BPOtnNE pcvpy ( 

PfV 

STOP 

PCV 

ENO 

PCV 


SU99CUTINE DNOrOEN IHXXITC^ STCPCHC* N» XNOW* DEI TXX> NDECI «N0 

CCHPCN NTPXJ#NH«E>^SEC»XN41Srl}l»YNaS>l^)rXNIl$»nl»SPOStmS*^9C9ND 
«O)>XCU;i»XFIlSI*t0(15)#YFIlSI>Za(X!)>ZF(l!l«NP0SITSII9)»l)H0(E73»9X0 
t, ISEPl»lF>OF»N*XlT»F£NC(iIT»*WUT. TVIOU»H¥F»F£OFl£»NTFCFT*T«FLONIFNO 
fST3l>090US)«ITC->ICPT*IT»ZEFClce*NEhTF«0>«RPAy(9T3»XOl FNO 

1 NIECE# OFO • FNO 

LOCICFL ZE#GL0e*NEk1##J «N0 

OINENSION XN0W(NCEC»1 DELUXtNDECi #N0 

NSS6X • NSEC-l »N0 

DO 30 KK • l.NAXITC NNO 

I TO • KK * FNC 

DC 30 I • 1»NI#0FT »N0 

lOFT • I 

N£>TF»J . ,T»te. #N0 

b#ITE I6>901CI IO#D( JI>0«1>NT##0 I PNO 

CALL ONEUrON IN.XKCbtDELI AX<S TOFCHCt NDEC ) PND 

C ■ ' PNO 

C PNO 

C , . PNO 

C . TNPLON NOW CONTAINS NOISE EFFECTS FPCP THE CONSTANT TPAJS# . PNO 

C . SC ADO IN THE NOISE FPON THE TPAJ. 8EING OPT. . PNO 

C . . PNO 

C . PNO 

C PNO 

CALL ACCR (II PNO 

IF (NTPCPT.EO.I) FETL.PN PNO 

C PNT 

C PND 

C . .PNC 

C . CYCLE THE ORDER OF THE TRAJS . PNO 

C . .PNO 

C ..PNO 

C K N D 

Nl • NTPCPT-I 

CfCTf^P • CPOa» PND 

OC IC • l#M PKD 

OPO(KI • OROtK>l) PMD 

IC CDNTI^UP PK’P 

OPD^NTRCPT) ■ CROTHP P^’O 

00 20 J • l#NSEGl PND 

CO 20 K ■ li3 PNO 

L • ( PNO 
If (KrEOtll yNCWttfi) • Xn(CROm#J#l) P^‘D 

If iK.ec.2) i^NcuMri) • YH<cpo(;)#j*n pno 

If <K.£C.3) }>NCV(l#U • Zf<fOROm#J>n P^o 

2C CGKTlNie PNO 

3C CCNtlKUf PKO 

WPITc (6»<9C2CI NIPCPr#f<*KUG Pf*D 

RETUPN PKO 

C PND 

«CIC fCfW/4 (/>5X#»*CPDEP Cf TPAJSt •*»15M2f2YM PNO 


2070 
20 PC 


1C 

20 

90 


10 
20 
30 
^0 
90 
60 
70 
PO 
90 
ICC 
lie 
120 
190 
19C 
190 
160 
17C 
IPO 
190 
200 
210 
220 
220 
290 
290 
26C 
27C 
2P0 
29C 
3CC 
310 
320 
33C 
39C 
350 
360 
370 
39C 
300 
900 
910 
42C 
930 
990 
650 
460 
970 
9 PC 
9 00 
900 
9 1C 


9QZ0 (/*5X,«RCUNC fCBJK PROCeSS CCHPUTI. E*CH OF THf tPPNF 5iC 

Mjs cpTiPi/eo "iiz#” sso 

END PNC 9*0 


SLSPOtTJHE *PBE nSIflM *C0 

CCPCN NT**i/*Nft*.-<FSCG»»mU»J3}»VFnS*13»»lFMJ»13»»5POSlT(19»i5C*DC 
I. J)» *0M5).S!>(l9V,YGMSJ*yFU5»tEC(15»*eFil9r»NPOSnsm>»PHO(J73)*0O 
l*ISI 61 »lPPPFiM*XIT»PeHCjlITr*Hen»TVra«.*»*lHP»P£CPli»NTPOPT>THPlONODP 
»S73)*0P0U5»»ITG,ICPT»IT»IIPClCefH£WTP«J»**F*ve973>lOI *00 

IKTECEP 0»0 *PP 

LCGIC*!. IERCLQB»N£WTP*J ACO 

C*U SECOND (TO) AC(' 

PO 10 K » 1,NN*P ACO 

£l • 10.*»t*RR*7(K»*)/lC,) *00 

E2 • lO^**(TNPlCN(KmQ.) * ado 

ENEN • fc2*JSIfiN*£l ADD 

TPPLfN(K) • 10.**t,CC10(EN£W) *00 

3 CONTINUE APR 

SETUPN ACD 

EM) *00 
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